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ABSTRACT 


The  electrical  slow  wave  activity  of  the  colon  was  studied  using 
in  vitro  muscle  strips  of  canine  colon.  Recordings  were  obtained 
using  both  extracellular  pressure  electrodes  and  intracellular  micro¬ 
electrodes.  Electrical  coupling  between  regions  was  studied  using 
multisite  measurements  on  each  side  of  the  intact  muscle  wall  and  on 
each  side  of  isolated  circular  and  longitudinal  muscle  layers.  Coupling 
and  intrinsic  oscillator  frequencies  were  also  studied  by  calculating 
frequencies  and  phase  angles  before  and  after  cutting  the  tissue 
between  recording  sites.  The  origin  of  the  slow  wave  and  the  roles 
played  by  the  individual  muscle  layers  in  the  generation  of  the  slow 
wave  were  studied  using  the  simultaneous  application  of  extracellular 
electrodes  to  both  sides  of  the  intact  muscle  wall  and  isolated 
muscle  layers.  The  above  points  were  also  examined  using  intracellular 
microelectrode  measurements.  Also,  the  i ntracel 1 ular-to-extracel 1 ul ar 
transduction  mechanism  was  examined  using  simultaneous  intracellular 
and  extracellular  measurements  and  intracellular  measurements  in  cells 
directly  below  simultated  pressure  electrodes.  Computer  programs 
were  written  to  manipulate,  process  and  display  experimental  results. 

The  results  may  be  summarized  as  follows: 

1.  The  electrical  slow  wave  originates  at  the  boundary  of  the  inner 
circular  layer  and  the  submucosa,  and  the  integrity  of  the  circular 
muscl e-to-submucosa  interface  is  essential  for  the  existence  of 
the  slow  wave. 

2.  The  slow  wave  is  not  generated  by  cells  of  the  longitudinal  muscle 
layer. 
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3.  The  mean  resting  potential  and  slow  wave  amplitude  of  cells  in 
the  circular  muscle  layer  was  -77.8  ±  6.8  mV  and  28.9  ±  5.5  mV 
respectively. 

4.  The  mean  resting  potential  of  longitudinal  muscle  cells  was  -52.6 
±  7.0  mV. 

5.  With  regard  to  electrical  coupling  between  oscillatory  regions  in 
the  circular  and  longitudinal  directions: 

(a)  Coupling  was  better  in  the  circular  direction  than  in  the 
longitudinal  direction.  The  mean  phase  angle  for  an  elec¬ 
trode  spacing  of  2.5  mm  in  the  circular  direction  was  29.0  ± 
23.2  as  compared  to  117.3  ±  51.0  degrees  in  the  longitudinal 
direction. 

(b)  No  gradient  of  intrinsic  oscillator  frequencies  existed  in 
the  longitudinal  direction  but  a  frequency  gradient  was 
present  in  the  circular  direction  with  respect  to  the 
mesenteric  border. 

6.  With  respect  to  the  intracel 1 ul ar-to-extracel 1 ular  transduction 
mechanism: 

(a)  The  results  do  not  support  the  generally  accepted  pressure 
electrode  hypothesis  proposed  by  Bortoff  (1967). 

(b)  Contrary  to  this  theory  the  application  of  a  pressure  elec¬ 
trode  resulted  in  a  reduction  in  the  intracellular  potential 
of  the  slow  wave  trough  but  did  not  change  the  intracel 1 ul ar 
potential  of  the  slow  wave  plateau. 
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CHAPTER  1 


INTRODUCTION 


1.1  MOTIVATION 

Electrical  signals  are  known  to  play  a  major  role  in  the  control 
of  many  physiological  systems.  The  majority  of  these  signals  may  be 
placed  into  one  of  two  categories.  Either  they  are  of  a  pulsatile 
nature  such  as  the  action  potentials  which  permit  communication  within 
the  nervous  system  or  they  are  cyclic  oscillatory  signals  of  the  type 
generated  by  the  sinoatrial  node  of  the  heart. 

It  has  been  known  since  1922  (Alverez  and  Mahoney)  that  the 
gastrointestinal  tract  is  electrically  active.  It  has  been  shown 
since  that  time  that  both  action  potentials  and  slow  oscillatory  type 
signals  generated  within  the  organs  themselves  play  a  role  in  the 
control  of  contractions  in  the  stomach  (Sarna,  1971),  small  intestine 
(Sarna,  1971),  and  colon  (Shearin  et.  al . ,  1979).  These  electrical 
signals  and  their  control  function  are  best  understood  in  the  stomach 
and  least  understood  in  the  colon. 

The  purpose  of  this  work  is  to  increase  the  understanding  of  the 
electrical  activity  and  its  control  function  in  the  human  colon.  The 
immediate  incentive  for  persuing  this  is  the  acquisition  of  informa¬ 
tion  which  will  lead  to  improved  clinical  treatment  of  colonic  dis¬ 
function. 
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1.2  SCOPE  OF  THE  PRESENT  WORK 

The  study  of  the  electrical  activity  of  the  colon  could  be  pursued 
in  a  number  of  different  ways.  One  could  study  the  intrinsic  cyclic 
electrical  activity  generated  by  the  smooth  muscle  of  the  colon  wall, 
the  neural  control  of  the  organ,  or  the  integration  of  both  the  intrinsic 
and  neural  signals  in  the  control  of  colon  function.  The  first  of 
these  approaches  was  adopted  by  this  present  work. 

The  intrinsic  cyclic  electrical  activity  of  the  gastrointestinal 
tract  has  been  referred  to  as  slow  waves,  electrical  control  potentials, 
and  basic  electrical  rhythm  (Bass,  1968;  Connor,  1979;  Szurszewski, 

1981).  The  term  slow  waves  will  be  used  in  this  thesis.  The  principal 
objectives  of  this  work  are  (1)  the  definition  of  the  origin  of  the 
slow  wave  activity  in  the  colon,  (2)  the  definition  of  the  role  played 
by  each  of  the  muscle  layers  in  the  generation  of  the  slow  wave,  and 
(3)  the  determination  of  the  basic  electrical  characteristics  of  this 
signal.  The  last  of  these  will  require  an  examination  of  the  intra¬ 
cellular  resting  potentials  and  slow  wave  amplitudes  as  well  as  an 
examination  of  the  phase  relationships  of  signals  recorded  from  different 
si tes . 

A  number  of  alternative  ways  of  attaining  these  objectives  were 
considered.  One  method  was  to  study  the  electrical  activity  in  vivo 
in  humans  using  surgically  implanted  electrodes.  An  alternative  to 
this  was  the  study  of  the  activity  in  vitro  from  organ  specimens 
removed  during  required  surgery.  Another  approach  was  the  study  of 
tne  colonic  electrical  activity  in  animals  using  either  in  vivo  or 
in  vitro  recording  techniques. 
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The  first  of  the  above  methods  was  seriously  considered  because 
it  would  examine  the  electrical  activity  of  the  human  colon  under 
close-to-normal  physiological  conditions.  However,  a  number  of 
difficulties  existed  with  this  approach.  The  first  of  these  difficulti 
was  in  finding  a  sufficient  number  of  subjects  for  the  study.  Elec¬ 
trodes  could  only  be  placed  in  patients  who  required  abdominal  surgery 
and  who  gave  informed  consent.  The  availability  of  such  subjects 
would  require  a  long  period  of  time  to  obtain  sufficient  data.  Since 
this  study  would  be  on  post  operative  patients,  many  of  whom  were 
being  treated  for  some  colon  disfunction,  the  normality  of  the  data 
might  be  suspect.  In  addition  to  this,  currently  available  in  vivo 
electrodes  produce  very  poor  quality  signals.  These  electrodes 
usually  consist  of  stainless  steel  or  platinum  wires  inserted  into  the 
colon  muscle  wall.  The  signals  obtained  from  such  electrodes  are  of 
very  low  amplitude  and  contain  large  noise  components  due  to  electrode 
movement. 

The  second  approach  considered  was  the  study  of  human  colonic 
muscle  in  vitro.  This  approach  has  advantages  in  comparison  to  in 
vivo  studies  because  it  does  not  have  to  contend  with  electrode  move¬ 
ment  artifacts  and  because  signals  obtained  from  in  vitro  electrodes 
are  usually  larger  in  amplitude  than  in  vivo  signals.  However,  this 
type  of  study  is  hindered  by  the  lack  of  availability  of  healthy 
human  tissue.  Specimens  are  only  available  from  patients  who  require 
colonic  surgery.  Therefore,  these  specimens  are  obtained  from 
abnormal  organs.  Because  the  surgeon's  primary  concern  is  for  his 
patient,  and  not  for  the  health  of  the  specimen,  the  viability  of 
such  specimens  would  be  questioned. 
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The  third  method  of  study  considered  was  the  in  vitro  and  in  vivo 
study  of  slow  wave  activity  in  an  animal  with  colon  and  colonic  functio 
similar  to  man.  The  major  advantage  of  this  approach  was  the  adequate 
supply  of  animals  and  animal  tissue.  However,  animal  in  vivo  studies 
have  electrode  and  movement  artifact  problems  similar  to  their  human 
counterparts.  In  vitro  animal  studies  have  the  disadvantages  (1)  that 
the  tissue  is  not  maintained  in  exactly  the  same  physiological  state 
as  it  would  be  in  vivo  and  (2)  that  the  tissue  under  study  is  not 
human  tissue.  Species  differences  could  make  the  results  inapplicable 
to  the  human  organ.  However,  studies  of  the  electrical  activity  of 
the  canine  stomach  and  small  intestine  have  greatly  increased  the 
understanding  of  the  electrical  activity  in  the  corresponding  human 
organs  (Sarna,  1971;  El-Sharkawy  et.  al . ,  1978;  Szurszewski,  1931). 

In  vitro  studies  have  the  advantages  of  (1)  a  ready  supply  of 
tissue,  (2)  few  movement  artifacts,  (3)  larger  signals  from  electrodes 
which  are  better  understood  than  in  vivo  electrodes,  and  (4)  the 
facility  to  permit  changes  in  the  ionic  composition  and  the  application 
of  hormonal  agents  at  the  measurement  site.  For  the  above  reasons, 
and  because  in  vitro  canine  studies  have  proven  very  useful  in  many 
other  gastrointestinal  studies,  the  research  described  in  this  thesis 
is  based  entirely  on  studies  of  in  vitro  strips  of  canine  colon. 

The  work  uses  both  extracellular  sil ver-sil verchloride  pressure  elec¬ 
trodes  (Kingma,  et.  al.,  1982)  and  intracel 1 ul ar  microel ectrodes . 
Because  the  relationship  between  extracellular  pressure  electrode 
signals  and  their  associated  intracel  1  ular  potential  variations  has 
not  been  adequately  defined,  it  was  necessary  to  investigate  the 
pressure  electrode  transduction  mechanism.  Also,  because  the 
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characteristics  of  sil ver-si 1 verchl oride  electrodes  were  not  known 
for  the  very  low  frequencies  (4  to  7  cycles/minute)  found  in  the  colon, 
it  was  necessary  to  measure  the  electrode  impedance  versus  frequency 
characteri sti cs . 

The  extracellular  studies  included  multisite  measurements  on 
strips  of  both  the  intact  colonic  muscle  wall  and  of  isolated  muscle 
layers.  Also,  they  included  simultaneous  measurements  on  opposite 
sides  of  the  intact  muscle  wall  as  well  as  on  opposite  sides  of  isolated 
muscle  layers. 

The  microelectrode  studies  involved  the  recording  of  intra¬ 
cellular  potentials  from  cells  in  both  muscle  layers  of  the  colonic 
wall.  These  measurements  were  also  used  to  compare  extracellular  and 
intracellular  signals  and  to  investigate  the  mechanism  responsible 
for  the  generation  of  the  extracellular  pressure  electrode  signal. 

1.3  THESIS  ORGANIZATION 

Because  of  the  interdisci  pi  inary  nature  of  this  work  it  is  expected 
that  it  will  be  read  by  both  engineers  and  medical  scientists.  To 
make  the  work  more  easily  understood  by  the  engineering  reader  Chapter 
2  provides  a  description  of  the  anatomy  of  the  colon.  Chapter  3 
provides  a  review  of  studies  of  electrical  slow  wave  activity  con¬ 
ducted  on  the  stomach,  small  intestine,  and  colon.  This  chapter 
also  includes  a  description  of  the  currently  accepted  theory  of 
pressure  electrode  signal  generation  and  of  slow  wave  modelling  using 
relaxation  oscillators.  Chapter  4  describes  the  research  objectives 
as  well  as  providing  a  detailed  description  of  the  experiments  conducted 
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to  attain  these  objectives.  Chapter  5  provides  the  results  of  all 
experiments  and  Chapter  6  includes  a  discussion  of  the  results  and  of 
the  conclusions  which  can  be  deduced  from  these  results. 


CHAPTER  2 


ANATOMY  OF  THE  COLON 


2.1  INTRODUCTION 

This  chapter  is  included  to  provide  the  reader  with  a  basic 
knowledge  of  the  anatomy  and  function  of  the  colon.  This  will 
facilitate  a  better  understanding  of  the  motives  and  significance 
of  the  research  described  in  this  thesis.  It  will  also  define 
physiological  terms  using  in  the  following  chapters. 

Much  of  the  data  provided  is  from  animals  other  than  man,  but 
where  possible  data  from  the  human  colon  is  included.  Since  the 
canine  colon  is  to  be  used  in  this  work  the  differences  and  similarities 
of  human  and  canine  colon  are  particularly  noted. 

2.2  GROSS  ANATOMY  AND  FUNCTION 

The  colon  or  large  intestine  is  the  portion  of  the  gastrointes¬ 
tinal  tract  extending  from  the  ileocecal  value  to  the  rectum.  The 
functions  of  the  colon  are  (1)  storage  of  fecal  matter  and  (2) 
absorption  of  water  and  el ectrolytes .  The  proximal  half  of  the 
colon  is  concerned  principally  with  absorption  and  the  distal  half 
with  storage.  The  caliber  of  the  organ  is  greatest  at  the  proximal 
end  and  its  diameter  reduces  towards  the  rectum.  As  shown  in  Figure 
2.1  there  is  an  alteration  of  fixed  and  mobile  parts  due  to  variations 
in  the  degree  of  fusion  to  the  dorsal  body  wall  by  means  of  a  membranous 
fold  called  the  mesentery  (Schofield,  1968). 

A  diagrammatic  cross  section  of  the  alimentary  canal  illustrating 
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structures  typical  of  the  canine  colon  is  shown  in  Figure  2.2.  As 
can  be  seen  from  this  figure  the  colon  wall  is  composed  of  a  number 
of  distinct  layers.  The  outer  lining  which  is  continuous  with  the 
mesentery  is  called  the  serosa.  The  muscle  wall  shown  immediately 
inside  the  serosa  consists  of  two  smooth  muscle  layers,  one  oriented 
in  the  longitudinal  direction  and  the  other  in  the  circular  direction. 
The  muscle  wall  facilitates  movement  and  mixing  of  bowel  contents  to 
improve  absorption  of  fluids.  Also  the  muscle  coat  is  responsible 
for  the  mass  movement  of  wastes  to  the  distal  end  of  the  colon  and  for 
final  evacuation  of  these  wastes.  As  seen  in  Figure  2.2  the  longitu¬ 
dinal  muscle  coat  of  the  canine  colon  is  continuous  around  the  complete 
circumference  of  the  organ.  However,  in  the  human  the  longitudinal 
muscle  coat  is  gathered  mainly  into  distinctive  bands  or  taenia  coli. 
The  mucosa  or  inner  lining  of  the  colon  is  responsible  for  the  absorp¬ 
tion  of  fluids  and  el ectrolytes .  Associated  with  the  mucosa  is  another 
thin  smooth  muscle  layer,  the  muscularis  mucosae.  This  muscle  layer 
causes  movement  of  the  mucosa  to  maximize  fluid  absorption.  There 
is  a  high  density  of  nerves  in  the  colon  wall.  The  majority  of  these 
nerves  is  found  in  two  nerve  networks  or  plexuses.  Auerbach's  plexus 
is  a  network  of  nerves  and  ganglia  located  between  the  longitudinal 
and  circular  muscle  layers.  Meissner's  plexus  is  located  in  the 
submocosa  which  is  a  connective  tissue  layer  between  the  mucosa  and 
the  circular  muscle  layer.  These  nerve  networks  form  a  sophisticated 
system  controlling  colon  function.  This  control  system  combines 
both  local  reflex  loops  in  the  colon  wall  and  extrinsic  inputs  from 
the  spinal  cord  and  higher  centers  to  control  absorption,  muscle 
contraction  and  mass  movement. 
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1  Attachment  of  the  human  colon  to  the  dorsal  wall 
via  the  mesentery.  (From  Schofield,  1968.) 
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2  Diagrammatic  cross  section  showing  the  structure  of 
the  alimentary  canal.  (From  Schofield,  1968.) 
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2.3  THE  MUSCLE  COAT 

There  is  much  species  variation  in  the  structure  of  the  smooth 
muscle  layers.  However,  all  species  studied  have  both  longitudinal 
and  circular  muscle  layers  composed  of  long  thin  spindle-like  cells 
ranging  in  diameter  from  4  to  6  pm  and  in  length  from  200  to  400  ym 
(Bennett  &  Rogers,  1967;  McGeachie,  1975;  Prosser  et.  al . ,  1960;  and 
Schofield,  1968).  A  connective  tissue  sheath,  the  epimysium,  completely 
surrounds  the  smooth  muscle  of  each  layer.  Thin  septa  extend  inward 
from  the  epimysium  to  form  the  perimysium  which  contains  collagenous 
and  elastic  fibres  as  well  as  fibroblasts,  capillaries,  and  nerves. 

An  example  of  this  arrangement  in  the  longitudinal  muscle  of  the  dog 
intestine  is  shown  in  Figure  2.3.  The  perimysium  divides  the  muscle 
into  bundles  of  smooth  muscle  cells.  This  formation  of  descrete  bundles 
has  been  observed  in  the  guinea-pig  and  cat  (Prosser  et.  al  . ,  1  960), 
the  dog  (Cajal ,  1933),  and  the  human  (Schofield,  1968).  Transverse 
sections  through  the  muscle  bundles  show  that  they  vary  in  outline 
from  rectangular  to  circular  and  are  20  to  200  ym  wide.  A  cross  section 
of  muscle  bundles  in  the  circular  muscle  layer  of  human  colon  is  shown 
in  Figure  2.4.  The  muscle  bundles  can  be  traced  in  serial  transverse 
sections  through  a  smooth  muscle,  before  they  lose  their  identity  by 
anastomosing  with  adjacent  bundles.  Because  it  is  difficult  to.  define 
where  one  bundle  ends  and  the  next  begins,  no  concensus  exists  on  the 
length  of  these  bundles. 

The  packing  of  smooth  muscle  cells  into  a  bundle  is  shown  in  cross 
section  in  Figure  2.5.  Each  cell  in  the  bundle  has  about  five  close 
neighbours.  These  cells  are  separated  by  a  protein-polysaccharide- 
filled  gap  of  50  to  80  nm.  The  three-dimensional  relationship  between 
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Figure  2.3  Diagram  of  a  longitudinal  section  through  the  smooth 
muscle  of  the  dog  intestine;  the  muscle  bundles  (b) 
and  the  perimysium  (P)  are  clearly  shown;  the 
perimysium  consists  of  collagenous  and  elastic  fibres 
as  well  as  capillaries  and  nerves;  at  irregular  inter¬ 
vals  adjacent  bundles  fuse  or  small  bundles  (I)  inter¬ 
connect  larger  bundles.  (From  Cajal ,  1933.) 
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Figure  2.4  Cross  section  of  the  human  colon  wall  showing  the 
smooth  muscle  bundles  in  the  circular  muscle  layer. 
(From  Schofield,  1968.) 
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Figure  2.5  The  arrangement  of  smooth  muscle  cells  in  a  muscle 
bundle  of  the  guinea-pig  vas  deferens.  Electron 
micrograph  montage  of  a  very  small  muscle  bundle 
(Bo)  in  the  longitudinal  muscle  coat.  Note  that 
the  muscle  bundle  is  clearly  separated  by  a  layer 
of  connective  tissue  (o) ,  the  perimysium,  from  the 
neighbouring  muscle  bundles  B1  and  B2.  Each  muscle 
cell  (M)  at  the  centre  of  the  bundle  is  surrounded 
by  about  five  other  muscle  cells;  a,  axons;  5, 
Schwann  cell.  (From  Bennett,  1972.) 
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muscle  cells  in  a  bundle  has  been  determined  from  studies  of  the  guinea- 
pig  (Bennett  &  Rogers,  1967)  and  the  mouse  (Taxi,  1965).  In  these 
animals  each  cell  is  closely  apposed  by  approximately  twelve  cells 
at  any  point  along  its  length.  Approximately  equal  numbers  of  cells 
overlap  at  either  end  of  any  given  cell. 

The  muscle  cells  within  the  bundles  are  connected  to  each  other 
at  points  along  their  lengths  by  three  different  types  of  structures 
called  tight  junctions,  desmosomes,  and  gap  junctions  (Staehelin  & 

Hull,  1978).  The  main  function  of  tight  junctions  is  to  enable  the 
group  of  cells  to  maintain  an  internal  environment  that  is  different 
from  the  external  one.  Desmosomes  are  the  mechanical  links  which 
hold  the  tissue  together.  Gap  junctions  are  believed  to  play  a  role 
in  cellular  chemical  and  electrical  communication.  A  broad  correla¬ 
tion  exists  between  the  appearance  of  gap  junctions  and  the  transmission 
of  action  potentials  with  a  tissue  (Dewey  &  Barr,  1962).  A  model  of 
a  gap  junction  is  shown  in  Figure  2.6.  These  junctions  consist  of  a 
2  nm  gap  between  the  apposed  outer  lamella  of  adjoining  cell  membranes 
which  is  bridged  by  an  hexagonal  array  of  structures  joining  the  mem¬ 
branes  together.  It  is  probable  that  these  structures  which  bridge 
the  2  nm  gap  provide  a  low  resistance  pathway  between  cells  (Brightman 
&  Reese,  1969). 

2.4  CONTROL  OF  CONTRACTIONS 

Muscle  contraction  is  generally  controlled  by  nerves.  In  skeletal 
muscle  the  effector  unit,  the  smallest  unit  controlled  by  the  nervous 
system,  is  the  single  muscle  fibre.  A  propagating  wave  of  depolariza¬ 
tion  (an  action  potential)  arrives  at  the  end  of  a  nerve  causing  the 
release  of  a  chemical  transmitter  substance.  For  skeletal  muscle 
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Figure  2.6  Model  of  a  Gap  Junction  depicts  the  structural  elements  that 
allow  the  exchange  of  nutrients  and  signal  molecules  between 
cells  without  loss  of  material  into  the  intercellular  space. 

The  communicating  "pipes'1  are  formed  by  pairs  of  abutting 
particles,  which  are  in  turn  composed  of  six  dumbbell -shaped 
protein  subunits  that  span  the  lipid  bi layer  of  each  cell  mem¬ 
brane.  The  channel  passing  through  the  cylindrical  particles 
is  about  20  angstroms  in  diameter,  limiting  the  size  of  the 
molecules  that  can  pass  through  it.  Unlike  the  tight  junction, 
fluids  and  tracers  in  the  intercellular  space  can  permeate 
the  gap  junction:  they  flow  around  the  pipes.  (From 
Staehelin  &  Hull,  1578.) 
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the  nerve  ending  is  in  close  proximity  to  the  muscle  fibre  it  controls. 
The  transmitter  substance  depolarizes  an  area  of  the  muscle  fibre 
membrane  to  a  threshold  level  initiating  an  action  potential  in  the 
muscle  fibre.  This  action  potential  causes  the  fibre  to  contract. 

In  smooth  muscle,  nerves  are  not  closely  associated  with  specific 
muscle  fibres.  The  effector  unit  is  not  the  single  muscle  cell.  It 
has  been  shown  by  Bennett  (1972)  that  it  is  difficult  or  impossible 
to  initiate  an  action  potential  in  a  single  muscle  cell.  The  response 
of  most  smooth  muscles  to  depolarization  of  their  surface  membranes 
is  a  passive  electrotonic  depolarization  (Kuriyama  &  Tomita,  1965 
and  Bennett  &  Merrillees,  1966).  In  smooth  muscle  the  effector  unit 
is  the  smooth  muscle  bundle.  It  has  been  shown  by  Bennett  (1967) 
and  Tomita  (1967)  that  a  propagating  action  potential  can  be  initiated 
if  a  large  area  of  the  muscle  is  simultaneously  depolarized,  either 
by  stimulating  the  nerves  to  the  muscle  or  by  passing  current  across 
the  muscle  with  large  extracel 1 ul ar  electrodes.  The  area  of  muscle 
which  must  be  simultaneously  depolarized  in  order  to  initiate  a  propa¬ 
gating  action  potential  indicates  that  the  effector  unit  is  the  muscle 
bundle.  The  nervous  system  causes  the  release  of  transmitter  substance 
in  the  general  area  of  a  muscle  bundle.  When  the  cells  in  the  bundle 
have  been  depolarized  to  a  threshold  level,  action  potentials  are 
initiated  and  the  bundle  contracts. 
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CHAPTER  3 


SLOW  WAVE  ACTIVITY 


3.1  INTRODUCTION 

The  majority  of  studies  of  slow  waves  have  been  done  on  the 
stomach  and  small  intestine.  This  chapter  will  first  review  what  is 
currently  known  about  the  slow  waves  occurring  in  these  organs.  It 
will  also  discuss  the  use  of  coupled  oscillator  models  to  simulate 
the  spatial  organization  of  the  slow  wave  activity.  Because  of  the 
importance  of  pressure  electrode  measurements  to  the  work  described 
in  this  thesis,  a  discussion  is  included  of  the  currently  accepted 
theory  of  how  these  electrodes  work.  Finally,  the  current  knowledge 
of  slow  wave  activity  in  the  human,  cat,  and  dog  colons  is  presented. 

3.2  SLOW  WAVES  IN  THE  STOMACH  AND  SMALL  INTESTINE 

3.2.1  Electrical  Characteri sti cs 

The  slow  waves  which  occur  in  the  smooth  muscle  wall  of  the  stomach 
have  a  single  frequency  of  one  to  three  cycles  per  minute  in  man 
in  vivo.  Intracellular  recordings  from  stomach  smooth  muscle  show 
an  electrical  complex  consisting  of  a  rapid  depolarization  followed 
by  a  plateau  upon  which  spikes  are  sometimes  imposed.  In  vitro  measure¬ 
ments  by  El -Sharkawy  et.  al .  (1978)  showed  a  cellular  resting  potential 
in  canine  and  human  stomach  averaging  -69  to  -74  mV  as  the  recording 
location  moved  from  the  orad  side  of  the  corpus  to  the  terminal  antrum. 
The  slow  wave  amplitude  increased  in  the  aboral  direction  from  21 
to  31  mV. 
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Slow  waves  in  the  small  intestine  of  the  human,  cat,  and  dog  have 
a  frequency  in  the  range  of  7  -  20  cycles/minute  (Connor,  1979).  Intracel¬ 
lular  micro-electrode  studies  from  the  small  intestine  of  the  rabbit  showed 
a  mean  resting  potential  -54.8  ±  0.3  mV  and  a  slow  wave  amplitude  of  17.9  ± 
0.2  mV  (El-Sharkawy  and  Daniel,  1975  a).  The  corresponding  measure¬ 
ments  for  the  cat  are  -67.0  ±  7.0  mV  and  30.0  ±  5.0  mV  (Connor  et.  al . , 
1977).  Characteristic  slow  wave  activity  from  the  cat  small  intestine 
is  illustrated  in  Figure  3.1. 

3.2.2  Relationship  to  Mechanical  Activity 

It  has  been  shown  (Daniel  et.  al . ,  1960;  Daniel  and  Chapman, 

1963)  that  spike  potentials  in  the  stomach  and  small  intestine  most 
frequently  occur  during  the  plateau  phase  of  the  slow  wave  as  shown 
in  Figure  3.1).  Although  contractions  can  occur  without  spike  poten¬ 
tials  (Daniel  and  Irwin,  1968;  Szurszewski,  1981),  these  potentials 
always  correspond  with  contractions.  Therefore  contractions  occur 
most  frequently  during  the  plateau  portion  of  the  slow  wave. 

Many  studies  employing  multisite  in  vivo  recordings  (Bass,  1968; 

Kelly  et.  al.,  1969;  Kelly  and  Code,  1971)  have  demonstrated  that 
slow  waves  are  responsible  for  the  highly  coordinated  contractile 
activity  of  the  stomach.  Slow  waves  are  phase  locked  over  the  stomach 
from  the  corpus  to  the  antrum,  and  the  frequency  and  relative  phase 
lag  of  these  waves  determine  the  frequency  and  velocity  of  peristaltic 
contractions  in  this  organ. 

The  relationship  of  the  small  intestine  electrical  activity  to 
the  role  of  the  intestine  musculature,  mixing  and  propelling  material, 
is  not  as  straightforward  as  in  the  stomach.  The  basic  element  of 
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Figure  3.1  Mi croel ectrode  recordings  from  the  musculature  of  cat  intestine. 
(A)  Two  examples  of  records  taken  from  longitudinal  layer 
cells.  (B)  Activity  in  circular  muscle  recorded  from  a 
1  mm  area  from  which  the  longitudinal  layer  had  been  dissected. 
Calibrations  are  the  same  in  all  records.  (From  Connor,  1979.) 
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contractile  activity  in  the  small  intestine  is  the  contraction  of  a 
ring  of  circular  muscle  fibers  a  few  centimeters  wide.  During  a  peris¬ 
taltic  wave,  a  ring  of  strong  contraction  travels  down  the  intestine  at 
the  velocity  of  the  slow  wave.  More  commonly,  rings  of  contraction 
occur  nearly  simultaneously  at  several  stationary  or  slowly  moving 
locations,  dividing  the  active  length  of  the  intestine  into  segments 
(Cannon,  1902).  The  contractions  recur  at  the  slow  wave  frequency  and 
it  has  been  suggested  that  the  wavelength+  of  the  periodic  slow  waves 
determines  the  segment  lengths  (Code  et.  al.,  1968). 


3.2.3  Site  of  Origin  of  Slow  Wave  Activity 


There  are  inconsistencies  in  the  findings  of  various  investigators 
from  studies  designed  to  determine  the  cel  1 ul ar  origin  of  slow  waves 
in  the  stomach  and  small  intestine.  These  differences  probably  arise 
because  of  the  different  ways  used  to  prepare  measurement  specimens. 

In  the  cat  small  intestine  the  longitudinal  and  circular  muscle 
layers  separate  easily  as  confirmed  by  light  and  electron  microscopy 
(Connor  et.  al.,  1974;  Prosser,  1974;  Taylor  et.  al.,  1977).  In  the 
rabbit,  the  dissection  techniques  employed  separate  the  musculature 
within  the  circular  muscle  layer  (Daniel  &  Sarna,  1978).  Studies  of 
isolated  circular  and  longitudinal  layers  from  the  cat  intestine  have 
shown  that  the  circular  layer  is  quiescent  while  the  longitudinal  layer 
generates  slow  waves  (Connor  et.  al.,  1977;  Kobayashi  et.  al.,  1966). 
These  findings  have  been  contrasted  with  data  taken  from  the  rabbit 
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intestine  by  (El-Sharkawy  and  Daniel,  1974)  and  (Daniel  and  Sarna, 

1978).  In  these  studies  slow  waves  were  found  in  both  isolated  circular 
and  longitudinal  muscle. 

In  the  stomach  there  is  agreement  that  the  slow  wave  activity 
originates  in  the  longitudinal  muscle  layer  (Szurszewski ,  1981). 

3.2.4  The  Ionic  Basis  for  Slow  Waves 

There  is  considerable  evidence  supporting  the  hypothesis  that 
an  electrogenic  Na-K  pumping  mechanism  is  responsible  for  part  of 
the  resting  potential  of  smooth  muscle  cells  in  the  gut  (Szurszewski, 
1981).  For  example,  the  measured  membrane  potential  in  guinea  pig 
taenia  coli  ranges  from  -51  to  -55  mV  while  the  membrane  potential 
calculated  from  the  Goldman  equation  (Goldman,  1943)  is  -35  mV 
(Casteels,  1970).  This  difference  in  measured  and  calculated  values 
is  caused  by  the  flow  of  pump  current  across  the  membrane  resistance. 

In  addition  to  the  Na-K  pump  there  is  evidence  of  other  metabolic 
pumping  mechanisms  responsible  for  the  active  transfer  of  Cl  and 
Ca  across  the  cell  membrane  (Brading,  1979;  Casteels  &  Van  Breeman, 
1975).  These  two  mechanisms  may  also  contribute  to  the  difference 
between  calculated  and  measured  membrane  potentials. 

It  is  generally  agreed  that  the  Na-K  pump  is  responsible  for 
slow  wave  generation  in  the  small  intestine  (Szurszewski,  1981),  but 
the  particular  ionic  mechanism  responsible  is  in  dispute.  There  are 
two  hypotheses:  The  first  states  that  the  slow  wave  results  from 
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the  oscillatory  activity  of  the  Na-K  pump.  The  pump  turns  off 
to  generate  the  positive  phase  of  the  slow  wave  and  on  again 
to  return  the  membrane  to  its  resting  potential.  Evidence 
in  support  of  this  theory  has  been  obtained  by  Connor  and  Prosser 
(1974)  and  Connor  et.  al .  (1974).  The  second  hypothesis  states  that 
slow  waves  result  from  cyclic  changes  in  membrane  permeability.  Studies 
supporting  this  hypothesis  have  been  reported  by  Mills  and  Taylor 
(1971)  and  El-Sharkawy  and  Daniel  (1975  b).  Additional  work  is  required 
to  determine  which  of  these  hypotheses  is  correct. 

3.3  MODELLING  WITH  COUPLED  RELAXATION  OSCILLATORS 

3.3.1  Introduction 

The  Hodgki n-Huxl ey  model  (Hodgkin  &  Huxley,  1952)  was  very  successful 
in  explaining  the  nerve  action  potential  in  terms  of  regenerative 
feedback  through  voltage  controlled  conductances.  This  supported  the 
use  of  relaxation  type  oscillators  in  explaining  cardiac  electrical 
activity  (Roberge  &  Nadeau,  1966).  It  was  therefore  natural  that  this 
type  of  oscillator  would  be  chosen  to  model  gastrointestinal  electrical 
activity  (Nelsen  &  Becker,  1968). 

Relaxation  type  oscillators  were  first  used  to  model  cardiac 
electrical  activity  (Van  der  Pol  &  Van  der  Mark,  1928).  The  equation 
describing  the  Van  der  Pol  oscillator  is  as  follows: 

x  -  a(l  -  x^)x  +  D x  =  0  3.1 

where  x  is  the  independent  variable  representing  the  amplitude  of 

the  oscillation,  a  is  the  damping  coefficient,  w  is  the  natural  frequency 


, 


in  radians  per  second,  and  x  and  x  are  the  first  and  second  derivatives 
of  x  respecti vely .  Fitzhugh  (1961)  used  Lienard's  transformation  to 
obtain  a  system  of  two  first  order  differential  equations  as  follows: 
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,  I  •  x  2  • 

then  y=— +xx-  x 

J  Oi 

using  equation  3.1 
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3.2 

3.3 


3.4 


Terms  were  then  added  to  equation  3.4  to  give 

12 

y  =  -— (cox-a  +  by)  3.5 

Also  from  equation  3.2 

3 

x  =  a(y  +  x  -  —)  3.6 


Equations  3.5  and  3.6  are  the  Fitzhugh  system  of  equations.  Sarna  et. 
al .  (1971)  further  generalized  these  equations  introducing  other 
coefficients  and  higher  powers  of  x  as  follows: 


2  3 

x  =  a(cy  +  fx  +  gx  +  hx  ) 

1  2  2  3 

y  =  -  —  (by  +  oj  x  +  cx  +  dx  -  a) 


3.7 

3.8 


The  smooth  muscle  of  the  gastrointestinal  tract  does  not  function 
as  a  syncytium  such  as  cardiac  muscle  where  the  electrical  activity 
of  the  complete  muscle  is  controlled  from  a  single  pacemaker  site. 
Instead,  any  small  region  of  the  smooth  muscle  wall  of  the  stomach, 
small  intestine  and  colon  generates  slow  wave  activity.  The  muscle 
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can  be  considered  to  consist  of  islands  of  electrically  oscillating 
tissue  electrically  coupled  to  other  surrounding  regions  which  are  also 
electrically  active.  As  a  consequence  of  this  the  electrical  activity 
of  the  stomach  and  small  intestine  has  been  modelled  using  an  array 
of  coupled  relaxation  oscillators  as  shown  in  Figure  3.2  where  the 
nth  oscillator  is  described  by  the  following  equations  (Sarna  et.  al . , 

1  971  ,  1  972  ): 


xn  =  k(3lyn  +  a2xn  +  a3xn  +  +  I  ) 


=  -  ¥(Vn  +  b2Xn  +  V,/  +  ’  b0> 


3.9 

3.10 


where  I  represents  the  resultant  of  all  other  oscillators  feeding 
into  the  nth  oscillator. 

When  any  two  of  these  oscillators  with  different  intrinsic  frequen¬ 
cies  are  coupled  together  a  number  of  different  results  are  possible 
depending  on  the  oscillator  and  coupling  parameters.  The  two  coupled 
oscillators  could  oscillate  at  a  frequency  (1)  higher  than  either  of 
the  intrinsic  frequencies,  or  (2)  equal  to  the  higher  of  the  two 
frequencies,  or  (3)  somewhere  between  the  two  intrinsic  frequencies. 

In  each  of  the  above  situations  the  two  oscillators  become  phase  locked 
with  a  fixed  phase  lag  between  the  two  oscillator  outputs.  If  the 
degree  of  coupling  is  small,  the  two  oscillators  could  oscillate  at 
different  frequencies  making  a  fixed  phase  impossible.  In  this  case 
the  oscillator  with  the  lower  intrinsic  frequency  may  have  its  frequency 
increased  due  to  the  coupling.  This  phenomena  is  referred  to  as 
"frequency  pulling". 
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Figure  3.2  Block  diagram  illustrating  the  arrangement  of  oscillators 
in  the  intestinal  model.  (From  Sarna,  1971.) 
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3.3.2  Modelling  of  Slow  Wave  Activity  in  the  Canine  Stomach 

Slow  wave  activity  in  the  canine  stomach  has  been  modelled  using 
a  two  dimensional  array  of  16  bidi rectional ly  coupled  oscillators  as 
shown  in  Figure  3.3  (Sarna  et.  al.,  1972).  The  slow  waves  in  the 
intact  organ  in  vivo  have  a  single  frequency  of  5.1  ±  0.52  cycles/minute. 

As  seen  in  Figure  3.4  the  intrinsic  oscillator  frequencies  decrease  from 
a  maximum  of  5.5  cycles/min.  in  the  corpus  to  2.5  cycles/min.  near  the 
pylorus.  In  the  intact  organ  and  in  the  model  the  slow  wave  signals 
were  shown  to  be  phase  locked  with  the  phase  angle  varying  from  70  -  100°/cm 
in  the  corpus  to  8  -  20°/crn  near  the  pylorus.  The  degree  of  coupling 
between  oscillating  units  increases  distally  to  a  maximum  near  the 
pyl  orus . 

3.3.3  Modelling  of  Slow  Wave  Activity  in  the  Canine  Small  Intestine 

Several  researchers  have  produced  models  of  the  slow  wave  activity 
of  the  canine  small  intestine  (Nelsen  &  Becker,  1968;  Diamant  et.  al., 

1970;  Sarna  et.  al.,  1971).  The  proximal  60  to  80  cm  of  the  small 
intestine  produces  slow  waves  of  constant  frequency  in  the  range  of 
18.5  to  21.0  cycles/min.  (Sarna  et.  al . ,  1971).  In  this  plateau  region 
of  constant  frequency,  slow  waves  are  phase  locked  with  a  phase  lag  of 
10  -  1 5°/cm.  Distal  to  the  plateau  region  the  frequency  gradually 
decreases  to  a  minimum  of  14  cyles/min.  in  the  ileum.  Coupling  in  the 
circular  direction  is  good  with  slow  waves  appearing  to  be  simultaneous 
at  points  around  the  circumference  of  the  organ.  Intrinsic  frequencies 
in  small  intestine  segments  decrease  exponentially  as  a  function  of 
distance  from  the  pylorus,  and  the  frequency  of  the  intact  plateau  region 
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Figure  3.3  Arrangement  of  oscillators  in  the  gastric  ECA  model. 

Oscillators  1  to  6  represent  the  ECA  of  the  greater  curva¬ 
ture;  oscillators  7  to  11  represent  the  ECA  of  the  midline; 
and  oscillators  10-13  represent  the  ECA  near  the  lesser 
curvature.  Numbers  inside  the  boxes  indicate  the  intrinsic 
frequencies  of  oscillators,  and  those  between  them  indicate 
coupling  factors.  (From  Sarna,  1971.) 
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Figure  3.4  Diagram  showing  the  intrinsic  frequency  gradients  along 

the  greater  curvature  ( — )  and  on  the  midline  (--)  of  the 
dog  stomach.  Numbers  refer  to  different  electrodes. 

( From  Sarna  ,  1 971 . ) 
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is  higher  than  the  highest  intrinsic  frequency  from  isolated  segments 
(Sarna  et.  al  . ,  1971).  The  model  developed  by  Sarna  et.  al .  (1971) 
consisted  of  a  series  string  of  bidirectional ly  coupled  oscillators 
with  intrinsic  frequencies  decreasing  distally  in  the  same  manner  as 
indicated  by  in  vivo  measurements.  Coupling  was  also  made  to  decrease 
as  a  function  of  the  oscillator  location  in  the  string.  The  model  was 
shown  to  exhibit  the  frequency  and  phase  lag  characteristics  of  the 
small  intestine.  In  addition,  the  model  was  shown  to  predict  the 
changes  in  slow  wave  activity  produced  by  transection  of  the  organ  at 
various  locations  along  its  length. 

3.4  INTRACELLULAR-TO-EXTRACELLULAR  SIGNAL  TRANSDUCTION 

3.4.1  Introduction 

Intracellular  mi croel ectrode  measurements  from  smooth  muscle  cells 
are  very  difficult  to  obtain  and  multisite  i ntracel  1  ul ar  recordings  have 
not  been  achieved  in  smooth  muscle.  For  these  reasons,  the  majority  of 
in  vitro  slow  wave  studies  has  used  pressure  electrodes  of  the  type  des¬ 
cribed  by  Bortoff  (1967).  A  large  portion  of  the  data  in  this  thesis 
was  obtained  from  pressure  electrode  measurements.  A  knowledge  of  the 
i ntracel  1  ul ar-to-extracel 1 ul ar  transduction  mechanism  is  necessary  to 
relate  these  extracel 1 ul ar  measurements  to  intracellular  phenomena. 

It  is  generally  believed  that  this  mechanism  results  from  the  generation 
of  an  injury  potential. 
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3.4.2  The  Injury  Potential  Theory 

In  1946  Graham  and  Gerard  demonstrated  that  the  impalement  of  muscle 
cells  by  large  microelectrodes  resulted  in  the  measurement  of  an  intra¬ 
cellular  potential  which  was  far  less  than  the  expected  membrane  potential. 
It  was  thought  that  an  injury  to  the  cell  caused  a  drop  in  the  intra¬ 
cellular  potential  and  thus  the  measured  potential  was  referred  to  as 
an  "injury  potential".  This  was  again  demonstrated  for  cardiac  muscle 
by  Woodbury  et.  al.  (1951)  who  showed  that  the  measured  potential 
would  fall  to  as  low  as  30%  of  the  correct  value.  Later  Gillespie 
(1962)  reported  a  detailed  study  of  this  phenomena  for  the  smooth 
muscle  cells  of  the  rabbit  large  intestine.  He  showed  that  in  addition 
to  measuring  a  reduced  resting  potential  at  the  site  of  the  injury, 
the  action  potential  at  this  site  was  reduced  by  the  same  factor.  He 
also  showed  that  mechanical  deformation  of  the  cell  due  to  stretching 
reduced  the  intracel  1  ul ar  potential.  From  these  observations  he  concluded 
that  mechanical  deformation  of  a  cell  by  an  electrode  too  large  to 
penetrate  it  would  cause  depolization  and  produce  an  injury  potential. 

He  developed  the  simple  model  shown  in  Figure  3.5  to  explain  his  results. 

In  this  model  R  is  the  normal  cell  membrane  resistence  and  vm  is  the 

m  Hi 

normal  i ntracel 1 ul ar  potential.  is  the  shunt  resistance  existing 
between  the  interior  and  exterior  of  the  electrode,  and  R^  is  a  variable 
membrane  leakage  resistance  created  by  the  injury.  In  the  extreme 
case  Rl  will  be  zero.  For  this  situation  the  measured  potential  v  is 
given  by  the  following  equation: 


v 


3.11 
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Figure  3.5 


The  injury  potential  generated  by  an  extracel 1 ul ar  electrode, 
and  an  equivalent  circuit  of  the  extracellular  signal 
generation.  (From  Gillespie,  1962.) 
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Of  course,  if  the  i ntracel 1 ul ar  potential  has  a  time  varying  component, 
such  as  a  slow  wave,  then  v  will  contain  a  proportional  time  varying 
component. 

Bortoff  (1967)  proposed  that  this  same  theory  was  applicable  to 
the  smooth  muscle  of  the  small  intestine  when  deformed  by  large  extra¬ 
cellular  pressure  electrodes. 

3.4.3  An  Extension  of  Gillespie's  Model 

Extracel 1 ul ar  pressure  electrodes  are  much  bigger  than  the  cells 
which  are  being  studied.  Figure  3.6(a)  schematically  shows  the  typical 
measurement  situation.  Figure  3.6(b)  is  the  equivalent  circuit.  The 
components  of  this  circuit  no  longer  represent  the  parameters  of  a 
single  cell  as  in  Gillespie's  model.  They  represent  the  equivalent 
parameters  of  many  cells  having  the  same  characteri sti cs .  It  is  assumed 
that  the  membrane  potential  vm  is  composed  of  a  resting  potential  Vm 
and  a  slow  wave  component  vsw>  From  Gillespie's  work  it  is  assumed 
that  the  cells  deformed  by  the  electrode  are  depolarized  by  the  creation 
of  a  membrane  leakage  resistance  for  each  of  these  cells.  in  Figure 
3.6(b)  represents  the  equivalent  leakage  resistance  of  all  cell  mem¬ 
branes  deformed  by  the  electrode.  It  is  assumed  that  the  internal 
resistances  of  the  cells  and  of  the  extracel 1 ul ar  fluid  of  the  tissue 
bath  are  both  negligible.  R^  is  the  equivalent  shunt  leakage  resistance 
existing  between  the  measurement  surface  of  the  electrode  and  the 
extracellular  fluid  of  the  tissue  bath.  Rc  is  the  equivalent  coupling 
resistance  existing  between  the  injured  and  normal  cells.  In  the 
extreme  case  the  cells  under  the  electrode  will  be  completely  depolarized 
and  Rl  will  approach  zero.  In  this  case  the  measured  potential  v  is 
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Figure  3.6  (a)  A  schematic  drawing  of  a  pressure  electrode  measurement, 

(b)  An  equivalent  circuit. 
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given  by: 


3.12 


Since  vm  =  vm  +  vsw>  the  measured  potential  will  have  a  time  varying 
component  proportional  to  v  ,  that  is 


sw 


-R 

s 


3.13 


Vac  R  +  R  +  R  vsw 
m  c  s 


where  v__  is  the  time  varying  component  of  v. 

No  direct  comparisons  of  simultaneous  intracellular  and  extra¬ 
cellular  pressure  electrode  measurements  have  been  reported  to  support 
this  theory.  However,  part  of  the  work  described  in  this  thesis  shows 
that  the  theory  described  above  does  not  directly  apply  to  the  smooth 
muscle  of  the  colon.  An  alternative  theory  will  be  proposed  in  this 
thesis  to  explain  the  generation  of  pressure  electrode  signals. 

3.5  SLOW  WAVE  ACTIVITY  IN  THE  COLON 

3.5.1  Introduction 

Less  is  understood  about  the  function  of  slow  waves  in  the  colon 
than  in  either  the  stomach  or  small  intestine.  Slow  waves  have  been 
recorded  from  the  colonic  smooth  muscle  of  the  dog  (El -Sharkawy , 

1978),  the  cat  (Christensen,  1975),  the  dog  (Vanasin  et.  al . ,  1974; 
Bowes  et.  al . ,  1  977)  and  the  human  (Vanasin  et.  al  . ,  1  974; 

Duthie,  1975;  Sarna  et.  al.,  1980;  Chambers  et.  al.,  1981). 

However  many  inconsistencies  occur  in  these  reports  either  due  to 
species  differences  or  differences  in  measurement  procedures. 
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3.5.2  Slow  Waves  in  the  Human  Colon 

In  vitro  studies  of  human  slow  wave  activity  are  hindered  by  the 
difficulty  of  obtaining  healthy  tissue  specimens  and  in  vivo  studies 
are  infrequently  attempted  because  of  the  post  operative  risks  to 
patients.  There  is  lack  of  agreement  on  human  colonic  slow  wave  activity 
with  respect  to  its  omnipresence,  its  frequency,  and  frequency  gradients. 
Taylor  et.  al  .  (1  975)  and  Snape  et.  al .  (1976)  reported  that  the  slow 
wave  is  i ntermi ttently  present  in  the  in  vivo  human  colon.  However, 
Provenzale  and  Pisano  (1971)  found  that  it  was  omnipresent.  A  number 
of  researchers  reported  two  distinct  slow  wave  frequencies  one  of  3 
and  one  of  12  cycles/minute  (Vanasin  et.  al . ,  1974;  Duthie,  1975; 

Stoddard  et.  al . ,  1979;  Sarna  et.  al.,  (1980).  However, 

Chambers  et.  al .  (1981)  report  a  single  frequency  of  2-9  cyles/ 
minute.  Vanasin  et.  al .  (1  974)  reported  the  existence  of  two 
frequency  gradients  with  the  frequency  decreasing  distally 
in  the  ascending  colon  while  increasing  from  transverse  colon 
to  rectum.  However,  Chambers  et.  al .  (1981)  report  that  no  frequency 
gradient  could  be  detected. 

No  data  has  been  published  concerning  the  i ntracel 1 ul ar  slow  wave 
characteristics  or  concerning  the  typical  phase  angle  existing  between 
oscillating  regions  in  the  human  colon. 

3.5.3  Slow  Waves  in  the  Cat  Colon 

The  electrical  activity  of  the  cat  colon  has  been  extensively 
studied  by  Christensen  et.  al .  (1969).  In  this  animal  the  slow  wave 
has  been  shown  to  originate  in  the  circular  muscle  layer  and  to  be 
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nonexistant  in  the  longitudinal  layer  (Christensen  et.  al.,  1969).  This 
electrical  activity  is  phase  locked  67  percent  of  the  time  in  the  long¬ 
itudinal  direction  and  95  percent  of  the  time  in  the  circular  direction. 

The  mean  frequency  in  isolated  muscle  segments  was  determined  to  be 
3.12  ±  1.62  cycles/minute  (Christensen  &  Hauser,  1971  a;  1971  b).  The 
phase  lag  between  recording  sites  varies  linearly  with  electrode  spacing 
and  the  size  of  the  oscillating  regions  is  believed  to  be  less  than 
0.5  mm  (Christensen  &  Rasmus,  1972).  The  complete  organ  was  studied 

in  vitro  by  Christensen  et.  al .  (1974).  From  this  study  it  was  found 
that  a  frequency  gradient  exists  in  the  aboral  direction  with  an  intrinsic 
frequency  of  7.3  cycles/minute  at  the  distal  end  and  1.9  cycles/minute 
at  the  cecum.  In  addition  to  the  slow  wave,  migrating  spike  bursts 
were  also  recorded.  The  spike  bursts  occurred  every  74.3  seconds  with 
a  mean  duration  of  26.6  seconds  and  propagated  distally  over  the  complete 
organ . 

3.5.4  Slow  Waves  in  the  Canine  Colon 

Less  has  been  reported  on  the  canine  colon  than  on  either  the 
human  or  cat  colon,  and  some  of  the  reports  have  been  contradictory. 

Vanasin  et.  al  .  (1  974)  reported  the  existence  of  a  range  of  slow  wave 
frequencies  with  a  frequency  gradient  existing  from  the  cecum  to,  the 
rectum.  However  Shearin  et.  al.  (1979)  found  no  significant  gradient 
of  intrinsic  frequencies.  Coupling  has  been  found  to  be  better  in  the 
circular  direction  than  in  the  longitudinal  direction  (El -Sharkawy, 

1978;  Shearin  et.  al.,  1979).  All  researchers  have  found  the  slow  waves  to 
be  omni-present  with  frequencies  in  the  range  from  4.1  to  6.2  cycl es/mi nuxe. 
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Shearin  et.  al  .  (1979)  determined  mean  phase  angles  between  recording 
sites  using  visual  and  computer  analysis.  However,  these  results  are 
to  be  questioned  because  there  was  not  sufficient  agreement  between 
the  visual  and  computer  results. 

There  are  a  number  of  questions  concerning  slow  waves  in  the  canine 
colon  which  have  not  been  answered  in  the  literature.  The  site  of  origin 
of  the  slow  wave  is  not  known.  What  is  the  role  played  by  each  of  the 
muscle  layers  in  the  generation  of  slow  waves?  Are  slow  waves  present 
in  the  cells  of  both  layers?  In  addition  to  these  points,  the  intra¬ 
cellular  characteristi cs  of  the  slow  wave  in  the  canine  colon  have  not 
been  reported.  The  work  described  in  this  thesis  addresses  each  of 
the  above  points. 

It  would  have  been  most  desirable  to  study  human  colonic  muscle. 
However,  the  uncertainty  concerning  the  availability  of  healthy  specimens 
prevented  this.  As  stated  previously  the  dog  has  proven  to  be  a  most 
suitable  animal  for  the  study  of  stomach  and  small  intestinal  slow  wave 
activity  (Sarna  et.  al  . ,  1971;  1  972  )  and  a  base  of  knowledge  had  been 
obtained  from  previous  studies  in  our  laboratory  on  the  canine  colon 
(Shearin  et.  al.,  1979;  Bowes  et.  al.,  1977).  Thus,  the  dog  was  the 
obvious  animal  on  which  to  base  the  present  research. 
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CHAPTER  4 


OBJECTIVES,  METHODS,  AND  MATERIALS 
4.1  INTRODUCTION 

From  initial  preliminary  studies  on  the  canine  colon  it  was  deter¬ 
mined  that  the  distal  colon  was  completely  unpredictable  in  its  ability 
to  produce  slow  wave  activity.  Consequently,  this  work  is  based  on 
specimens  taken  from  the  proximal  one  half  to  two  thirds  of  the  canine 
colon.  The  distal  portion  of  the  colon  will  require  future  study. 

Healthy  mongrel  dogs  weighing  between  8  and  15  kilograms  were 
anesthetized  using  pentobarbital,  and  the  required  sections  of  the 
colon  were  surgically  removed.  These  sections  were  immediately  placed 
in  oxygenated  Krebs-Ringer  solution^  and  cut  open  either  along  the 
mesenteric  or  antimesenteri c  border.  These  sections  were  carefully 
washed  in  three  separate  containers  of  Krebs-Ringer  solution  to  insure 
the  complete  removal  of  fecal  contents.  The  sections  were  then  placed 
in  a  constant  temperature  bath  of  Krebs-Ringer  solution  and  the  mucosa 
and  muscularis  mucosa  were  removed  using  sharp  dissection.  The  sections 
were  then  trimmed  and  cut  into  circularly  or  longitudinally  oriented 
strips  (10  mm  x  40  mm)  for  specific  experiments. 

The  objectives  of  the  research  are  as  follows: 

1.  The  determination  of  the  anatomical  origin  of  colonic  slow  wave 
acti vi ty . 

2.  The  definition  of  the  role  played  by  each  of  the  muscle  layers  in 
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the  generation  of  slow  wave  activity. 

3.  A  more  accurate  determination  of  the  physical  and  electrical  charac¬ 
teristics  of  slow  wave  activity.  This  will  include  an  examination 
of  i ntracel 1 ular  potentials  to  determine  the  resting  potentials 
and  slow  wave  amplitudes.  It  will  also  include  an  examination  of 
slow  wave  frequencies,  of  frequency  gradients,  of  the  size  of  oscillating 
regions,  and  of  the  relative  phase  angles  existing  between  regions 
in  the  circular  and  longitudinal  directions. 

An  additional  objective  resulting  from  the  extensive  use  of  pressure 
electrodes  is  to  determine  the  relationship  between  pressure  electrode 
signals  and  the  i ntracel 1 ul ar  potential  variations  of  cells  under  these 
el ectrodes . 

Extracel 1 ul ar  pressure  electrodes  are  essential  to  much  of  the 
work  in  this  thesis.  This  chapter  first  includes  a  description  of  the 
development  and  study  of  the  pressure  electrode  used  in  this  work. 

This  is  followed  by  a  detailed  description  of  extracellular  and  intra¬ 
cellular  experiments  and  of  the  computer  processing  used  to  analyze 
the  results. 

4.2  THE  DESIGN  AND  STUDY  OF  A  NEW  PRESSURE  ELECTRODE 

The  majority  of  slow  wave  studies  have  been  completed  using  extra¬ 
cellular  pressure  electrodes  (Bortoff,  1967;  1975;  Christensen  et.  al . , 

1969;  Shearin  et.  al.,  1979).  This  is  due  to  the  fact  that  intracellular 
measurements  are  very  difficult  to  obtain  in  comparison  to  extracel 1 ul ar 
measurements.  Intracel 1 ul ar  measurements  are  much  more  time  consuming, 
require  more  sophisticated  measurement  apparatus,  and  require  more 
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tedious  measurement  procedure.  In  addition  to  this,  multisite  intra¬ 
cellular  measurements  are  very  difficult  to  make  whereas  the  number  of 
simultaneous  extracel 1 ul ar  recordings  is  limited  only  by  the  size  of 
the  tissue  and  the  availability  of  recording  channels. 

Nevertheless,  extracellular  measurements  have  a  number  of  limita¬ 
tions.  The  most  important  of  these  is  that  no  information  can  be 
obtained  about  the  resting  (dc)  potentials  existing  in  cells  under  the 
electrode.  A  second  is  that  at  present  it  is  impossible  to  calibrate 
the  electrode  signals  in  terms  of  actual  potential  variations  inside 
the  cells.  The  pressure  electrode  theory  outlined  in  the  previous 

chapter  has  not  yet  been  proven,  and  even  if  it  were  correct,  it  does 

0 

not  provide  the  necessary  parameters  to  relate  i ntracel 1 ul ar  and  extra¬ 
cellular  potentials.  In  addition  to  the  above  problems,  extracellular 
electrodes  are  limited  by  noise  and  low  frequency  drift.  Most  extra¬ 
cellular  pressure  electrodes  utilize  a  si  1 ver-si 1 verchl ori de  layer  as 
the  signal  transducer.  The  resistance  and  capacitance  of  this  type  of 
electrode  increase  as  an  inverse  power  of  frequency.  This  characteri sti c 
could  cause  severe  distortion  of  signals  such  as  slow  waves  which 
have  very  low  frequency  components.  Typical  slow  wave  frequencies  in 
the  dog  colon  are  in  the  range  of  4  to  7  cycles/min.  (0.067  to  0.117 
Hz)  in  vitro. 

4.2.1  Electrode  Description 

Our  work  on  the  colon  started  with  a  glass  pore  electrode  of  the 
type  described  by  Christensen  &  Hauser  (1971  a).  However,  this  electrode 
was  inadequate  because  of  its  high  impedance,  its  large  offset  voltage 
and  its  1 ow-frequency  drift.  Therefore,  the  development  of  a  better 


. 


42 


electrode  was  initiated.  Si  1 ver-si 1 verchl ori de  electrodes  have  been 
used  in  many  physiological  applications  (Geddes,  1972)  and  have  been 
shown  to  have  good  impedance,  noise,  and  drift  characteri sties  (Geddes 
&  Baker,  1968).  These  facts  led  to  the  design  of  an  electrode 
using  a  chlorided  silver  wire  as  its  transducing  element. 

The  electrode  is  shown  in  Figure  4.lt  It  consists  of  a  1  mm 
diameter  silver  wire  inside  a  rigid  glass  capillary  with  an  outside 
diameter  of  2  mm.  During  electrode  application  the  capillary  provides 
a  seal  between  the  conducting  extracellular  fluid  and  the  active  elec¬ 
trode  surface  ensuring  that  the  shunt  resistance  R$  is  very  large. 

As  determined  by  equation  3.13  in  Chapter  3,  the  larger  the  magnitude 
of  Rs  the  larger  the  recorded  signal  and  as  a  consequence  the  better 
the  si gnal - to-noi se  ratio.  The  surface  area  of  the  electrode  was 
maximized  by  roughening,  and  it  was  optimally  chlorided  to  minimize 
electrode  impedance  and  noise  characteri sti cs .  The  wire  was  fixed 
inside  the  capillary  with  epoxy  resin. 

4.2.2  Electrode  Impedance  Minimization 

As  stated  above,  the  effective  surface  area  of  the  electrode  was 
increased  before  it  is  chlorided.  This  was  accomplished  by  roughening 
the  end  surface  of  the  silver  wire.  It  was  emperically  determined  that 
a  maximum  increase  in  effective  surface  area  was  obtained  when  the 
surface  was  roughened  with  an  aluminium  oxide  abrasive  with  a  mesh  size 
of  240.  It  was  difficult  to  measure  the  effective  surface  area,  but  by 
comparison  of  the  capacitance  of  a  roughened  electrode  to  that  of  a 
polished  electrode  it  was  estimated  that  the  area  was  increased  by  a 
factor  of  five  or  better. 


+Thi s  electrode  design  was  originally  proposed  by  Prof.  Y.J.  Kingma. 
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1mm  dia.  silver  wire 
2mm  dia.  glass  capillary 


ep^xy  resin 


AgCI  layer 


Figure  4.1  Pressure  Electrode  Details 


Reports  from  the  literature  vary  greatly  with  regard  to  the  optimum 
chloriding  levels  required  to  produce  minimum  impedance  sil ver-si 1 ver- 
chloride  electrodes  (Geddes,  1972).  The  applicability  of  these  reports 
to  our  electrodes  is  questionable  because  these  reports  have  been  for 
much  larger  electrodes  and  for  frequencies  of  10  Hz  and  greater.  Thus 
it  was  necessary  to  determine  optimum  chloriding  levels  and  to  measure 
the  electrode  impedance  characteri sties  at  low  frequencies. 

Figure  4.2  is  a  plot  of  impedance  versus  chloriding  level  for  five 
different  frequencies  from  0.1  Hz  to  10  KHz.  The  chloriding  level  is 
indicated  as  a  product  of  chloriding  current  amplitude  multiplied  by 
the  application  time  in  seconds.  For  frequencies  at  10  Hz  and  greater 
these  results  agree  with  data  published  by  Geddes  (1972).  For  0.1  Hz 
the  optimum  chloriding  level  is  3  mA-seconds.  This  was  chosen  as  the 
optimum  level  for  the  electrodes. 

4.2.3  Electrode  Impedance  Measurement 

To  measure  the  impedance  of  the  electrode  it  was  placed  in  a  0.1 
Normal  saline  solution  at  a  fixed  distance  from  a  large  sil ver-si 1 ver- 
chloride  reference  electrode  of  known  negligible  impedance.  The  impedan 
across  the  two  electrodes  was  measured  using  two  different  methods  of 
measurement  for  which  the  frequency  ranges  overlapped.  For  frequencies 
above  5  Hz  a  vector  impedance  meter  was  used.  For  frequencies  in  the 
range  from  0.1  Hz  to  20  Hz  the  impedance  was  determined  by  passing  a 
sinusoidal  current  of  peak  density  25  yA/cm  through  the  electrode¬ 
reference  pair. 

The  noise  amplitude  of  the  electrode  was  also  measured  with  respect 
to  a  large  sil ver-sil verchloride  reference  electrode  in  a  0.1  Normal 
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Figure  4.2  Electrode  impedance  as  a  function  of  chloriding  levels 
and  frequency. 
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sal i ne  bath . 

4.3  EXTRACELLULAR  STUDIES 

The  study  of  the  oscillating  regions  of  the  colon  and  of  their 
slow  waves  requires  the  simultaneous  recording  of  electrical  signals 
from  a  number  of  sites  on  the  muscle  wall.  As  stated  previously,  to 
obtain  multisite  measurements  it  was  necessary  to  use  the  extracellular 
electrodes  previously  described. 

Two  sets  of  experiments  were  designed;  one  set  to  examine  coupling 
characteri sti cs  between  various  oscillating  regions  and  the  other  to 
examine  the  role  played  by  each  of  the  muscle  layers  in  the  generation 
of  slow  waves.  The  first  set  of  experiments  is  summarized  in  Figure  4.3. 
This  set  of  experiments  is  composed  of  two  groups.  One  group  of  experi¬ 
ments,  conducted  to  measure  average  relative  phase  angles  as  a  function 
of  electrode  spacing,  is  described  in  sections  4.3  and  4.3.2.  The 
second  group,  conducted  to  study  frequencies  and  phase  angles  before 
and  after  tissue  cuts,  is  described  in  section  4.3.3. 

The  second  set  of  experiments,  designed  to  study  the  role  played 
by  the  individual  muscle  layers,  is  also  composed  of  two  groups  of 
measurements  as  shown  in  Figure  4.4.  The  first  group,  isolated  layer 
studies,  is  described  in  section  4.3.4,  while  the  second  group,  the 
dual  chamber  studies,  are  described  in  section  4.3.6. 

4.3.1  Multisite  Electrode  Studies  on  the  Intact  Muscle  Wall 

This  set  of  experiments  consists  of  two  groups  of  measurements 
using  four  to  eight  extracellular  electrodes.  The  first  group  is  a 
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Figure  4.3  Multisite  electrode  experiments  designed  to  examine  coupled 
oscillator  characteristics . 
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Figure  4.4  Experiments  designed  to  determine  the  origin  of  the  slow 
wave  si gnal . 


.  •- 


49 


study  of  the  relative  phase  angles  in  the  longitudinal  and  circular 
directions  as  the  inter-electrode  spacing  is  varied  from  2.5  mm  to 
7.5  mm.  The  second  group  is  a  study  of  signal  frequencies  and  relative 
phase  angles  before  and  after  the  tissue  has  been  cut  between  recording 
si tes . 

The  first  objective  of  this  set  of  experiments  was  the  generation 
of  a  data  base  of  frequencies  at  and  relative  phase  angles  and  between 
recording  sites  on  the  transverse  and  proximal  colon.  This  data  would 
be  used  to  determine  the  size  of  groups  of  cells  which  oscillate  together, 
to  determine  if  a  frequency  gradient  exists,  and  to  determine  differences 
in  coupling  between  the  longitudinal  and  circular  directions.  The 
second  objective  was  to  determine  if  the  region  of  the  mesentery  has  a 
role  with  respect  to  slow  wave  generation. 

4.3.2  Phase  Angles  as  a  Function  of  Electrode  Spacing 

Six  experiments  were  conducted  to  study  the  relative  phase  angle 
as  a  function  of  the  spacing  between  recording  sites.  Strips  of  intact 
colonic  muscle  wall  were  prepared  as  described  in  Section  4.1  and 
placed  in  a  constant  temperature  bath  as  shown  in  Figure  4.5  with  the 
mucosal  side  facing  up.  The  strips  were  cut  either  in  a  circular  or 
longitudinal  orientation.  The  specimen  was  superfused  with  oxygenated 
Krebs-Ringer  solution  at  a  rate  of  6  ml/minute  and  the  temperature 
maintained  at  37.5°C. 

The  combined  electrode  holder  and  force  transducer  shown  in 
Figure  4.5  was  used  to  apply  from  two  to  four  electrodes  to  the  tissue 
with  an  average  force  of  0.059  newtons  (o  grams)  per  electrode.  The 
electrode  spacing  was  set  at  2.5  mm,  5.0  mm,  or  7.5  mm. 
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Figure  4.5  The  constant  temperature  bath  and  extracellular  measurement 
apparatus . 
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In  these  six  experiments  twelve  sets  of  measurements  were  made 

in  the  circular  direction  and  eighteen  sets  of  measurements  in  the 
longitudinal  direction. 

In  addition  to  these  six  experiments  data  was  obtained  from  the 
basal  recordings  of  seventeen  experiments  conducted  for  other  purposes. 

All  of  the  data  were  recorded  on  a  Type  R  Beckman  polygraph  as 
well  as  on  an  fm  tape  recorder  for  later  computer  processing  as  des¬ 
cribed  in  section  4.6. 

4.3.3  The  Cutting  Experiments 

In  a  system  of  poorly  coupled  oscillators  as  exists  in  the  canine 
colon  it  is  difficult  to  quantify  the  degree  of  coupling.  Also  it 
is  difficult  to  determine  if  a  gradient  exists  in  the  intrinsic  frequencies 
of  individual  oscillators.  In  a  tightly  coupled  system,  as  exists  in 
the  stomach,  separation  of  oscillators  by  cutting  has  proven  very 
effective  in  determining  the  frequency  gradient  as  well  as  in  providing 
insight  into  the  type  and  degree  of  coupling  (Sarna,  1971). 

Twenty  cutting  experiments  were  conducted  on  intact  muscle  segments. 
Ten  of  these  experiments  were  on  muscle  strips  oriented  in  the  longitu¬ 
dinal  direction  and  ten  were  on  strips  oriented  in  the  circular  direc¬ 
tion.  For  each  experiment,  circular  and  longitudinal  strips  were 
prepared  as  previously  described  in  Section  4.1  and  mounted  in  the 
constant  temperature  bath  of  Figure  4.5.  After  a  period  of  30  minutes, 
measurements  from  three  to  four  sites  were  recorded  on  fm  tape  and 
on  a  polygraph  recorder.  These  basal  recordings  werp  usually  1 5  to 
30  minutes  in  duration.  Upon  completion  of  the  basal  recording  the 
specimen  was  cut  between  two  recording  sites  without  moving  either 
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the  electrodes  or  the  specimen.  Of  the  ten  experiments  on  circularly 
oriented  strips,  five  utilized  cuts  on  the  mesenteric  border  and  the 
remainder  cuts  1  cm  from  the  mesentery.  Thirty  minutes  after  the 
cutting  procedure  a  second  set  of  recordings  was  obtained. 

The  records  from  "before"  and  "after"  cutting  were  analyzed  visually 
and  using  the  computer  programs  described  in  section  4.6.  Also  the 
computer  was  used  to  calculate  the  average  coherence  (Bendat  and 
Piersol,  1966)  for  a  band  of  descrete  frequencies  around  the  signal 
frequency.  The  coherence  function  is  defined  as  the  Fourier  Transform 
of  the  cross-correl a ti on  of  two  signals  divided  by  the  product  of  the 
auto-correl ation  Fourier  Transforms.  The  average  coherence  is  a 
measure  of  the  correlation  and  coupling  between  adjacent  channels. 

For  completely  coupled  signals  the  coherence  value  is  unity,  and  for 
completely  independent  signals  it  is  zero. 

4.3.4  Studies  of  Isolated  Muscle  Layers 

To  determine  the  role  played  by  each  of  the  smooth  muscle  layers 
in  the  generation  of  the  slow  wave  signal,  it  was  necessary  to  determine 
the  electrical  activity  of  each  of  the  isolated  muscle  layers.  To  do 
this,  a  dissection  technique  was  developed  for  separating  the  longitu¬ 
dinal  and  circular  muscle  layers.  In  developing  this  technique  it 
was  discovered  that  a  very  thin  layer  of  submucosa  could  be  easily 
removed  from  the  mucosal  side  of  the  circular  muscle  layer.  Conse¬ 
quently,  the  electrical  activities  of  the  intact  muscle  and  isolated 
circular  layer  were  recorded  with  and  without  this  submucosal  layer. 
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4.3.5  Separation  and  Study  of  Longitudinal  and  Circular  Muscle  Layers 

To  separate  the  circular  and  longitudinal  muscle  layers  a  strip 
of  intact  colon  wall  was  mounted  with  its  serosal  side  facing  down  in 
a  constant  temperature  bath  of  oxygenated  Krebs-Ringer  solution.  The 
mucosa  was  removed  using  sharp  dissection  and  the  tissue  was  trimmed 
to  form  a  1  cm  by  2  cm  rectangle  with  edges  parallel  to  the  circular 
and  longitudinal  directions.  The  tissue  was  turned  over  and  at  an  edge 
parallel  to  the  circular  muscle  direction  the  two  layers  were  teased 
apart  using  iris  scissors.  When  an  edge  of  longitudinal  muscle 
(sufficiently  large  to  be  gripped  with  forceps)  was  developed,  small 
strips  of  longitudinal  muscle  2  to  4  mm  in  width  were  gently  peeled 
back  leaving  behind  the  circular  muscle  layer.  This  procedure  was 
continued  until  all  of  the  longitudinal  layer  had  been  removed.  Then 
strips  of  isolated  circular  muscle  or  longitudinal  muscle  were  pinned 
to  the  bottom  of  the  constant  temperature  bath  for  study  with  extra¬ 
cellular  electrodes  as  illustrated  in  Figure  4.5. 

Seven  experiments  using  ten  different  specimens  were  conducted 
on  isolated  longitudinal  and  circular  muscle.  From  each  experiment 
histology  specimens  were  prepared  to  confirm  that  the  separation  was 
in  the  plane  between  the  circular  and  longitudinal  layers.  For  this 
purpose  portions  of  the  intact  wall  and  isolated  layers  were  placed 
immediately  into  gl  ucoal  dehyde  and  later  mounted  and  stained  using 
hematoxyl i n-eosi n  or  trichrome  stains.  Figure  4.6  shows  three  slides 
typical  of  the  results  obtained  using  hematoxyl i n-eosi n  stain.  Figure 
4.6a  shows  the  intact  muscle  wall  with  the  longitudinal  layer  near 
the  left  margin  of  the  slide.  The  longitudinal  muscle  cells  lie 
approximately  parallel  to  the  plane  of  the  picture.  Figures  4.6b 
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Figure  4.6  Histology  slides  showing  (a)  the  intact  muscle  wall, 
(b)  the  isolated  lonitudinal  layer,  and  (c)  the 
isolated  circular  layer.  Magnification:  x62.5. 
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and  4.6c  show  the  longitudinal  and  circular  layers  respectively  after 
separation.  For  the  majority  of  histology  specimens  the  separation 
occurred  in  the  plane  between  the  two  layers. 

4.3.6  Dual  Chamber  Studies 

The  dual  chamber  studies  were  conducted  to  further  investigate  the 
role  of  each  of  the  muscle  layers  in  the  generation  of  slow  wave  activity. 
These  experiments  were  based  on  the  simultaneous  recording  of  the 
amplitude  and  phase  relationships  of  the  slow  wave  signals  on  both 
sides  of  the  intact  muscle  wall  and  isolated  layers. 

The  dual  chamber  constant  temperature  bath  shown  in  Figure  4.7 
was  designed  to  facilitate  this  set  of  experiments.  This  bath  consisted 
of  two  electrically  isolated  chambers,  each  separately  supplied  with 
oxygenated  Krebs-Ringer  solution.  Both  chambers  were  maintained  at 
the  same  constant  temperature  by  a  surrounding  electrically  insulated 
water  bath. 

A  two  centimeter  square  section  of  colonic  muscle  was  prepared  in 
the  manner  previously  described  in  section  4.1  and  mounted  between  the 
two  chambers.  A  1  cm  diameter  window  permitted  contact  o^  opposite 
sides  of  the  tissue  with  the  Krebs-Ringer  solution  of  the  chambers. 
Petroleum  jelly  was  used  to  insure  good  electrical  insulation  in  all 
areas  except  for  the  1  cm  windows.  Before  each  experiment  the  impedance 
between  the  chambers  was  checked  to  ensure  the  absence  of  leaks  creating 
a  short  circuit  condition. 

A  total  of  twelve  experiments  were  conducted.  Six  of 
these  experiments  were  conducted  on  specimens  of  intact  muscle 
wall  containing  both  circular  and  longitudinal 
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Figure  4.7  The  dual  chambered  measurement  system. 


57 


muscle  layers.  The  other  six  experiments  were  conducted  on  isolated 
longitudinal  and  circular  muscle. 

4.4  INTRACELLULAR  STUDIES 

4.4.1  Introduction 

Extracellular  slow  wave  measurements  would  be  of  maximum  benefit 
if  they  could  be  directly  related  to  variations  in  i ntracel 1 ul ar 
potentials.  This  relationship  was  uncertain  because  information  was 
not  available  on  the  intracellular  potentials  of  canine  colonic  smooth 
muscle.  Therefore  it  was  necessary  to  conduct  a  number  of  microelec¬ 
trode  intracellular  studies.  These  studies  were  also  required  to 
investigate  the  transduction  mechanism  responsible  for  the  generation 
of  the  extracel 1 ul ar  signal. 

Mi croel ectrode  measurements  on  smooth  muscle  are  difficult  to 
perform  because  of  the  small  size  of  smooth  muscle  cells  and  because 
the  muscle  is  continually  active.  For  successful  penetration  of  this 
type  of  cell  the  microel ectrode  must  be  less  than  one  micron  in  diameter. 
Such  electrodes  are  very  fragile  and  are  easily  broken  by  muscle  move¬ 
ment  or  system  vibration. 

A  measurement  system  was  designed  which  minimized  system  vibration 
and  electrode  breakage.  Twenty-four  i ntracel 1 ul ar  experiments  were 
conducted  on  thirty-eight  different  specimens  from  the  transverse  and 
proximal  colon.  Eight  of  these  studies  included  measurements  on  the 
intact  muscle  wall,  five  included  studies  of  isolated  longitudinal  and 
circular  muscle  layers,  eight  were  studies  of  the  i ntracel 1 ul ar-to- 
extracel 1 ul ar  transduction  mechanism,  and  three  recorded  extracellular 
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and  i ntracel 1 ul ar  signals  simultaneously.  Basal  measurements  of  intra¬ 
cellular  potentials  from  the  mucosal  side  of  the  intact  muscle  wall  were 
obtained  for  all  but  five  of  the  experiments. 

4.4.2  The  Measurement  System 

The  intracellular  measurement  system  is  shown  in  Figure  4.8.  It 
consisted  of  a  vibration  free  table,  the  constant  temperature  tissue 
bath  originally  used  for  extracel 1 ular  studies,  a  microelectrode  manipulato 
and  holder,  a  hydraulic  manipulation  system  for  precise  control  of 
electrode  movement,  and  a  battery  powered  el ectrometer-i nput  amplifier. 

The  output  from  the  amplifier  was  simultaneously  recorded  on  a  storage 
oscilloscope,  on  an  fm  tape  recorder,  and  on  a  paper  chart  recorder. 

The  complete  system  with  the  exception  of  the  recording  devices  and  the 
hydraulic  control  micrometer  was  housed  inside  a  Faraday  cage. 

The  vibration- free  bench  consisted  of  a  horizontal  concrete  platform 
(4  ft  x  6  ft  x  4  inches)  supported  by  two  legs  made  from  concrete  blocks. 
The  Faraday  cage  was  placed  on  this  platform.  Inside  the  cage  and  of 
slightly  smaller  dimensions  was  a  second  concrete  platform  supported 
by  four  air-filled  rubber  tennis  balls  and  two  balls  of  steel-wool. 

The  rubber  balls  provided  decoupling  from  building  vibrations  and  the 
steel-wool  balls  were  adjusted  to  provide  suitable  damping. 

The  mi croel ectrode  manipulator  and  holder  was  of  the  standard 
type  designed  for  mi croel ectrode  measurements  but  a  hydraulic  control 
system  was  added  to  isolate  the  experimenter  from  the  measurement  system 
and  to  improve  the  precision  of  the  microelectrode  manipulation  system. 

The  details  of  the  measurement  system  are  shown  in  Figure  4.9.  The 
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FARADAY  SHIELD 


Figure  4.8  The  intracellular  measurement  system. 
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hydraulic  system  consisted  of  two  oil  filled  precision  syringes  coupled 
by  means  of  an  oil  filled  polyethylene  tube.  The  ratio  of  diameters 
for  the  two  syringes  was  25.  The  smaller  of  the  two  syringes  was 
driven  by  a  micrometer  system  which  lay  outside  the  Faraday  cage.  The 
control  micrometer  permitted  a  slow  precise  movement  of  the  microelectrode 
with  one  unit  movement  of  the  micrometer  corresponding  to  an  electrode 
movement  of  0.1  micron. 

The  microel ectrode  amplifier  was  a  Medistor  Model  A-35L  electro- 

g 

meter-input  amplifier  with  an  input  impedance  of  5  x  1  0  ohms  and  an 
RMS  noise  output  of  20  y  V . 

The  electrodes  were  pulled  from  1.2  mm  borosilicate  glass  tubing 
and  they  were  filled  with  3  M  KC1 .  It  was  determined  that  an  electrode 
impedance  of  between  50  and  100  megohms  was  optimum  for  this  work. 

As  shown  by  Firth  and  Defilice  (1971)  the  tip  diameter  of  such  an 
electrode  is  in  the  range  from  0.1  to  0.5  microns. 

4.4.3  Circular  Muscle  Measurements 

Ten  experiments  were  conducted  on  the  mucosal  side  of  intact  muscle 
strips  (with  the  mucosa  removed)  to  determine  the  best  procedures  for 
obtaining  intracellular  measurements  and  to  obtain  statistically  sig¬ 
nificant  information  on  the  basal  potentials  in  the  circular  muscle 
1 ayer. 

The  tissue  was  prepared  as  described  in  section  4.1.  However  the 
final  specimen  was  trimmed  to  the  much  smaller  size  of  0.2  cm  by  1.0  cm 
to  minimize  problems  caused  by  muscle  contraction.  This  was  done  such 
that  the  specimen  sides  were  parallel  to  the  circular  and  longitudinal 
directions.  The  way  in  which  a  specimen  is  mounted  in  the  tissue  bath 
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is  of  critical  importance  for  the  success  of  an  experiment.  The  optimum 
mounting  arrangement  is  shown  in  Figure  4.10(a).  The  tissue  is  pinned 
in  four  locations  such  that  it  is  not  under  tension  and  so  that  a  line 
joining  the  two  closest  pins  is  parallel  to  the  short  dimension  of  the 
specimen.  The  optimum  measurement  location  is  then  at  the  midpoint 
between  these  two  pins.  Care  must  be  taken  to  ensure  that  the  tissue 
is  not  pinched  or  distorted  from  its  normal  shape  by  the  pins. 

In  addition  to  the  above  ten  experiments,  basal  recordings  on  the 
mucosal  side  of  the  intact  muscle  wall  were  obtained  from  thirteen  other 
experiments . 

4.4.4  Longitudinal  Muscle  Measurements 

The  measurement  of  i ntracel 1 ul ar  potentials  in  the  longitudinal 
layer  proved  to  be  a  much  more  difficult  task  than  the  measurement 
of  such  potentials  in  the  circular  muscle  layer.  The  submucosal  layer 
provided  easy  access  to  the  circular  muscle  cells  without  problems  due 
to  electrode  breakage.  However,  the  serosa,  inspite  of  the  fact  that 
it  is  a  very  thin  layer,  proved  to  be  a  major  obstacle  preventing  access 
to  longitudinal  cells.  Unsuccessful  attempts  were  made  to  remove  the 
serosa  without  causing  significant  damage  to  the  longitudinal  muscle. 

The  technique  which  proved  most  successful  utilized  very  thin  cross- 
sections  of  colon  wall  cut  parallel  to  the  longitudinal  direction. 

These  sections  were  approximately  0.1  cm  thick  and  were  mounted  as 
shown  in  Figure  4.10(b). 
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Figure  4.10  Tissue  mounting  configurations. 
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4.5  STUDIES  OF  THE  INTRACELLULAR-TO-EXTRACELLULAR  TRANSDUCTION  PROCESS 

4.5.1  Introduction 

As  reported  in  Chapter  3,  the  generation  of  an  "injury  potential" 
has  been  given  as  the  reason  for  the  existence  of  the  the  extracellular 
pressure  electrode  slow  wave  signals.  However  no  research  has  been 
published  which  proves  that  this  is  correct.  The  work  done  by  Gillespie 
(1962)  produced  evidence  which  led  Bortoff  (1967)  and  others  to  believe 
an  "injury  potential"  was  responsible  for  the  extracel  1  ul ar  signal. 
However,  Gillespie's  work  related  to  large  mi croel ectrodes  and  individual 
smooth  muscle  cells.  The  pressure  electrodes  used  in  this  present 
work  are  many  times  larger  than  the  smooth  muscle  cells  responsible  for 
the  signal  generation.  In  fact,  the  2  mm  diameter  electrode  has  of 
the  order  of  five  to  ten  million  cells  directly  under  it.  One  cannot 
assume  the  theory  developed  for  individual  cells  applies  in  this 
situation.  Consequently  eleven  experiments  were  conducted  to  inves¬ 
tigate  the  relationship  between  the  extracel 1 ul ar  signals  and  the  intra¬ 
cellular  potentials  of  cells  directly  under  the  pressure  electrodes. 

The  objective  of  this. set  of  experiments  was  twofold:  firstly, 
to  demonstrate  that  the  extracellular  slow  wave  signal  is  a  function 
of  intracellular  variations  of  the  same  shape  and  frequency;  and 
secondly,  to  show  that  an  injury  potential  of  the  type  described  by 
Gillespie  (1962)  exists  in  cells  under  pressure  electrodes. 

4.5.2  A  Comparison  of  Intracel 1 ul ar  and  Extracellular  Electrode 

Potential  s 


Three  experiments  were  conducted  to  simultaneously  record  and 
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compare  extracellular  and  intracellular  potentials.  This  was  difficult 
to  accomplish  because  the  extracellular  electrode,  in  deforming  the 
tissue,  created  an  increased  contractile  activity  around  the  extra¬ 
cellular  measurement  site.  This  tissue  movement  greatly  reduced  the 
number  of  successful  penetrations.  Also  the  smaller  tissue  specimens 
required  by  the  i ntracel 1 ul ar  measurements  were  not  sufficiently  large, 
to  completely  surround  the  extracellular  electrode,  as  in  the  extra¬ 
cellular  measurements  of  section  4.3.  Consequently  the  extracel 1 ul ar 
signals  obtained  were  smaller  in  amplitude,  and  the  waveform  not  as 
constant  from  cycle  to  cycle.  The  extracellular  and  intracellular 
electrode  holders  prevented  the  two  measurement  sites  from  being  closer 
than  5  mm. 

4.5.3  Injury  Potential  Measurement 

Eight  experiments  were  conducted  to  determine  if  an  "injury 
potential"  was  associated  with  the  pressure  electrode  recordings.  The 
apparatus  shown  in  Figure  4.11  was  designed  to  simulate  the  pressure 
electrode.  It  was  constructed  from  the  2  mm  diameter  glass  used  to 
construct  the  extracellular  electrodes.  The  glass  was  cemented  to  a 
metal  holder  to  which  a  force  was  applied.  This  simulated  electrode 
was  first  applied  to  the  mucosal  side  of  the  tissue  as  shown  in 
Figure  4.11  with  a  known  force  of  1  to  5  grams  (.01  to  .05  newtons). 
Then  the  microelectrode  was  positioned  in  the  center  of  tne  cylinaer 
and  gradually  advanced  until  cell  penetration  was  obtained. 
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Figure  4.11  Apparatus  simulating  the  application  of  the  pressure 
electrode  to  an  intact  muscle  specimen. 
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4.6  COMPUTER  ANALYSIS 

The  data  from  all  experiments  were  digitized  and  processed  on  a 
Hewlett-Packard  21  MX/E  computer  system.  This  system  shown  in  Figure 
4.12  consisted  of  the  21  MX/E  computer,  a  5  megabyte  disc  drive,  an 
IBM  compatible  digital  magnetic  tape  unit,  an  eight  channel  analog- 
to-digital  converter,  a  three  channel  digital-to-analog  converter 
connected  to  a  CRT  display,  a  high  speed  paper  tape  reader,  and  a  high 
speed  printer  terminal.  Programs  were  available  for  high  speed  analog- 
to-digital  conversion  and  descrete  Fast  Fourier  Transformation  (Cooley 
&  Tookey,  1969)*.  Additional  programs  were  written  to  perform  the 
following  functions: 

1.  To  transfer  data  files  to  and  from  the  computer  and  its 
storage  and  input/output  devices. 

2.  To  create  and  apply  digital  filters. 

3.  To  do  auto  and  cross-correl ations  of  selected  data  channels. 

4.  To  determine  the  statistical  parameters  of  cycl e-to-cycl e 
frequencies  and  phase  angles. 

5.  To  determine  the  coherence  function  for  selected  channels. 

6.  To  display  and  plot  selected  data  files. 

The  programs,  descriptions,  and  examples  of  typical  inputs  and 
outputs  are  contained  in  Appendix  A. 


*Both  programs  written  by  Dr.  Z.  Koles,  University  of  Alberta. 
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Figure  4.12  The  Hewlett-Packard  21  MX/E  computer  system. 
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CHAPTER  5 


EXPERIMENTAL  RESULTS 


5.1  INTRODUCTION 

In  this  chapter  the  results  obtained  from  the  extracellular  and 
intracel 1 ul ar  experiments  described  in  Chapter  4  are  presented.  The 
organization  of  this  chapter  parallels  that  of  Chapter  4.  The  results 
from  the  electrode  studies  are  presented  first.  Following  this  are  th 
outcomes  of  the  extracellular  pressure  electrode  experiments.  These 
results  are  divided  into  two  parts,  the  results  providing  information 
on  coupled  oscillator  characteristics  and  the  results  which  provide 
evidence  relating  to  the  origin  of  the  slow  wave  and  the  role  played 
by  the  individual  muscle  layers  in  the  generation  of  this  activity. 

The  reader  is  referred  again  to  Figures  4.3  and  4.4  to  see  how  the 
extracel 1 ul ar  experiments  and  results  have  been  organized.  The  extra¬ 
cellular  information  is  followed  by  the  i ntracel 1 ul ar  results  and  then 
by  the  results  of  experiments  investigating  the  relationship  between 
extracel 1 ul ar  and  intracellular  signals. 

5.2  ELECTRODE  CHARACTERISTICS 

From  the  impedance  versus  frequency  measurements  on  the  silver- 
si  1  verchl  oride  pressure  electrodes  it  was  found  that  the  impedance 
magnitude  increases  inversely  with  decreasing  frequency  down  to  0.01 
Hz.  The  resistance  and  reactance  both  vary  as  an  inverse  power  of 
frequency  and  are  close  to  being  equal  at  all  frequencies.  The  real 
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component  of  the  impedance  consists  of  two  parts,  a  constant  component 
and  a  component  which  varies  with  frequency.  Figure  5.1  is  a  plot  of 
the  real  and  reactive  components  of  a  typical  electrode.  The  equations 
describing  the  electrode  resistance  and  reactance  are: 

R  =  400  +  1820  f'0'36  ohms  5.1 

X  =  1200  f  ®’3®  ohms  5.2 

The  impedance  characteri sti cs  of  the  electrodes  changed  continually 
from  the  time  of  chloriding.  During  the  first  24  hours  after  chloriding 
the  impedance  characteristics  were  not  measured.  However,  after  this 
period  the  electrode  impedance  increased  daily  and  changed  such  that 
it  no  longer  varied  as  a  constant  power  of  frequency.  A  plot  of  the 
real  and  reactive  components  of  the  impedance  of  a  twelve-day-old 
electrode  is  shown  in  Figure  5.2. 

The  amplitude  of  noise  produced  by  these  electrodes  was  very  small. 
Using  a  high  gain  amplifier  and  a  Beckman  polygraph  recorder,  the  noise 
amplitude  was  found  to  be  of  the  order  of  0.2  yV  peak-to-peak .  Figure 
5.3  is  a  recording  of  typical  electrode  noise  compared  to  the  noise 
of  an  unchlorided  silver  electrode. 

5.2.1  An  Electrode  Model 

The  electrode  impedance  may  be  modelled  using  a  modification  of 
the  model  proposed  by  Warburg  (1899).  The  model  is  shown  in  Figure  5.4. 
This  modification  is  necessary  in  order  that  the  model  better  represent 
an  electrode  characteri sti c  as  demonstrated  by  our  experiments,  that 
is  the  bulk  resistance  of  the  electrode  consists  of  two  components. 
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Figure  5.1  The  real  and  reactive  components  of  the  impedance  of  a 
typical  electrode  plotted  as  a  function  of  frequency. 
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Figure  5.2  Real  and  reactive  components  of  the  impedance  of  a  twelve- 
day-old  electrode. 
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Figure  5.3  Electrode  Noise 

(a)  Normal  Electrode 

(b)  Unchlorided  Electrode 
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One  component  Rgp  is  the  bulk  resistance  of  the  fluid  surrounding  the 
electrode  tip,  and  the  other  R^^  is  the  bulk  resistance  of  the  silver- 
chloride  layer.  Measurement  of  the  electrode  impedance  at  high  frequency 
before  and  after  chloriding  has  yielded  average  values  of  R^p  and  R^^ 
of  300  and  100  ohms  respecti vely .  Also  from  high  frequency  measurements 
the  double  layer  incremental  capacitance  has  an  average  value  of 

_3 

15  x  10  microfarads.  The  charge  transfer  resistance  of  the  electrode 
can  be  calculated.  From  Cobbold  (1974),  R  is  defined  as  follows: 


R  =  -51 

Kt  zFI 


5.3 


o 


where  R  is  the  gas  constant,  T  is  the  temperature  in  °K,  z  is 
the  valence,  F  is  Faraday's  constant, 
and  I  is  the  exchange  current  of  the  electrode 

If  T  =  300°K  and  the  exchange  current  is  the  same  as  reported  for 

Ag/Ag+  in  0.1  Normal  saline  (Cobbold,  1974);  the  calculated  value  of 

R  is  0.72  ohms. 

From  plots  such  as  Figure  5.1,  R^  and  can  be  calculated.  For 
the  electrode  of  Figure  5.1 

R  ,  =  1029  f-0*36  ohms  5-4 

d 

and  C.  =  102.1  x  10"6  f-0’64  farads.  5.5 
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5.2.2  Consideration  of  Signal  Distortion  Due  to  Electrode  Characteri sti cs 

One  reason  for  studying  the  electrode  impedance  was  the  concern 
that  the  electrode  might  create  signal  distortion  due  to  the  very  low 
frequency  components  of  the  signal  and  because  si  1 ver-si  1  verchl ori de 
electrode  impedances  greatly  increase  with  decreasing  frequency.  To 
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determine  the  significance  of  signal  distortion  caused  by  the  electrode, 

the  electrode  model  is  considered  in  the  slow  wave  measurement  circuit 

of  Figure  5.5.  The  equivalent  source  resistance  R  of  the  slow  wave 

eg 

was  measured  and  found  to  be  in  the  range  from  20  to  to  40  to.  In  the 

circuit  a  mean  value  of  30  to  is  used.  Using  the  values  for  R.  and 

d 

0^  qiven  in  the  previous  section,  the  equivalent  impedance  Z  of  the 
electrode  is  defined  by  the  following  equation: 


z  =  400.72  +  1  029f  +  j ( 0 . 4 8  +  27.1  x  10"6f  +  38.8  x  10~3f°-64) 

eg  1  -  45.2  x  10"9f  +  j ( 67 . 9  x  10"9f  +  97.0  x  1 0“6f°* 


The  slow  wave  signal  will  have  no  significant  frequency  components 
above  100  Hz.  Thus,  when  considering  the  effect  of  the  electrode 
impedance  on  the  slow  wave  signal,  Z  may  be  approximated  by  the 


following  equation: 

Z 

eg 


400.72  +  1  02 9 f  +  j(0.48  +  38.3  x  1 0_3f° * 64 ) 

1  +  j0 


5.7 


From  the  circuit  of  Figure  5.5,  the  potential  v  at  the  input  to  the 
amplifier  is  given  by 


v  = 


2x10 


2.03  x  10  +  Z 


sw 


5.8 


eg 


where  the  input  impedance  R.  of  the  amplifier  is  2  Mft.  Substituting 
the  approximate  value  of  Z^  from  equation  5.7  into  equation  5.8  gives 


2.0  x  10 


2.03  x  106  +  400  +  1 029f  +  j ( 0 . 48  +  38.3  x  10"3f°'64) 


5.9 


For  frequencies  below  100  Hz  the  imaginary  component  of  this  expression 
is  not  significant  and  the  phase  shift  produced  by  the  electrode  is 


- 


. 
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Figure  5.5  Slow  wave  measurement  circuit. 
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very  close  to  zero  degrees.  Therefore,  for  frequency  components  from 

zero  to  100  Hz,  the  measured  voltage  v  ranges  from  0.94  v  to  0.99 

sw 

vsw  with  no  si9nifican1:  distortion  due  to  phase  shift.  Thus,  if  the 
pressure  electrode  theory  described  in  Chapter  3  is  correct,  the  intra¬ 
cellular  slow  wave  will  be  faithfully  reproduced  with  only  amplitude 
seal ing. 

5.3  EXTRACELLULAR  RESULTS 

5.3.1  Phase  Angles  as  a  Function  of  Electrode  Spacing 

As  described  in  Chapter  4,  section  4.3.2,  six  experiments  were 
conducted  with  from  two  to  four  electrodes  applied  to  one  side  of 
the  intact  muscle  wall  with  the  objective  of  studying  the  frequency 
and  phase  relationships  of  the  oscillating  regions  of  the  colon. 

In  addition,  multisite  recordings  were  obtained  from  seventeen 
experiments  conducted  for  other  purposes.  Figure  5.6  shows  a  typical 
recording  from  four  electrodes  applied  to  the  mucosal  side  of  the 
intact  muscle  wall.  The  electrodes  were  oriented  in  the  circular 
direction  with  a  spacing  of  2.5  mm.  Similar  recordings  were  obtained 
from  the  serosal  side  of  the  intact  wall. 

The  data  from  the  above  twenty-three  experiments  were  digitized 
and  processed  using  fast  fourier  transform,  zero  crossing,  and  cross¬ 
correlation  programs.  For  the  recordings  of  Figure  5.6  the  output  of 
the  fast  fourier  transform  program  "PSPT"  is  shown  in  Figure  5.7. 
Normalized  signal  power  is  plotted  versus  frequency  in  cycles/minute. 
Figure  5.8  is  the  print-out  of  the  program  "PHSE"  which  utilizes  zero 


r ,  .  r 
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Figure  5.6  A  recording  from  four  pressure  electrodes  applied  to  the 
mucosal  side  of  the  intact  muscle  wall.  The  electrodes 
are  oriented  in  the  circular  direction  with  a  spacing  of 
2.5  mm . 
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82 


crossing  techniques  to  compute  cycl e-by-cycl e  frequencies  and  relative 
phase  angles  of  two  selected  channels.  The  frequency  of  each  channel 
and  the  relative  phase  angle  in  degrees  is  plotted  versus  the  cycle 
number.  Figure  5.9  is  the  output  of  the  cross  correlation  program 
"XCORR"  when  applied  to  the  recordings  of  Figure  5.6.  The  plot  gives 
the  normalized  cross-correlation  function  versus  the  delay  time  in 
seconds.  Program  listings,  flow  charts,  and  typical  input/output  data 
for  each  of  these  programs  is  given  in  Appendix  A. 

As  indicated  in  section  4.3.2  of  Chapter  4,  only  six  of  the  nine 
experiments  produced  useful  data. 

The  data  from  the  above  experiments  are  summarized  in  Table  5.1. 
Each  entry  in  the  table  is  an  average  based  on  a  cycl e-by-cycl e  analysis 
of  a  ten  to  fifteen  minute  recording  analyzed  visually  as  well  as  with 
the  zero-crossing  program  "PHSE"  and  the  cross-correl ation  program 
"XCORR". 

5.3.2  Interpretation  of  Results 

The  following  two  observations  can  be  deducted  from  Table  5.1: 

1.  There  is  great  variability  of  phase  angles  for  each  of  the 
electrode  spacings. 

2.  The  average  phase  angle  in  the  circular  direction  is  less  than 
that  in  the  longitudinal  direction. 

These  two  points  require  further  clarification.  In  Chapter  3 
it  was  postulated  that  the  extracel 1 ul ar  slow  wave  signal  is  produced 
by  the  distortion  of  cells  under  and  around  the  periphery  of  the 
pressure  electrode.  The  area  encompassed  by  these  cells  is  not  known 
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Table  5.1  Phase  angle  between  recording  sites 
as  a  function  of  electrode  spacing. 
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as  a  function  of  electrode  spacing  (cont. 
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and  will  be  a  function  of  a  particular  tissue  and  of  the  force  with 
which  the  electrode  is  applied.  Because  of  this  fact,  the  distance 
between  regions  from  which  different  slow  waves  originate  is  unknown. 

For  the  2  mm  diameter  electrodes  with  center-to-center  spacing  of  2.5 
mm,  the  spacing  between  oscillating  regions  could  range  from  less  than 
0.5  mm  to  more  than  4.5  mm.  Also,  more  than  two  oscillating  regions 
could  exist  between  the  electrodes. 

The  above  situation,  plus  the  fact  that  the  visual  analysis  and 
the  zero-crossing  program  measure  angles  with  a  modulus  of  360  degrees, 
make  it  impossible  to  be  certain  about  the  phase  angles  existing  between 
oscillating  regions. 

However  some  insight  may  be  gained  from  these  measurements 
regarding  the  probable  order  of  magnitude  of  both  the  relative  phase 
angle  and  the  spacing  between  regions.  If  the  spacing  of  oscillating 
regions  in  any  particular  direction  is  much  smaller  than  the  effective 
electrode  spacing,  more  than  two  regions  exist  between  the  electrodes. 

The  resulting  phase  angle  will  be  a  product  of  the  average  phase  angle 
between  regions  and  the  number  of  regions  less  one.  If  this  product 
is  180  degrees  or  greater,  doubling  the  electrode  spacing  will  not 
necessarily  result  in  a  doubling  of  the  measured  angle,  because  the 
angle  is  determined  with  respect  to  the  modulus  360  degrees.  However, 
if  the  spacing  between  regions  is  of  the  same  size  as  the  effective 
electrode  spacing,  the  phase  angle  would  increase  in  proportion  to 
the  increase  in  electrode  spacing. 

In  the  circular  direction  the  average  phase  angles  for  2.5  mm  and 
5.0  mm  electrode  spacings  are  29.0  ±  23.2  and  68.3  ±  48.0  degrees 
respectively.  A  doubling  of  the  electrode  spacing  results  in  the  average 
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phase  angle  changing  by  a  factor  of  2.36.  This  is  evidence  that  the 
average  phase  angle  per  millimeter  is  of  the  order  to  11.60  degrees  in 
the  circular  direction.  However,  in  the  longitudinal  direction,  doubling 
the  electrode  spacing  results  in  the  mean  phase  angle  changing  from 

117.3  -  51.0  to  168.0  ±  79.4  degrees.  This  increase  of  only  43  percent 
indicates  either  (1)  a  very  large  phase  lag  between  oscillating  regions, 

(2)  a  spacing  between  oscillating  regions  much  smaller  than  the  electrode 
spacing,  or  (3)  a  combination  of  both  these  factors. 

5.3.3  The  Results  from  the  Cutting  Experiments 

The  results  from  the  twenty  cutting  experiments  may  be  divided 
into  two  parts  (1)  the  results  from  longitudinally  oriented  muscle 
specimens  and  (2)  the  results  from  the  circularly  oriented  specimens. 

The  results  from  the  longitudinally  oriented  muscle  strips  are  summarized 
in  Table  5.2.  The  following  facts  can  be  obtained  from  this  table: 

1.  The  average  frequency  of  the  slow  waves  before  cutting  was  5.6  ± 

1.2  cycles/minute  and  there  was  no  significant  change  in  frequencies 
at  recording  sites  on  either  side  of  the  cut.  (Paired  Student  t 
test,  p  <  . 06 . ) 

2.  For  an  electrode  spacing  of  2.5  mm  the  average  phase  angles  between  the 
sites  to  be  decoupled  was  121.7  ±  49.8  degrees.  After  cutting  the  phase 
angle  continually  varied  and  neither  visual  analysis  nor  program  "PHSE" 
could  determine  a  meaningful  average  phase  angle  in  80  percent  of  the 
results.  In  these  cases  "INC"  meaning  "not  coupled"  is  shown  in  the  table. 

3.  The  average  coherence  before  cutting  is  0.60  ±  0.24,  and  the  average 
after  cutting  is  0.23  ±  0.21.  There  is  a  significant  drop  of  61 
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percent  in  the  coherance  figure  as  a  result  of  the  decoupling  caused 
by  cutting. 

The  results  may  be  interpeted  as  follows: 

1.  There  is  no  frequency  gradient  in  the  longitudinal  direction  since 
the  frequencies  are  not  changed  by  the  cutting. 

2.  The  average  phase  angle  is  121.7 ±  49.8  degrees.  This  supports  the 
figure  of  117.3  ±  51.0  degrees  computed  previously.  Prior  to  these 
experiments  it  was  considered  a  possibility  that  a  phase  angle  might 
exist  because  of  measurement  technique  and  not  due  to  coupling 
between  regions.  This  hypothesis  is  rejected  by  virtue  of  the 

fact  that  an  average  phase  angle  could  not  be  computed  in  80  percent 
of  the  observations  after  cutting. 

3.  The  results  show  that  the  average  coherence  may  be  used  as  a  gross 
indicator  of  the  degree  of  coupling.  They  do  not  indicate  a  means 
of  relating  coherence  magnitude  directly  to  the  degree  of  coupling. 
However,  Table  5.2  shows  that  the  amount  by  which  this  factor  is 
less  than  one  is  an  indicator  of  the  statistical  variation  which 
will  exist  in  the  relative  phase  angle. 

The  data  from  strips  oriented  in  the  circular  direction  are  sum¬ 
marized  in  Table  5.3.  These  results  are  divided  into  two  sections, 
for  experiments  involving  cuts  on  the  mesenteric  border  and  for 
experiments  involving  cuts  1  cm  from  the  mesentery. 

The  results  are  summarized  as  follows: 

1.  The  mean  frequency  before  cutting  was  7.75  -  1.2  cycles/minute 
and  the  mean  phase  angle  was  70.2  ±  52.4  degrees.  After  cutting 
either  on  or  off  the  mesentery  no  meaningful  average  angle  could 
be  computed.  This  situation  is  indicated  as  "NC"  in  the  table. 


, 


Table  5.3  Cutting  experiment  results  for  circularly  oriented  strips. 
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2.  After  cutting  on  the  mesentary  there  is  a  drop  in  frequency  at  all 
recording  sites  but  not  a  significant  difference  in  frequency 
between  sites. 

3.  After  cutting  1  cm  from  the  mesentery  there  is  a  drop  in  frequency 
at  all  recording  sites,  but  the  drop  at  sites  decoupled  from  the 
mesentery  because  of  the  cut  is  much  greater. 

4.  In  all  cases  the  coherence  value  dropped  significantly  as  a  result 
of  cutti ng . 

These  results  may  be  interpreted  as  follows: 

1.  Coupling  is  again  shown  to  be  better  in  the  circular  direction  than 
in  the  longitudinal  direction,  as  indicated  by  the  much  smaller 
average  phase  angle. 

2.  There  is  a  frequency  gradient  in  the  circular  direction  in  both 
directions  away  from  the  mesentery.  This  is  indicated  by  the  drop 
in  frequency  of  regions  decoupled  by  cutting  on  either  side  of 
the  mesentery. 

3.  The  mesenteric  border  plays  a  role  in  determining  slow  wave  frequency. 
This  is  indicated  by  the  above  as  well  as  by  the  fact  that  cutting 

on  the  mesentery  resulted  in  a  drop  in  frequency  at  all  recording 
si tes . 

4.  Again  it  is  apparent  that  the  coherence  figure  is  a  good  indicator 
of  the  statistical  variation  in  relative  phase  angle,  and  conse¬ 
quently  an  indicator  of  the  degree  of  coupling. 

5.3.4  Results  from  the  Studies  of  Isolated  Muscle  Layers 

With  the  development  of  a  suitable  dissection  procedure  it  was 

possible  to  study  isolated  longitudinal  and  isolated  circular  muscle. 


. 
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The  electrical  activity  from  isolated  longitudinal  muscle  was  recorded 
in  four  experiments.  Slow  wave  activity  was  not  observed  in  any  of  these 
experiments.  Figure  5.10  is  a  typical  record  from  four  extracellular 
electrodes  applied  to  the  serosal  side  of  the  isolated  longitudinal 
1 ayer . 

Seven  extracel 1 ular  experiments  were  conducted  on  the  isolated 
circular  muscle  layer.  In  five  out  of  seven  experiments  slow  wave 
activity  was  present  on  both  sides  of  the  isolated  circular  muscle  layer. 
In  two  of  the  seven,  when  the  dissection  had  been  particularly  difficult, 
no  slow  wave  was  observed.  In  each  experiment  the  signal  obtained  from 
the  longitudinal  side  of  the  circular  muscle  was  lower  in  amplitude  and 
more  irregular  in  shape  from  cycle  to  cycle.  It  is  presumed  that  this 
was  due  to  cell  damage  caused  by  the  dissection  process.  Figure  5.11(a) 
is  a  typical  record  obtained  from  four  electrodes  applied  to  the 
mucosal  side  of  the  isolated  circular  muscle  layer  while  Figure  5.11(b) 
is  a  record  obtained  from  the  serosal  side  of  that  same  specimen. 

5.3.5  The  Submucosal  Layer 

As  described  in  Chapter  4,  in  preparation  for  the  in  vitro  measure¬ 
ments  the  mucosa  and  muscularis  mucosa  are  removed.  Upon  the  removal 
of  these  two  layers  a  thin  layer  of  submucosa  remains  attached  to  the 
intact  muscle  wall.  It  was  discovered  that  this  thin  layer  of  sub¬ 
mucosa  could  be  easily  and  uniformly  peeled  from  the  circular  muscle 
layer.  In  fact,  the  removal  of  this  layer  is  accomplished  much  more 
easily  than  the  separation  of  longitudinal  and  circular  muscle  layers. 
However,  what  was  most  unexpected  was  that  slow  wave  activity  in  the 
intact  wall  and  in  isolated  circular  muscle  ceased  when  this  layer 
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Figure  5.10  Four  pressure  electrode  signals  from  the  serosal  side  of 
isolated  longitudinal  muscle. 
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Figure  5.11(a)  Slow  waves  from  the  mucosal  side  of  isolated  circular 

muscl e . 
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Figure  5.11(b)  Slow  waves  from  the  longitudinal  side  of  isolated 

circular  muscle. 
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was  removed. 

In  five  experiments  using  twelve  different  specimens  the  normal 
slow  wave  was  abolished  in  all  cases  by  complete  removal  of  the  sub¬ 
mucosal  layer.  For  only  two  of  the  specimens  was  there  anything  which 
could  be  interpreted  as  cyclic  electrical  activity. 

Histology  specimens  were  prepared  for  all  of  these  experiments. 
Figure  5.12  is  typical  of  the  histology  results  showing  the  circular 
muscle  layer  before  and  after  the  removal  of  the  submucosal  layer  and 
also  showing  the  isolated  submucosal  layer.  All  three  slides  have 
been  treated  with  Trichrome  stain.  As  can  be  seen  the  submucosal  layer 
consists  mainly  of  colagen,  appearing  as  a  green  color  on  the  slides. 
Also  a  small  layer  of  cells  which  look  similar  to  circular  muscle  cells 
is  present  in  the  submucosal  layer.  All  attempts  to  differentiate 
these  cells  using  different  histology  stains  have  been  unsuccessful. 

In  the  above  five  experiments  normal  slow  wave  activity  was  first 
recorded  on  the  intact  specimen.  Then  the  submucosal  layer  was  removed 
from  half  of  the  circular  muscle.  Recordings  were  made  from  the  intact 
half  of  the  tissue  as  well  as  from  the  area  where  the  submucosa  had 
been  removed.  Normal  slow  wave  activity  was  recorded  from  the  intact 
half.  However,  slow  wave  activity  was  found  in  the  other  half  as  well 
but  with  amplitude  diminishing  as  a  function  of  distance  from  the  intact 
portion  of  the  tissue.  Figure  5.13  shows  a  typical  recording  from  four 
electrodes  before  and  after  the  removal  of  half  of  the  submucosal  layer. 
Electrode  four  is  closest  to  the  intact  tissue  and  electrode  one  is 
furthest  away. 

Finally  the  remaining  half  of  the  submucosal  layer  was  removed  and 
as  shown  by  Figure  5.14,  all  slow  wave  activity  was  abolished. 


■« 


. 


' 


97 


Figure  5.12  Histology  slides  showing  (a)  the  circular  muscle  before 
the  removal  of  the  submucosal  layer,  (b)  the  circular 
muscle  after  the  removal  of  the  submucosal  layer,  and 
(c)  the  submucosal  layer.  Magnification:  x62.5. 
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Figure  5.13(a)  Slow  waves  recorded  from  the  mucosal  side  of  the  intact 

muscle  wall  before  removal  of  the  submucosa. 
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Figure  5.13(b)  Slow  waves  recorded  from  the  mucosal  side  of  the  intact 

muscle  wall  after  removal  of  the  submucosa  from  half 
of  the  specimen. 
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Figure  5.14  A  recording  from  the  mucosal  side  of  the  intact  muscle  wall 
after  complete  removal  of  the  submucosal  layer. 
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No  electrical  activity  was  found  in  the  isolated  submucosal  layer. 
Also,  upon  removal  of  the  submucosa  it  was  impossible  to  produce  any 
cyclic  activity  similar  to  slow  waves  through  external  electrical 
stimulation  of  the  circular  muscle. 

5.3.6  Results  from  the  Dual  Chamber  Studies 

Twelve  dual  chamber  studies  were  conducted  as  described  in  Chapter 
4,  section  4.3.6.  Six  of  these  experiments  were  conducted  on  the  intact 
muscle  wall  and  six  were  conducted  on  isolated  longitudinal  and  isolated 
circular  muscle.  For  the  intact  muscle  segments  the  signals  recorded 
on  the  mucosal  side  were  always  3  to  5  times  larger  in  amplitude  than 
those  on  the  serosal  side.  A  typical  recording  is  shown  in  Figure 
5.15.  Channel  1  is  the  signal  recorded  from  the  serosal  side  of  the 
tissue,  channel  2  is  the  signal  from  the  mucosal  side,  and  channel  3 
is  the  difference  of  channels  1  and  2.  These  signals  were  recorded 
on  an  fm  magnetic  tape.  They  were  later  studied  at  fifteen  times  normal 
speed  using  a  storage  oscilloscope  to  measure  the  relative  phase  angles 
of  the  serosal  signal  versus  the  mucosal  signal.  With  the  time  scale 
expanded  fifteen  times,  the  relative  phase  angle  of  channel  1  versus 
channel  2  was  too  small  to  be  measured. 

It  has  been  shown  by  Bortoff  (1965)  that  an  increased  concentration 
of  30  m  M  KC1  causes  depolarization  of  small  intestine  smooth  muscle 
and  abolishes  slow  waves  in  that  organ.  In  a  set  of  four  experiments, 
an  increased  concentration  of  30  m  M  KC1  applied  to  the  chamber  in 
contact  with  the  mucosal  side  significantly  reduced  the  amplitude  of  the 
slow  wave  on  both  sides.  This  occurred  in  four  out  of  four  experiments 
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Figure  5.15  Typical  signals  from  an  intact  muscle  specimen  in  the  dual 
chambered  bath. 

No.  1  is  the  serosal  signal,  No.  2  is  the  mucosal  signal, 
and  No.  3  is  the  difference  of  Nos.  1  and  2. 
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and  a  representative  recording  is  shown  in  Figures  5.16  (a)  and  (b). 
Washing  with  normal  Krebs-Ringer  solution  reversed  the  effect,  and  slow 
waves  returned  as  shown  in  Figure  5.16(c).  As  shown  by  Figure  5.16(d) 
an  increased  concentration  of  30  m  M  KC1  on  the  longitudinal  side  had 
no  measurable  effect. 

In  six  experiments  the  circular  and  longitudinal  layers  were 
separated  and  studied  independently.  Normal  slow  waves  were  recorded 
from  both  sides  of  the  isolated  circular  muscle  layer.  The  amplitude 
of  the  signal  recorded  on  the  serosal  side  was  of  the  same  order  of 
magnitude  as  that  recorded  on  the  mucosal  side.  Again  the  average  phase 
angle  of  one  side  with  respect  to  the  other  was  too  small  to  be  measured. 
A  typical  recording  is  shown  in  Figure  5.17. 

Increasing  the  KC1  concentration  to  30  m  M  on  the  mucosal  side 
greatly  reduced  the  slow  wave  amplitudes  measured  on  both  sides  of  the 
circular  muscle  layer.  Again  the  effect  was  reversed  by  washing  with 
normal  Krebs-Ringer  solution.  Increasing  the  KC1  concentration  on  the 
serosal  side  of  the  circular  muscle  layer  produced  no  immediate  effect, 
but  abolished  slow  waves  after  a  twenty  to  thirty  minute  delay. 

The  isolated  longitudinal  muscle  layer  produced  no  slow  wave 
activity.  However,  another  type  of  signal  was  observed.  This  signal, 
shown  in  Figure  5.18,  consisted  of  bursts  of  activity  with  a  mean 
frequency  of  23.5  cycles/minute  and  an  average  of  19.25  cycl es/burst . 

The  bursts  occurred  1.82  times/minute  and  were  not  always  present. 

Their  occurrence  appeared  to  be  a  function  of  the  tension  with  which 
the  specimen  was  mounted,  since  it  could  be  abolished  by  reducing 
tissue  tension. 

Neither  this  activity  nor  the  normal  slow  wave  was  affected  by 
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Figure  5.16(a)  Dual  chamber  recording  with  normal  Krebs-Ringer  solution. 

Channels  1,  2,  and  3  are  as  in  Fig.  5.15. 
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KC I  concentration  changed  to  30  mM 


Figure  5.16(b)  KC1  concentration  changed  to  30  mM  on  the  mucosal  side 

of  the  bath. 

Channels  1,  2,  and  3  are  as  in  Fig.  5.15. 
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Figure  5.16(c)  Tissue  washed  with  norami  Krebs-Ringer  solution. 

Channels  1,  2,  and  3  are  as  in  Fig.  5.15. 
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KCI  concentration  changed  to  30  mM  on  the  longitudinal 
side  of  the  bath. 

Channels  1,  2,  and  3  are  as  in  Fig.  5.15. 
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Figure  5.17  Dual  chamber  study  of  isolated  circular  muscle. 
Channel  1  is  from  the  longitudinal  side, 

Channel  2  is  from  the  mucosal  side,  and 
Channel  3  is  the  difference. 
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Figure  5.18  Burst  activity  recorded  from  isolated  longitudinal  muscle. 
Channels  1,  2,  and  3  are  as  in  Fig.  5.17. 
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the  administration  of  10  and  10  M  tetrodotoxin  to  either  chamber 
of  the  bath. 

Another  measurement  of  interest  was  that  of  the  potential  of  the 
reference  electrode  in  the  chamber  adjacent  to  the  circular  muscle  with 
respect  to  the  reference  electrode  in  the  other  chamber.  Records  of 
two  such  measurements  are  shown  in  Figure  5.19.  Trace  5.19(a)  was 
recorded  in  the  presence  of  an  intact  muscle  specimen  while  trace 
5.19(b)  was  recorded  in  the  presence  of  an  isolated  circular  muscle 
specimen.  The  signals  associated  with  both  types  of  specimens  were 
cyclic  with  the  same  repetition  frequency  as  the  slow  wave.  However, 
the  signals  associated  with  isolated  circular  specimens  were  much 
smaller  in  amplitude  than  those  associated  with  intact  muscle  specimens. 

5.4  INTRACELLULAR  RESULTS 

5.4.1  Circular  Muscle  Measurements  on  the  Intact  Muscle 

As  described  in  Chapter  4,  section  4.4.3,  ten  i ntracel 1 ul ar  experi¬ 
ments  were  conducted  on  the  mucosal  side  of  intact  muscle  strips.  In 
eight  of  these  experiments  successful  cell  penetrations  were  achieved 
and  maintained  for  up  to  twenty  minutes.  Two  experiments  were  not 
successful  in  achieving  cell  penetrations  of  sufficient  duration  to  be 
included  in  the  results.  A  typical  intracellular  recording  is  shown 
in  Figure  5.20.  Both  channels  1  and  2  represent  the  same  i ntracel 1 ul ar 
potential.  However,  channel  1  is  calibrated  to  permit  measurement  of 
the  dc  i ntracel 1 ul ar  resting  potential  while  channel  2  has  the  dc 
offset  removed  to  permit  a  higher  amplifier  gain  and  a  more  exact  display 
of  the  slow  wave  itself. 
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Figure  5.19  Recordings  of  the  potential  of  the  circular  reference 

electrode  with  respect  to  the  longitudinal  reference  elec¬ 
trode  in  the  presence  of  (a)  an  intact  muscle  specimen,  and 
(b)  an  isolated  circular  muscle  specimen. 
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Figure  5.20  A  typical  i ntracel 1 ul ar  recording  from  a  circular  muscle 

cell.  Channel  1  is  calibrated  to  show  the  resting  potential 
and  Channel  2  is  calibrated  to  better  show  the  slow  wave 
amplitude. 
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In  addition  to  the  above  eight  experiments,  intracel 1 ul ar  recordings 
from  the  mucosal  side  of  the  intact  muscle  wall  were  obtained  from 
thirteen  other  experiments. 

As  shown  in  Table  5.4,  233  successful  penetrations  were  obtained 
in  the  circular  muscle  cells  of  38  specimens.  The  mean  values  of  the 
resting  potential  and  the  slow  wave  in  the  circular  muscle  layer  of 
proximal  and  transverse  canine  colon  were  calculated  to  be  -77.8  ± 

6.8  mV  and  28.9  ±  5.5  mV  respectively. 

5.4.2  Measurements  on  Isolated  Circular  Muscle 

Five  experiments  were  conducted  on  isolated  muscle  layers.  These 
layers  were  prepared  as  described  in  Chapter  4  and  were  mounted  in  the 
same  manner  as  the  intact  muscle  specimens.  Table  5.5  is  a  summary  of 
the  results  from  these  experiments.  The  mean  resting  potential  and 
the  slow  wave  amplitude  on  the  mucosal  side  of  the  isolated  circular 
muscle  layer  were  -73.9  ±  5.9  mV  and  30.7  ±  4.9  mV  respecti vely .  These 
values  are  not  statistically  different  from  those  obtained  for  intact 
muscle  strips  (Student's  t  Test  p  <  .05).  However,  on  the  longitudinal 
side  of  the  isolated  circular  muscle  layer  the  resting  potential  and 
slow  wave  amplitude  were  very  different.  From  twenty-eight  measurements 
the  mean  resting  potential  and  slow  wave  amplitude  were  computed  to  be 
-54.3  ±  5.3  mV  and  6.6  mV  ±  2.0  respectively.  The  explanation  for  this 
difference  may  be  due  to  cell  damage  caused  by  the  dissection  process. 

In  two  experiments  intracellular  measurements  were  obtained  from 
the  circular  muscle  layer  after  the  removal  of  the  submucosa.  These 
measurements  confirmed  the  absence  of  the  slow  wave  and  verified  the 
extracel 1 ul ar  results  of  section  5.3.5.  The  mean  resting  potential 
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Table  5.5  A  summary  of  intracel lular  results  from 

isolated  longitudinal  and  circular  layers. 
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of  the  circular  muscle  cells  after  the  removal  of  the  submucosal  layer 
was  found  to  be  -72.4  ±  19.9  mV. 

5.4.3  Longitudinal  Muscle  Measurements 

A  total  of  thirty  successful  penetrations  in  the  longitudinal 
muscle  cells  of  six  different  specimens  were  obtained.  In  none  of 
these  measurements  was  slow  wave  activity  observed.  A  typical  recording 
from  a  longitudinal  cell  penetration  is  shown  in  Figure  5.21.  The 
data  from  these  measurements  are  summarized  in  Table  5.4.  The  mean 
resting  potential  for  longitudinal  muscle  cells  was  -52.6  ±  7.0  milli- 
vol ts . 

The  results  from  the  five  experiments  on  isolated  muscle  layers 
also  confirmed  the  absence  of  slow  waves  in  isolated  longitudinal 
muscle.  Fifteen  measurements  from  three  different  specimens  of  isolated 
longitudinal  muscle  produced  a  mean  resting  potential  of  -54.3  ±  6.6 
millivolts.  This  is  not  statistically  different  from  the  value  computed 
for  longitudinal  cells  of  the  intact  muscle  wall. 

Figure  5.22  presents  another  interesting  result  obtained  from  two 
specimens  of  isolated  longitudinal  muscle.  This  recording  shows  a 
high  frequency  burst  of  activity  not  unlike  that  previously  observed 
in  the  extracellular  experiments  (section  5.3.6).  This  activity  was 
not  present  in  all  specimens.  The  mean  resting  potential  during  this 
activity  was  -52.5  ±  9.9  mV.  The  mean  frequency  of  the  activity  was 
16.81  ±  2.6  cycles/minute.  The  repetition  rate  of  the  bursts  was  not 
always  regular,  but  when  it  was,  the  mean  duration  of  each  burst  was  23.6 
seconds  and  the  mean  repetition  rate  was  0.71  cycles/minute. 

The  interesting  aspect  of  this  data  is  that  it  can  be  thought  of 


5(1 


117 


• - 30  sec - * 

iiiiiiiiiiiiiiiiiiiiiiiiimmiiiiiiiiiiimiiinii]iiiiii7i§fliiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!ii 
0  volts  — - j 

\ - 

-100  mV —  ’ 


4olmV 


Figure  5.21  A  typical  intracellular  recording  from  a  longitudinal  cell. 
Channels  1  and  2  are  as  in  Fig.  5.20. 
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Figure  5.22  Burst  activity  recorded  i ntracel 1 ul arl y  from  a  longitud¬ 
inal  cell  penetration. 
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as  the  mixing  (interaction)  of  two  i ntracel 1 ul ar  signals.  Frequently, 
extracel 1 ul arl y  recorded  slow  waves  produced  signals  similar  to  that 
shown  in  Figure  5.22.  However,  this  extracellular  phenomena  had  been 
assumed  to  be  the  result  of  the  beating  of  oscillator  cell  groups  of 
slightly  different  frequency  (Kingma  et.  al . ,  1979).  The  observations 
from  this  experiment  create  the  possibility  that  the  waxing  and  waning 
of  extracel 1 ul arl y  recorded  slow  waves  could  be  occurring  at  the 
cel  1 ul ar  1 evel . 

5.5  RESULTS  FROM  STUDIES  OF  THE  INTRACELLULAR-TO-EXTRACELLULAR 

TRANSDUCTION  PROCESS 

5.5.1  A  Comparison  of  Simultaneous  Intracellular  and  Extracellular 
Measurements 

As  was  described  in  Chapter  4  section  4.5.2,  three  experiments  were 
conducted  to  simultaneously  record  and  compare  extracellular  and  intra¬ 
cellular  potentials. 

A  typical  recording  from  this  set  of  experiments  is  shown  in 
Figure  5.23.  Channels  1  and  2  show  the  intracellular  signal  and  Channel 
3  the  extracel  1  ul  ar  signal. 

A  study  of  the  results  from  these  three  experiments  showed  clearly 
that  the  extracel 1 ul ar  waveform  is  a  function  of  the  i ntracel  1  ul  ar 
variations  in  cells  up  to  5  mm  away  from  the  extracel 1 ul ar  recording 
site.  The  phase  angle  between  the  two  sites  ranged  from  zero  to 
fifteen  degrees  but  the  variation  in  extracel 1 ul ar  waveform  ^hape 
made  it  difficult  to  obtain  a  mean  value  for  the  phase  angle.  In 
addition,  the  limited  number  and  short  duration  of  these  measurements 
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Figure  5.23  Simultaneous  recording  of  i ntracel 1 ul ar  and  extracellular 
si gnal s . 
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preclude  the  computation  of  a  statistically  significant  value  for 
the  phase  angle. 

5.5.2  Injury  Potential  Measurement 

Eight  experiments  were  described  in  Chapter  4  section  4.5.3  which 
were  conducted  to  determine  if  an  "injury  potential"  is  responsible 
for  the  extracel 1 ular  pressure  electrode  signal.  Figure  5.24  is  a 
typical  recording  from  one  of  these  eight  experiments.  The  top  two 
tracings  show  the  i ntracel 1 ul ar  potential  recorded  from  a  circular 
muscle  cell  under  the  simulated  pressure  electrode  when  no  force  was 
applied.  The  bottom  two  tracings  show  the  intracel 1 ul ar  potential 
under  the  same  pressure  electrode  at  the  same  location  but  with  a 
one  gram  force  applied. 

The  results  of  the  eight  experiments  are  given  in  Table  5.6. 

V  and  V  .  indicate  values  of  the  maximum  and  minimum  i ntracel 1 ul ar 
max  min 

potentials  which  occur  at  the  slow  wave  plateau  and  trough  respecti vel y . 

To  investigate  the  relationship  between  the  intracellular  poten¬ 
tials  and  the  applied  force  of  the  extracel 1 ul ar  electrode,  the  mean 
potentials  associated  with  each  application  force  are  given  in  Table 
5.7.  One  of  the  eight  experiments  showed  only  a  small  change  in  slow 
wave  amplitude  as  a  result  of  the  pressure  electrode.  The  remaining 
seven  experiments  showed  a  significant  reduction  in  slow  wave 
amplitude  from  a  mean  value  of  28.9  ±  4.0  mV  to  a  mean  value  of  14.43 
±  4.3  mV  with  the  application  of  forces  ranging  from  0.5  to  1  gram. 

As  the  force  of  the  pressure  electrode  increased  the  slow  wave  trough 
(or  resting  potential)  increased  from  a  mean  of  -76.4  ±  6.8  mV  to  a 
mean  of  -59.7  ±  7.0  mV.  However,  the  level  of  the  slow  wave  plateau 
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Figure  5.24  Intracellular  potentials  in  the  presence  of  an  extracel 1 ul ar 
pressure  electrode. 

Channel  1  is  calibrated  for  dc  potentials  and 
Channel  2  has  the  dc  component  removed. 


Table  5.6  Results  from  experiments  to  determine  the  affect 
of  extracellular  electrode  force  on  intracellular 
potential . 


EXP. 

NO. 

FORCE(GMS) 

N 

^max  (mV) 

^min  (mV) 

SLOW  WAVE 
AMPLITUDE  (mV) 

60 

0 

21 

-52 . 4±1 0. 9 

-76 . 9±1 0.4 

26.3i3.8 

1 

1 

-70.0 

-70.0 

0.0 

61 

0 

18 

-43 . 0±6 . 3 

-73 . 8±7 . 6 

27 . 7i7 . 3 

1 

4 

-45 . 0±4 . 4 

-70. 0±7 . 1 

25. 3i8. 1 

2 

5 

-54 . 0±1 5 . 2 

-79 . 8±8. 7 

25 .8i7 . 3 

3 

2 

-43. 5± . 71 

-65 . 0±7 . 1 

21 . 5i6. 4 

63 

0 

9 

-51 . 1 ±4. 6 

-74 . 0±3 . 5 

23 . 3i3 . 7 

1 

3 

-47 . 0±0 

-63 . 3±2 . 9 

1 4 . 7i2 . 1 

2 

5 

-44.4±4.0 

-54 . 0±6 . 5 

13. 3i2 . 2 

3 

3 

-47 . 5±4. 1 

-54. Oil .7 

9 . 3il  .2 

4 

3 

-47 . 3±2 . 5 

-58 . 5i4. 9 

10. OiO. 0 

5 

2 

-40. 0±0 . 0 

-50. OiO. 0 

10. OiO. 0 

68 

0 

5 

-45 . 0±1 .4 

-81 . 0i3 . 8 

33 . 2i3. 3 

1 

2 

-41 .0±1 .4 

-51 . 0i5 . 7 

1 0. 5i4. 9 

0 

3 

-44 . 7±4 . 0 

-75. OiO. 0 

30. 3i4 . 0 

1 

2 

-42 . 0±2 . 8 

-48 . 0i5 . 7 

6 . 0i2 . 8 

69 

0 

3 

-40. 3±6 . 4 

-82 . 3i2 . 5 

40. 0i6 .1 

SPEC . 1 

0.5 

1 

-60.0 

-75.0 

15.0 

0 

3 

-43 . 3±1 4 . 4 

-71 . 3i9. 3 

28 . 0i5 . 3 

SPEC. 2 

0 

1 

-39.0 

-76.0 

34.0 

0.5 

2 
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-57. OiO 

1 3 . 5i0. 71 

1 

2 

-45 . 0±0. 0 

-50. OiO 

5 . 5i0. 7 

0 

1 
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-96.0 

24.0 

73 

0 

13 

-53. 1 ±8. 3 

-79. 6i5 . 6 

26 . 7i5 . 6 

1 

2 

-47 . 5±3 . 5 
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19. Oil  .4 

2 

3 
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-64 . 7i0. 6 

1 4 . 7i0. 6 

4 

2 
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-60. 0i7 . 1 
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0 

1 
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74 

0 
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-43 . 3±1 0. 0 
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36 . 3±1 . 6 

0.5 

3 
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-73 .3i2 .4 

29.8i4.2 

75 

0 

1 

-35.0 

-80.0 

45.0 

SPEC . 1 

0.5 

3 

-55 . 7±4 . 0 

-74. Oil .7 

1 9.0±4. 0 

1 

5 

-42 . 8±7 . 0 

-59.0i5. 5 

15.0+3.1 

0 

4 

-36. 3±4. 8 

-70. 5i6. 4 

35 . 5i3. 7 

SPEC. 2 

0 

1 

-35.0 

-70.0 

35.0 

0.5 

1 

-35.0 

-70.0 

34.0 
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Table  5.6  Results  from  experiments  to  determine  the  affect 
of  extracellular  electrode  force  on  i ntracel 1 ul ar 
potential  (cont.). 


EXP. 

NO. 

FORCE (GMS) 

N 

^max  (mV) 

^min  (mV) 

SLOW  WAVE 
AMPLITUDE  (mv) 

1 

1 

-47.0 

-65.0 

18.0 

2 

1 

-45.0 

-65.0 

20.0 

3 

1 

-43.0 

-68.0 

25.0 

SPEC. 3 

0 

3 

-44. 3±1 .2 

-70. 7±1 .2 

26 . 7±2 . 9 

0.5 

2 

-41 . 5±2 . 1 

-73 . 5±4 . 9 

31 .0±1 .4 

1 

2 

-45 . 0±0. 0 

-62.0±0. 0 

1 7 . 0±0 

2 

1 

-45.0 

-70.0 

25.0 

2.5 

1 

-50.0 

-65.0 

15.0 

3 

1 

-45.0 

-60.0 

15.0 

4 

2 

-49. 9±1 .4 

-63. 5±4 . 9 

14. 5±3. 5 

5 

2 

-47 . 5±3 . 5 

-57 . 5±3 . 5 

10.0±0 

0 

4 

-45 . 3±2 . 1 

-70. 8±1 .5 

26 . 5±2 . 6 
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Table  5.7  Summary  of  the  results  of  Table  5.6  arranged  in 
order  of  increasing  electrode  force. 


FORCE(GM) 

NO.  OF 
PENETRATIONS 

Vmax  (mV) 

^min  (mV) 

S.W.  AMPLITUDE  (mV) 

0 

83 

-47 . 9±7 . 9 

-76.4±6.8 

28 . 9±4 . 0 

.5 

12 

-44.8±5.2 

-63 . 6±1 3 . 1 

1  9.6±8.4 

1 

20 

-45 . 7±3 . 8 

-58. 6±3 . 9 

12.7  ±2 . 5 

2 

10 

-46.2±2.8 

-59 . 9±4 . 6 

1 5 . 6±1 .6 

3 

5 

-46.1 ±3.2 

-58. 0±1 .3 

13. 5± . 9 

4 

7 

-47.8+2.6 

-60. 4±5 . 6 

1 2 . 0±2 . 6 

5 

4 

-43.8±2.5 

-53.8±2. 5 

lO.OiO.O 
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increased  only  a  very  small  amount  from  -47.9  ±  7.9  mV  to  -45.8  ± 

3.8  mV  for  the  full  range  of  pressure  electrode  forces.  In  some 
experiments  the  plateau  level  remained  unchanged.  Figure  5.25  is 
a  plot  of  slow  wave  amplitudes  (normalized  with  respect  to  basal 
values)  versus  pressure  electrode  application  force  in  grams.  These 
graphs  show  that  the  relationship  between  applied  force  and  slow  wave 
amplitude  is  nonlinear.  However,  for  a  force  of  1  gm  or  greater  the 
mean  slow  wave  amplitude  is  reduced  to  60  percent  or  less  of  its  basal 


value. 


SLOW  WAVE  AMPLITUDE  NORMALIZED  W.R.T.  BASAL  AMPLITUDE 


127 


1.25  r- 


1.00 


0  75  - 


0.50 


0  25- 


2  3 

ELECTRODE  FORCE  (GRAMS) 


KEY 

EXPERIMENT  60 
EXPERIMENT  61 
EXPERIMENT  63 
EXPERIMENT  68 
EXPERIMENT  69  a 
EXPERIMENT  69  b 
EXPERIMENT  73 
EXPERIMENT  74 
EXPERIMENT  75  a 
EXPERIMENT  75  b 
EXPERIMENT  75  c 
—  MEAN 


Figure  5.25  Normalized  slow  wave  amplitude  versus  extracellular  electrode 
application  force. 


CHAPTER  6 


CONCLUSIONS 


6.1  INTRODUCTION 

The  results  of  Chapter  5  give  a  better  understanding  of  the  origin, 
i n tracel 1 ul ar  characteristics ,  and  coupling  of  canine  colon  electrical 
activity.  The  results  also  give  some  insight  into  the  role  played 
by  each  of  the  muscle  layers  in  the  generation  of  this  activity  and 
into  the  mechanism  for  intracell ul ar-to-extracel 1 ul  ar  transduction 
of  slow  wave  signals.  In  this  chapter  each  of  these  areas  is  examined 
and  conclusions  from  the  results  of  the  previous  chapter  are  given. 
Intracel 1 ul ar  results,  single-sided  extracel 1 ul ar  recordings,  and 
double-sided  extracellular  recordings  are  compared  to  present  infor¬ 
mation  in  support  of  each  of  the  conclusions. 

Following  this,  electrical  models  are  presented  of  slow  wave 
generation  at  the  cellular  level  and  of  the  i ntracel 1 ul ar-to-extra- 
cellular  transduction  mechanism.  A  model  is  also  presented  to  show 
the  roles  played  by  each  of  the  muscle  layers  in  the  generation  and 
coupling  of  slow  wave  activity. 

The  chapter  concludes  with  an  indication  of  the  direction  of 
further  work. 

6.2  SLOW  WAVE  ORIGIN 

The  following  conclusions  may  be  drawn  concerning  the  origin 
of  slow  wave  activity: 
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1.  The  slow  wave  originates  in  a  layer  of  cells  at  the  boundary  of 
the  circular  muscle  layer  and  the  submucosa. 

2.  This  layer  is  essential  for  the  generation  of  slow  wave  activity 
in  the  circular  muscle  layer. 

3.  Slow  waves  are  not  generated  in  the  cells  of  the  longitudinal 
muscle  layer. 

These  conclusions  are  supported  by  the  results  of  the  single-sided  and 
double-sided  studies  of  the  intact  muscle  wall  and  isolated  muscle 
layers  as  well  as  by  i ntracel 1 ul ar  measurements. 

The  isolated  layer  studies  of  Chapter  5,  section  5.3.4  provide 
information  relating  to  the  origin  of  the  slow  wave  activity.  In  five 
of  seven  experiments,  slow  waves  were  recorded  extracel 1 ul arl y  from 
both  sides  of  the  isolated  circular  muscle  layer.  However,  in  none  of 
four  experiments  were  slow  waves  recorded  extracel 1 ul arly  on  either 
side  of  the  isolated  longitudinal  muscle  layer.  This  led  to  the 
suspicion  that  the  slow  wave  origin  was  either  in  the  circular  muscle 
layer  or  its  associated  submucosal  layer. 

As  was  described  in  Chapter  5,  section  5.3.5,  the  submucosal  layer 
could  easily  be  removed  from  the  circular  muscle  layer.  In  five 
experiments  using  twelve  different  specimens  the  normal  slow  wave 
in  the  circular  muscle  layer  was  abolished  by  the  complete  removal 
of  the  submucosal  layer.  However,  when  the  submucosal  layer  was 
removed  from  only  part  of  the  specimen,  the  slow  wave  was  abolished 
only  in  the  region  from  which  the  submucosa  had  been  removed.  As  seen 
in  Figure  5.13,  the  amplitude  of  the  slow  wave  in  the  region  devoid 
of  the  submucosa  was  a  function  of  the  distance  from  the  intact  portion 
of  the  tissue.  The  circular  muscle  layer  thus  requires  the  presence 
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of  the  submucosal  layer  to  produce  slow  wave  activity. 

The  dual  chamber  studies  provided  additional  information  about 
the  slow  wave  origin.  The  pertinent  results  are  as  follows: 

1.  Slow  waves  recorded  simultaneously  on  both  sides  of  the  intact 
muscle  wall  had,  in  every  experiment,  larger  amplitudes  on  the 
mucosal  side  of  the  specimen. 

2.  An  increased  concentration  of  30  m  M  KC1  applied  to  the  mucosal 
side  of  an  intact  muscle  specimen  in  all  cases  significantly 
reduced  slow  wave  activity  as  recorded  on  both  sides  of  the  specimen. 
However,  the  same  increased  concentration  applied  to  the  serosal 
side  produced  no  effect.  Thus  the  mucosal  side  of  the  circular 
muscle  layer  is  of  particular  significance  with  regard  to  the 
generation  of  the  slow  wave. 

3.  Slow  waves  were  recorded  on  both  sides  of  isolated  circular  muscle 
and  were  always  of  approximately  the  same  amplitude. 

4.  Slow  waves  were  never  recorded  on  either  side  of  isolated  longi¬ 
tudinal  muscle. 

The  most  convincing  evidence  relating  to  the  origin  of  the  slow 
wave  are  the  results  of  the  intracel 1 ul ar  experiments  present  in  Chapter 

5.  section  5.4.  Slow  waves  were  recorded  from  233  penetrations  of  cells 
located  throughout  the  circular  muscle  layers  of  37  intact  specimens. 

Slow  waves  were  also  recorded  throughout  the  isolated  circular  muscle 
layer  when  its  submucosa  was  intact.  However,  slow  waves  were  not 
recorded  from  the  circular  muscle  layer  after  the  submucosa  had  been 
removed.  In  28  penetrations  of  longitudinal  muscle  cells  in  the  intact 
wall,  slow  wave  activity  was  never  recorded.  This  was  also  observed 
for  isolated  longitudinal  muscle. 
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To  recapitulate,  slow  waves  do  not  exist  inside  longitudinal 
muscle  cells,  but  they  originate  in  the  circular  muscle  layer  at  the 
boundary  with  the  submucosa.  Cells  within  the  submucosa  or  special 
cells  near  the  inner  surface  of  the  circular  muscle  layer  interact 
with  the  cells  of  the  circular  muscle  layer  to  generate  slow  waves. 

Two  of  the  more  likely  explanations  of  this  interaction  are: 

(1)  That  pacemaker  cells  within  the  submucosal  layer  or  in 
the  immediately  adjacent  circular  muscle  layer  provide  an  elec¬ 
trical  stimulus  to  the  circular  muscle  layer, 
or  (2)  That  cells  in  the  submucosal  layer  provide  a  chemical  input 
to  the  circular  layer  changing  the  resting  potential  to  cause 
slow  wave  generation  within  the  inner  circular  layer. 

Slow  waves  recorded  extracel 1 ul arly  on  the  longitudinal  side  of 
the  intact  muscle  wall  are  the  result  of  passive  current  flow  from  the 
circular  muscle  layer  through  the  longitudinal  muscle.  This  is  confirmed 
by  the  lower  amplitude  slow  wave  recorded  on  the  longitudinal  side  of 
the  specimen  in  the  dual  chamber  experiments.  The  potentials  generated 
by  the  circular  muscle  layer  are  attenuated  by  the  resistance  of  the 
longitudinal  layer.  This  is  also  supported  by  results  from  the  measure¬ 
ment  of  the  potential  of  the  longitudinal  reference  electrode  with 
respect  to  the  circular  reference  electrode  as  described  in  Chapter 
5,  section  5.3.6.  In  the  presence  of  an  intact  muscle  specimen  this 
potential  varied  at  the  same  frequency  as  the  slow  wave.  However, 
these  variations  were  an  order  of  magnitude  smaller  with  isolated 
circular  muscle  in  the  tissue  bath.  Slow  wave  currents  flowing  from 
the  circular  layer  through  the  longitudinal  layer  cause  the  potential 
of  the  circular  side  to  vary  with  respect  to  the  serosal  side.  However, 


. 
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in  the  esse  of  isolated  circular  muscle  the  slow  wave  currents  from  the 
circular  muscle  cells  flow  directly  to  each  reference  electrode  passing 
through  only  the  low  resistance  of  the  tissue  bath  on  both  sides  of 
the  tissue. 

6.3  INTRACELLULAR  SLOW  WAVE  CHARACTERISTICS 

Prior  to  this  work,  intracellular  recordings  from  the  canine 
colon  have  not  been  reported  in  the  literature.  Consequently,  this 
work  has  produced  the  first  comprehensive  intracellular  study  of  canine 
slow  wave  characteristics. 

As  reported  in  Chapter  5,  section  5.4,  233  successful  penetrations 
were  obtained  from  the  circular  muscle  layers  of  37  different  specimens 
of  intact  muscle  wall  from  the  proximal  and  transverse  colon.  There 
was  no  variation  in  cellular  resting  potential  or  slow  wave  amplitude 
throughout  the  circular  muscle  layer  or  over  the  regions  of  the  colon 
studied.  The  mean  resting  potential  and  slow  wave  amplitude  were 
-77.8  ±  6.8  mV  and  28.9  ±  5.5  mV  respecti vely .  From  the  results  of 
five  experiments  (described  in  Chapter  5,  section  5.4.2)  conducted  on 
the  isolated  circular  muscle  layer,  the  mean  resting  potential  and 
slow  wave  amplitude  on  the  mucosal  side  of  this  layer  were  found  to 
be  -73.9  ±  5.9  mV  and  30.7  ±  4.9  mV  respectively .  This  is  not  signif- 
igantly  different  (Student's  t  test,  p  <  .05)  from  the  corresponding 
results  from  the  intact  muscle  wall.  Therefore,  the  separation  of 
the  layers  does  not  affect  slow  wave  generation.  However,  the  intra¬ 
cellular  results  confirmed  that  the  removal  of  the  submucosal  layer 
abolished  the  slow  wave  in  the  circular  muscle  layer  (Chapter  5, 
section  5.4.2). 
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The  intracellular  results  from  the  circular  muscle  layer  also 
confirmed  that  the  extracel 1 ul arl y  recorded  slow  wave  signal  is  an 
accurate  representation  of  the  intracellular  potential  variations 
of  the  cells  surrounding  the  electrode. 

Intracellular  measurements  were  also  obtained  from  the  longitudinal 
muscle  cells  in  both  the  intact  muscle  wall  and  the  isolated  longitu¬ 
dinal  layer.  Because  the  longitudinal  layer  is  much  thinner  than  the 
circular  layer  and  because  the  serosa  was  an  impediment  to  electrode 
penetration,  the  number  of  successful  penetrations  in  the  longitudinal 
layer  were  less  than  that  in  the  circular  layer.  Nevertheless,  28 
successful  penetrations  were  obtained  in  six  different  specimens  of 
intact  muscle  wall  and  15  successful  penetrations  were  obtained  in 
three  different  specimens  of  isolated  longitudinal  muscle.  Slow  waves 
were  not  observed  in  any  of  these  recordings.  The  mean  resting  poten¬ 
tials  for  intact  and  isolated  longitudinal  muscle  were  -52.6  ±  7.0  mV 
and  -54.3  ±  6.6  mV.  Again,  the  process  of  muscle  separation  has  not 
affected  the  intracellular  potentials. 

Another  interesting  result  obtained  from  the  intracellular  experi¬ 
ments  was  the  recording  of  a  higher  frequency  burst  activity  observed 
i ntracel 1 ul arl y  in  longitudinal  muscle  cells  (Chapter  5,  section  5.4.3 
and  Figure  5.21).  This  activity  was  also  observed  extracel 1 ul arl y 
as  reported  in  the  results  of  the  dual  chamber  studies  (Chapter  5, 
section  5.3.6)  and  has  been  reported  by  El-Sharkawy  (1979).  The 
intracellular  results  of  section  5.4.3  indicate  that  this  burst  activity 
is  an  i ntracel 1 ul ar  phenomena  and  not  the  result  of  the  mixing  of 
signals  from  different  cell  groups. 
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6.4  ELECTRICAL  COUPLING  WITHIN  THE  COLON 

The  data  from  23  single-sided  extracellular  experiments  was 
processed  on  the  HP  21MX/E  computer  (as  described  in  Chapter  5,  section 
5.3.1)  to  determine  the  relative  phase  lag  between  recording  sites. 

The  results  from  this  processing  confirmed  that  electrical  coupling 
between  oscillating  regions  in  the  colon  is  very  poor.  There  was 
great  variability  in  the  relative  phase  angles  of  slow  waves  recorded 
from  sites  separated  by  as  little  as  2.5  mm.  However,  from  the  large 
amount  of  data  processed  it  was  determined  that  the  mean  phase  angle 
between  recording  sites  separated  by  2.5  mm  was  29.0  ±  32.2  degrees 
in  the  circular  direction  and  117.3  ±  51.0  in  the  longitudinal  direc¬ 
tion.  A  conclusion  to  be  deduced  from  this  data  is  that  coupling  is 
much  better  in  the  circular  direction  than  in  the  longitudinal 
direction. 

The  results  from  the  cutting  experiments  (described  in  Chapter 
5,  section  5.3.3)  also  support  the  hypothesis  that  coupling  between 
regions  is  better  in  the  circular  direction  than  in  the  longitudinal 
direction.  In  addition  to  this,  these  results  show  that  no  gradient 
of  natural  frequencies  exists  in  the  longitudinal  direction  but  that 
such  a  gradient  exists  in  the  circular  direction  with  respect  to  the 
mesentery.  This  indicates  (as  stated  in  section  5.3.3)  that  the 
mesentery  plays  a  role  in  determining  slow  wave  frequency. 

6.5  PRESSURE  ELECTRODE  SIGNAL  TRANSDUCTION 

The  currently  accepted  theory  for  the  generation  of  the  pressure 
electrode  slow  wave  signal  is  that  it  results  from  an  injury  potential 
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caused  by  cell  distortion.  This  theory  was  reviewed  in  Chapter  3  and 
extended  to  explain  the  situation  when  many  cells  are  electrically 
coupled.  However,  descrepanci es  exist  between  the  theory  described 
in  Chapter  3  and  the  results  of  the  transduction  experiments  described 
in  Chapter  5,  section  5.5.2. 

The  extension  of  Gillespie's  model  (Gillespie,  1952)  as  described 
in  Chapter  3  proposes  that  cell  distortion  causes  membrane  depolariza¬ 
tion  resulting  from  a  drop  in  membrane  resistance.  This  drop  in 
membrane  resistance  was  represented  in  the  model  of  Figure  3.6  by  the 
variable  resistor  .  However,  this  drop  in  resistance  would  result 
in  the  reduction  of  both  the  dc  and  the  ac  components  of  the  intra¬ 
cellular  potential.  Both  ac  and  dc  components  of  the  intracellular 
potential  of  the  cells  of  the  rabbit  large  intestine  were  reduced 
by  the  same  factor  when  these  cells  were  injured  by  large  microelec¬ 
trodes  (Gillespie,  1962).  The  results  of  Chapter  5  include  only  the 
data  from  what  were  considered  to  be  successful  cellular  penetrations. 
However,  on  many  occasions  in  all  i ntracel 1 ul ar  experiments,  signals 
were  obtained  from  what  were  recognized  as  "unsuccessful  penetrations". 
These  unsuccessful  penetrations  were  identified  by  the  fact  that  the 
resting  potential  and  slow  wave  amplitude  were  not  maintained  constant 
during  the  penetration.  Immediately  after  the  penetration  the  resting 
potential  increased  slowly  towards  zero  potential.  In  these  cases 
the  slow  wave  amplitude  was  reduced  in  proportion  to  the  reduction 
in  the  magnitude  of  the  resting  potential.  The  value  of  the  plateau 
potential  of  the  slow  wave  was  also  changed  by  the  same  factor.  These 
"unsuccessful"  penetrations  produced  results  consistant  with  Gillespie's 
injury  potential  theory.  However,  the  results  from  the  study  of  the 
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pressure  electrode  mechanism  (Chapter  5,  section  5.5.2)  do  not  support 
the  injury  potential  hypothesis. 

In  seven  of  eight  experiments  the  application  of  a  pressure  elec¬ 
trode  with  a  force  of  0.5  to  1  gm  caused  a  significant  reduction  in 
slow  wave  amplitude  from  28.9  ±  4.0  mV  to  14.4  ±  4.3  mV.  However, 
this  reduction  occurred  without  a  change  in  the  slow  wave  plateau 
potential.  The  mean  plateau  potentials  before  and  after  the  pressure 
electrode  application  are  -47.9  ±  7.9  mV  and  -45.8  ±  3.8  mV  respec- 
ti vely. 

The  reduction  in  slow  wave  amplitude  results  almost  entirely 
from  an  increase  in  the  cellular  resting  potential.  The  "injury 
potential"  theory  requires  a  proportional  reduction  of  both  resting 
and  slow  wave  plateau  potentials. 

This  leads  to  the  following  conclusions:  First,  if  an  "injury 
potential"  is  created  by  the  pressure  electrode,  it  is  not  the  result 
of  a  total  reduction  in  membrane  resistance.  Second,  the  results  imply 
the  existence  of  separate  mechanisms  for  the  maintenance  of  the  cellular 
resting  potential  and  for  the  initiation  of  the  slow  wave  plateau 
phase.  Third,  the  distortion  caused  a  reduction  in  the  cellular  resting 
potential  but  did  not  interfere  with  the  ionic  mechanism  responsible 
for  the  initiation  of  the  slow  wave  plateau.  The  latter  two  of  these 
conclusions  require  further  clarification. 

6.6  A  PROPOSED  THEORY  FOR  SLOW  WAVE  GENERATION  AND  PRESSURE  ELECTRODE 

SIGNAL  GENERATION 


Since  it  was  not  one  of  the  objectives  of  this  thesis  to  examine 
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the  ionic  mechanisms  responsible  for  the  slow  wave,  the  above  conclusions 
are  implied  by  the  results  but  are  not  proven.  The  two  generally 
accepted  theories  with  respect  to  slow  wave  generation  involve  either 
modulation  of  the  Na-K  pump  or  modulation  of  the  membrane  impedance 
associated  with  such  a  pump.  In  either  case  if  the  slow  wave  were 
the  result  of  a  relaxation  oscillator  involving  a  membrane  voltage 
threshold,  it  would  be  expected  that  a  change  in  resting  potential 
would  result  in  a  change  in  slow  wave  frequency.  This  did  not  happen 
when  the  membrane  resting  potential  was  altered  by  the  application  of 
the  simulated  pressure  electrode.  This  is  an  indication  that  the  start 
of  the  slow  wave  was  not  the  result  of  a  threshold  phenomena  involving 
either  of  the  ions  associated  with  the  Na-K  pump.  The  time  of  the 
slow  wave  upswing  and  the  slope  of  the  intracellular  potential  prior 
to  the  upswing  were  unchanged  with  the  change  in  resting  potential. 

The  indication  is  that  some  other  ion  (which  has  little  effect  on  the 
membrane  potential)  is  responsible  for  the  slow  wave.  The  concentra¬ 
tion  of  such  an  ion  would  move  towards  a  threshold  value  while  the 
cellular  potential  is  at  the  resting  value.  When  the  threshold  concen¬ 
tration  is  attained,  this  ion  could  initiate  the  slow  wave  plateau 
phase  by  either  a  direct  effect  on  the  membrane  impedance  or  by 
inhibition  of  the  Na-K  pump.  During  the  plateau  phase  the  concentra¬ 
tion  of  this  ion  would  return  to  a  subthreshold  level.  El-Sharkaway 
and  Daniel  (1975  b)  and  Connor  (1979)  postulated  that  the  Ca  ion 
might  fulfill  such  a  role  in  the  stomach  and  small  intestine. 

The  plateau  level  of  the  slow  wave  is  not  significantly  changed 
by  the  pressure  electrode  application.  This  indicates  that  the  ions 
and  membrane  conductances  associated  with  the  plateau  are  not  altered 
by  cell  distortion.  This  work  provides  no  evidence  regarding  the  identity 
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of  such  ions  or  why  such  conductances  remain  unchanged. 

However  the  results  provide  support  for  the  hypothesis  that  slow 
waves  are  the  result  of  a  modulation  of  the  membrane  impedance.  The 
results  are  consistent  with  the  idea  that  the  membrane  impedance  consists 
of  two  parallel  components,  a  constant  impedance  which  is  primarily 
responsible  for  the  resting  potential  and  another  variable  impedance 
which  causes  the  membrane  potential  variations  called  slow  waves. 
Distortion  due  to  the  pressure  electrode  reduces  the  constant  component 
of  the  membrane  impedance  but  does  not  affect  the  variable  component 
and  consequently  does  not  affect  the  level  of  the  slow  wave  plateau. 

The  generation  of  the  cellular  resting  potential  and  the  slow 
wave  potential  variations  can  be  modelled  as  shown  in  Figure  6.1. 

This  model  is  not  designed  to  include  action  potential  generation. 

Rm  represents  the  normal  membrane  resistance  in  the  absence  of  an  action 
potential.  I  is  a  constant  current  source  representing  the  total 
current  produced  by  the  electrochemical  gradients  and  the  Na-K  pump. 

A  time  varying  current  source  models  the  slow  wave  potential 
varation. 

The  pressure  electrode  causes  cell  distortion  which  results  in  an 
increase  in  the  resting  potential.  The  i ntracel 1 ul ar  potential  rises 
to  a  resting  value  correspondi ng  to  the  level  of  the  slow  wave  plateau. 
This  is  illustrated  in  Figure  5.24.  Similarly,  damage  due  to  the 
dissection  process  produces  the  same  results.  This  was  verified  by 
the  data  from  the  intracellular  measurements  on  isolated  circular 
muscle  presented  in  section  5.4.2.  These  results  show  that  on  the 
longitudinal  side  of  the  isolated  circular  muscle  layer,  the  plateau 
potential  of  the  slow  wave  is  -47.7  mV,  corresponding  to  that  found  in 
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Intracellular  Fluid 


Figure  6.1  Model  of  the  smooth  muscle  membrane. 
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the  circular  layer  of  the  intact  muscle  wall.  However,  the  slow  wave 
amplitude  was  reduced  to  a  mean  of  6.6  mV  indicating  a  reduction  in  i  . 

The  generation  of  the  extracellular  pressure  electrode  signal  may 
be  explained  using  Figure  6.2.  The  current  source  i ^  represents  the 
equivalent  slow  wave  generator  of  all  cells  distorted  by  the  pressure 
electrode.  Rs  is  the  shunt  resistance  between  the  active  surface  of 
the  pressure  electrode  and  the  reference  electrode  in  the  tissue  bath. 

R.  and  R  are  the  equivalent  internal  and  coupling  resistances  of  the 

*1-  C 

cells  from  which  the  recording  is  obtained.  Assuming  an  infinite 
input  impedance  for  the  amplifier,  the  recorded  voltage  v(t)  as  deter¬ 
mined  from  the  circuit  of  Figure  6.2(b)  will  be 


6.1 


If  in  the  limit  V  (t)  is  reduced  to  zero  as  indicated  by  the  abolition 

s  w 

of  slow  waves  under  the  pressure  electrode,  then 


6.2 


/  ,  \  _  aw _ a 

vltj  2R  +  2R  .  +  R 


m  ^  c 


However  for  the  normal  cell,  the  slow  wave  component  vsw(t)  the 
in tracel 1 ul ar  potential  is  given  by 


v  (t)  =  -i  (t)  R 
sw  sw  m 


6.3 


Substituting  this  into  equation  6.1  gives 


2R  +  2R.  +  R 


R 

s 


•  v  (t)  =  K  •  v  (t) 
sw  sw 


6.4 


m  ^  c 


where  K  is  a  constant. 
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Figure  6.2  Extracellular  pressure  electrode  slow  wave  generation. 
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Therefore  the  extracellular  slow  wave  signal  is  proportional  to 

the  intracellular  potential  variation  v  (t).  If  the  amplifier  input 

impedance  is  finite,  then  in  the  above  expression  must  be  replaced 

by  the  value  of  the  parallel  combination  of  R  and  R.  . 

s  1  n 

6.7  AN  EXAMINATION  OF  THE  ROLES  PLAYED  BY  THE  INDIVIDUAL  MUSCLE  LAYERS 

IN  THE  GENERATION  OF  SLOW  WAVE  ACTIVITY 

6.7.1  The  Circular  Layer 

Evidence  presented  in  this  work  indicates  that  a  special  layer  of 
cells  at  the  boundary  of  the  submucosa  and  the  circular  muscle  layer 
is  essential  for  the  existence  of  slow  waves  in  the  circular  muscle 
layer.  In  the  absence  of  the  special  layer  the  circular  muscle  cells 
have  i ntracel 1 ul ar  potentials  which  remain  constant  at  approximately 
-72  mV.  Therefore  the  biological  oscillators  of  the  circular  muscle 
layer  are  either  a  monostable  type  requiring  an  input  trigger  signal 
from  the  special  layer  or  they  require  a  chemical  or  hormonal  input 
to  change  oscillator  parameters  to  cause  oscillations.  If  the  oscil¬ 
lators  are  of  the  monostable  type,  the  trigger  signal  from  the  special 
layer  must  be  generated  by  astable  oscillators  within  that  layer  or 
within  the  circular  layer  of  the  boundary  with  the  submucosa.  A  model 
showing  a  possible  arrangement  of  these  oscillators  is  shown  in  Figure 
6.3.  The  physiological  units  represented  by  the  astable  oscillators 
could  be  individual  cells  or  groups  of  tightly  coupled  cells  oscillating 
in  unison. 

The  results  from  the  cutting  experiments  (Chapter  5,  section 
5.3.3)  show  that  the  oscillators  of  the  colon  have  a  frequency  gradient 


. 
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with  respect  to  the  mesentery.  The  frequency  of  slow  wave  oscillations 
is  reduced  significantly  when  the  tissue  is  separated  from  the  mesentery. 

The  model  of  Figure  6.3  shows  coupling  to  exist  between  the 
oscillators  of  the  circular  muscle  layer.  This  was  shown  to  be  correct 
because  partial  removal  of  the  special  layer  did  not  result  in  total 
abolition  of  the  slow  wave  in  the  regions  devoid  of  the  special  layer. 

The  amplitudes  of  slow  waves  recorded  in  these  regions  diminished  with 
distance  from  the  portion  of  the  tissue  with  the  special  layer  intact. 

6.7.2  The  Longitudinal  Layer 

The  cells  of  the  longitudinal  layer  do  not  actively  participate 
in  slow  wave  generation.  The  evidence  obtained  from  the  work  described 
in  this  thesis  indicates  that  this  layer  acts  as  a  passive  conductor 
for  electrical  activity  generated  in  the  circular  muscle  layer.  The 
lack  of  a  measurable  phase  lag  across  this  layer  indicates  that  the  layer 
is  almost  purely  resistive. 

In  the  model  of  Figure  6.3  the  longitudinal  muscle  layer  is  repres¬ 
ented  by  the  resistors  R^  and  R  .  The  resistance  R^  represents  the 
coupling  which  exists  between  the  circular  and  longitudinal  muscle 
layers.  R^  represents  the  distributed  resistance  of  the  longitudinal 
layer.  Although  R^  is  shown  connected  to  a  separate  output  of  the 
monostable  oscillator,  it  is  possible  that  it  could  be  connected  to 
the  normal  output  coupling  the  oscillator  to  an  adjacent  oscillator. 

In  either  case,  R^  and  R^  are  factors  which  influence  electrical 
coupling  in  both  the  longitudinal  and  circular  directions. 
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Special  Layer 


Figure  6.3  A  model  showing  a  possible  arrangement  of  oscillators  in 
the  colon  wal 1 . 
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6.8  FUTURE  WORK 

The  quest  for  knowledge  is  an  unending  search.  As  a  consequence 
of  this  it  is  expected  that  the  research  described  in  this  thesis  will 
present  a  number  of  questions  to  be  dealt  with  in  future  work.  The 
more  interesting  of  these  questions  pertain  to  the  metabolic  mechanisms 
responsible  for  slow  wave  generation.  The  hypothesis  presented  in 
Section  6.6  should  be  tested.  This  can  be  accomplished  by  conducting 
intracellular  measurements  using  voltage  clamp  or  chemical  clamp 
techni ques . 

The  knowledge  gained  from  the  canine  colon  studies  has  significance 
only  if  it  can  be  related  to  the  electrical  activity  of  the  human  colon. 
At  present,  i ntracel 1 ul ar  studies  of  human  colon  specimens  have  not 
been  very  successful.  Using  the  procedures  and  equipment  developed 
for  the  canine  studies,  future  measurements  on  human  colon  tissue 
should  be  more  successful.  Procedures  should  be  developed  to  insure 
the  viability  of  human  specimens.  With  such  procedures  the  small 
number  of  human  specimens  will  be  optimally  used  and  procedures  for 
successful  cellular  penetrations  will  be  developed.  Then  the  knowledge 
gained  from  the  canine  studies  can  be  applied  to  the  study  of  normal 
and  abnormal  human  tissue  with  the  final  objective  being  the  treatment 


of  disease. 


. 
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APPENDIX  A-l  INTRODUCTION 

This  appendix  provides  detailed  information  on  each  of  the  major 
programs  developed  specifically  to  process  the  data  obtained  from  this 
research.  These  programs  permit  the  transfer  of  data  between  various 
computer  storage  devices,  display  and  plot  graphical  results,  and 
process  the  data.  The  data  processing  is  used  to  eliminate  noise 
and  to  provide  signal  statistics,  power  spectra,  and  the  correlation 
between  data  channels. 

On  many  large  computer  systems,  such  as  the  Amdahl  370  system, 
packaged  programs  are  usually  available  to  perform  the  functions 
described  above.  Because  of  the  difficulties  in  obtaining  high  speed 
analog-to-digital  conversion  on  the  university  Amdahl  370  system,  a 
Hewlett-Packard  21  MX/E  minicomputer  had  been  previously  obtained  to 
digitize  the  data  from  human  and  canine  colon  studies.  It  was  more 
convenient  and  cheaper  to  use  the  minicomputer  for  the  required  data 
processing.  Therefore  the  programs  described  in  this  appendix  were 
developed. 

Because  of  the  high  speed  and  random  nature  of  data  access  from 
the  computer  disc,  any  data  file  requiring  processing  is  initially 
transferred  from  magnetic  tape  to  the  disc. 

This  transfer  is  accomplished  using  program  FLOCT  described  in 
Appendix  A-2.  All  other  programs  obtain  their  input  data  from  the 
disc  and  many  place  their  processed  results  back  on  the  disc.  Programs 
TDSPT  and  PLOT,  described  in  Appendices  A-3  and  A-4  respectively,  are 
utility  programs  to  display  unprocessed  or  processed  results  in  graphical 
form  on  the  system  CRT  display  or  on  the  high  speed  teleprinter. 

Appendix  A-5  describes  programs  FILTR  and  DCAC.  Program  FILTR  applies 


. 


a  band  pass  digital  filter  to  the  data  file  on  disc  and  program  DCAC 
removes  the  dc  component  of  the  data  record.  Appendices  A-6  to  A-8 
describe  programs  PSPT,  PHSE,  and  XCORR  which  determine  data  power 
spectra,  phase  and  frequency  statistics,  and  cross  correlation  data 
between  channels.  Program  MAIN,  described  in  Appendix  A-9,  is  a 
master  program  which  permits  the  researcher  to  select  a  menu  of  any 
or  all  of  the  above  programs  for  application  to  a  particular  data 
fil  e. 

Each  of  the  following  appendices  provides  a  brief  description  of 
each  program,  a  list  of  input  parameter  and  subroutine  requi rements , 
a  flow  chart  of  the  more  complex  programs,  and  a  complete  listing 
of  each  program  and  its  required  subroutines. 
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APPENDIX  A-2 

PROGRAM  FLOCT  (£ile  LOCate  and  Transfer) 

Descri ption 

Program  FLOCT  is  designed  to  transfer  digital  data  files  from 
magnetic  tape  to  the  HP  21  MX/E  disc.  The  program  provides  the  user 
with  the  options  (1)  of  defining  the  record  by  tape  location  or  file 
name,  (2)  of  transferring  all  or  any  portion  of  a  data  file,  and  (3) 
of  reducing  the  effective  sampling  rate  by  transferring  on  every  n'th 
sample.  Prior  to  the  transfer  the  data  is  appropriately  filtered  to 
prevent  aliasing  errors.  The  filtering  is  accomplished  with  a  fourth 
order  recursive  digital  filter.  Two  forms  of  the  data  file  are  trans¬ 
ferred  to  the  disc  using  the  subroutines  DWRI  and  STPT.  Both  of 
these  subroutines  are  also  described  later  in  this  section.  Subroutine 
DWRI  writes  the  data  file  onto  disc  without  change  while  STPT  scales 
that  data  and  changes  it  into  a  form  suitable  for  output  to  the  CRT 
display.  A  flow  chart  of  program  FLOCT  is  provided  in  Figure  A-2- 1 
and  the  program  is  given  in  Program  Listing  A-2- 1 . 

Subroutines  Required: 

1 .  DWRI  (Disc  WRIte) 

2.  STPT  (STore  PloT) 

Required  User  Input: 

1 .  File  name  or  +N  or  -N 

a.  The  program  examines  the  first  six  characters  of  a  file  name 
to  locate  a  file. 
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b.  A  space  followed  by  +N  or  -N  will  advance  or  rewind  the  tape 
N  files.  The  space  cannot  be  omitted. 

c.  If  the  record  cannot  be  found  the  program  will  search  to  the 
end  of  the  tape,  rewind,  and  print  "FILE  NOT  FOUND".  The  user 
may  then  input  a  new  file  name  or  tape  location. 

2.  Following  "OK  TO  PROCEED" 

a.  Answer  "YES"  or  "NO".  Printing  "NO"  provides  the  option  of 
selecting  another  file. 

3.  Number  of  channels  in  the  data  file. 

4.  The  number  of  256  word  records  to  be  skipped  before  data  transfer 
be g i  ns  . 

a.  This  permits  transfer  of  any  portion  of  the  file. 

5.  The  number  of  256  word  records  to  be  transferred. 

6.  The  data  sampling  rate. 

7.  The  number  of  data  samples  to  be  skipped  during  the  transfer. 
(Process  every  N'th  sample.) 

a.  "N"  must  be  a  power  of  2.  If  it  is  not,  it  will  be  reduced 
to  the  nearest  such  power. 

Program  Output: 

1 .  File  name. 

2.  Actual  data  sampling  rate. 

3.  Number  of  256  word  records  transferred. 

4.  Error  messages : 

i .  Fi 1 e  not  found . 

ii.  End  of  file  (EOF)  encountered. 

iii.  Wrong  number  of  channels. 


iv.  Tape  read  error. 


i  08 


Figure  A-2-1  Flow  chart  of  FLOCT 
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Program  Listing  A-2-1  FLOCT 


PAUL  OutJl  (fIN4--Ph.Le.ASE  24177B —  JULY,  1971) 


0  0  0  1 

F  T  N  4  ,  l 

0  0  U  2 

0  0  03 
00  0  4 

0  0  05 

0  0  06 

0  0  07 

5 

0  0  08 

100 

0  0  09 
0010 

0  0  11 

1  0  1 

00  12 
0013 
0014 

') 

0015 

6 

0  0  16 
0017 
00  18 
0019 
0020 

0  021 

0  0  22 

21 

0  0  23 

22 

0  024 

0  0  25 

2  5 

1)  076 

0  07? 

0  028 

31 

0  0  29 
00(0 
0031 
(10  32 

0  0  3  3 

0  0  34 

24 

0035 

15 

0  0  56 

0  0  37 

1  0 

0038 

0  0  39 
004  0 

0  0  4  1 

0  0  42 

0  0  43 

1  1 

0  0  44 

0  0  45 

0  u  4  6 

25 

0  0  47 

0  0  48 
0049 

16 

0050 

1  7 

0  051 

104 

0  0  52 
0053 
0054 

105 

0  055 
0056 

106 

PROGRAM  )  I  QL I 

DIMENSION  IRAY  (i’u  Ul>  ,  JR AY ( 8 , 256 )  ,KRAY<256)  ,Y(256)  ,Yi(8,4)  , 

1  A<5)  ,D(5)  ,  INAM(3) ,  LAFjF ( 135)  , XI ( 8 , 8 ) ,  Y2( 8 ) , MAX (8,96) 

L  QIJIVAl  ENCE  <Y2(5),Y) 

CALL  F  XI  C(  17  ,  IF'TRK  ,  ILTRK  , 1SIZE) 

WRITE  (  1  ,100) 

FORMAT  ( / "  INPUT  FILL  NAME  UR  +  N  Ok-  N  TU  ADVANCE  OR  RIU1ND  N  FI! 
1  M ) 

READ< 1  ,101) ( I  NAM ( 1 )  ,  1-1 ,3) 

FORMA1 ( 3 A? ) 

IF( INAn(l) , NE.2FU  )CO  TO  7 
CO  TU  6 

I F (  I  NAM ( 1 )  . NE . 2H  )GO  TO  lb 
1NH-INAM(2)72 56 
INL=INAM(2)-INH*256 
IF  (  INI  .  EQ  .  4 OB )  GO  TU  21 
I NH- ( INH-60B ) ♦ 1 0 

I  NL  -  I  NL  60B 
IK INH+ 1NL)*2 
CO  TO  22 
IKINH-60H)*? 

IF  (  I  N AM (  1 )  .  EQ  .  21  h  )UO  IU  2  3 
I- -I  -1 

CAI  I  P  TAPE  (  1  2 , 1  ,  0  ) 

CAl  I  F  XI  L<  1  ,  1 141:,  |  RAY  ,211411) 

IF  (  II  UT  (  12)  .  1.  T  .  0  )Gll  TO  1004 
IF  ( I WR  I)S  ( 1  2 )  .  EQ  .  36  )  CO  TO  2r> 

CAl l  P TAPE ( 12,-1 , -2 ) 

CAI  I  E  XF.C  (  1 ,114B,  IRAY  ,2048) 

II  (  I  URDU  (.12)  .EQ  .  36)  GO  IU  2S 
CAl L  P  TAPE (12,2,0) 

CD  TO  31 

CAI  L  P  TAPE  (  12,2,0) 

CAl  L  F  XF  C  (  1  ,  1 1  4  B  ,  IRAY  ,2040) 

IF  < IEO 1  ( 1 2)  . LT . 0  >GO  TO  10  04 
II  (  I WRD.J (12)  .EQ.  36)00  TO  1  1 
CAI  l  P  T  AP  E  (  12,-1  ,  -  2  ) 

CAl  L  EXPl'K  1  ,11  4B  ,  IRAY  ,2048) 

IF  ( I IJRDS ( 1  2)  .  EQ  .  36  ) CO  TO  11 
CAI  I  P  I  APE  (12,2,0) 

CD  IU  10 


I  NAM ( 1  ) >CU 
I  NAM  (  2  )  >G(J 
I  NAM  (  3  )  )i;i) 


IO  17 


I F (  IRAY  < 1 )  . NL 
IF  (  1  R  A Y  <  2  )  .  Nl 
II  (  IRAKI)  .NL 
DO  16  1  - 1  ,36 
L  ABF ( 1  )  -  I R  AY  < 1 ) 

I  F  <  IRAY  <  I  )  .  EG) .  0  )GU 
CONI INUL 
OR  1 TE ( 1  ,10  4) 

FORMAT </"*FILENAMI  **"/> 
CAI  I.  f  xf  L(2 , 1  ,  IRAY  ,  I  ) 
WRITE  ( 1  .105) 

FORMAT ( / "OK  TO  PROCEED? 
R  F  A  I)  (  1  ,106)  IANS 
FORMA  T ( A2 ) 


T  0 
TO 
IO 


24 

24 

24 


Program  Listing  A-2-1 


FLOCT  (cont.) 


PAGE  0  0  02  FI  OCT 


(FTN4--RELE.ASE  24177B — JULY ,  1971) 


t 

i 


0  0  57 

IF  ( IANS  .  fcU  •  2HYE  >l.ll  TO  16 

OGSfl 

CAl  L  PTAPE  (12,-1  ,0) 

0059 

t;o  TO  5 

0  0  AO 

18 

URITE( 1  ,  1 07) 

0  06  1 

1  07 

FORMAT (/"NO.  OF  CHANNELS? 

"> 

0  0  6? 

R  E  AD  <  1  ,  ior)nc:ha 

0063 

1  013 

FURMAM  11  ) 

0  0  6  4 

l A8F ( 37 ) -NCHA 

0  0  65 

WRITE  (1,110) 

0  0  66 

109 

FORMAT (/"SKIP  N  RECORDS? 

) 

0067 

REALM  1  ,*  )  NRECS 

0068 

IF  (  NRF.CS  .  GT  .  96)NRECS  =  96 

0  069 

WR1 IE( 1 , 109) 

0  070 

1  1  0 

FORMAT (/  "NO .  OF  RECORDS? 

) 

0  071 

RE  AD( 1  , * )NSKP 

0  0  72 

WRITE .  (  1  ,111) 

0  073 

1  1  1 

FORMAT (/"DATA  SAMPI1NG  RATE 

PER  MINUTE?  _") 

0074 

R  E  A  D (  1  , * )SR 

0075 

WRITE  < 1  ,112) 

0  0  76 

112 

FORMAT (/"PROCESS  EVERY  N'TH 

SAMPLE?  _." ) 

0077 

READ (  1  ,  *  )NTM 

0  070 

TF ( NT  H . L  r . 64 ) GO  TO  19 

0079 

NTH  =  6  4 

0080 

1  9 

I F ( NT  H . G 1  . 1 ) GO  10  33 

0081 

NT  HM 

0082 

CO  TO  86 

0  0  63 

33 

I T  W  O  =  0 

0  084 

2  0 

NTH-  N 1  M/2 

0  (i  85 

itwom t  wn+i 

0  08t> 

1 F ( N I H . CE . 2 ) CO  10  70 

0  087 

NTH-2**1 TWO 

0  0  88 

26 

I ASR-SN/NT  H 

0  089 

ASR-IASR 

0  0  90 

WRITEd  ,  1 1  4  )  A  S  R 

0091 

1  14 

FORMAT(/*,ACTlJAL  SAMPt  INC  RATE  I S  "  ,  X  ,  F  7 , 1  ,  "  CYCI.F.S/M1N 

0  092 

L ABF ( 38 ) =1 ASR 

0  093 

CALI  F  Xt  C(  1  ,  1146 , 1KAY  ,  0  ) 

0  094 

I  F  (  N!>K  P  ,  E  Q  .  0  )  CO  TO  39 

0095 

DO  27  Id  ,  NSK  P 

0096 

CAl  1  1  XI  C(  1  ,  1 1  4Ft ,  1  R AY  ,  0  ) 

0  0  97 

IF  (  IF  HI  (  12)  ,LT  .  0  )CD  TO  1001 

0  098 

27 

CONTI  Nlll 

0  0  99 

39 

DO  36  K-  1,256 

0100 

Y  (  K  )  -  0 

oi  n  i 

36 

CONT l NOE 

0102 

DO  41  KM  .32 

0  10  3 

X  1  (  K  )  =  ()  . 

01  04 

Y  1  (  K  )  -  0  . 

0105 

4  1 

CONT  I  NO! 

0106 

DO  4  5  KM  ,8 

0  107 

DO  45  1=1,96 

0106 

MAX ( K ,1)0 

0  109 

45 

CON  1 I  NOE 

0110 

SR  =  SR /6 0  . 

0111 

FC  =  SR / < N 1 H*  3  ) 

0112 

CAl  t  MITF  (SR,  EC,  A, 
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Program  Listing  A-2-1  FLOCT  (cont.) 


PACK  00  03  FLOCl  (FTN4~~NF.I  EASE  24  177b  —  Jill.  Y ,  1971) 


0113 

DD  2D  1=1,  NKECF. 

0114 

CALI.  1  XE  C(  1  ,  1 1  4K ,  IK  AY  ,20  48) 

0  115 

IF ( I  EOF ( 1 2 )  . L  T , 0  ) GO  18  1001 

0  116 

IF  ( I URDR ( i 2 )  ,  EQ  .  0  ) GO  T(l  1002 

0  117 

IF  (  1 WR  DS ( 1 2  > /256  *256 , NF  ,  1 WRDS ( 12 )  ) CO  TU  1002 

U  1  18 

NC-IL)RD5(12)/256 

0119 

IF (  NC  .  NE" .  NCHA )  GO  TO  1003 

l)  1 2  0 

29 

OO  30  J=1,NCHA 

0  1  2 1 

DO  30  K - 1 ,256 

0  122 

J  R  A  Y  (.1  ,K  )  -  I  R  A  Y  (  ( K  - 1  )  *NC  +  J  ) 

0  123 

DO  3  0  I  =5  ,  a 

0  124 

XI  (  J  ,  1  )  =  J  R  A  Y  C  .1  ,  L  4  ) 

0  1  25 

30 

CONT 1 NUF 

0  1  26 

DO  28  J - 1 , NCHA 

0127 

DO  43  L= 1  ,4 

0  1  28 

Y  2  C  l.  >  --  Y  1  ( J  ,  L ) 

0  1  29 

43 

CONTI NUF 

0130 

DO  34  K  =  1  ,  4 

0131 

Y ( K ) - A ( 1 ) # JR AY  <  J  ,K)+A(2)+X1(J  ,3+K)+A(3)*Xl(J, 2+K )  + 

0  1  32 

1 

A<  4  )  *X1  <  J  ,  1  +  K  >+A(5)*Xi  (.1  ,K  )+B(  1  >*Y2(3-«K  )  + 

0133 

L. 

1  B(2)*Y2<2-*K  >  +  BC3)*Y2<  i  +  K  )+B<4)*Y2  <l<  > 

01  54 

34 

CONTINUE 

0  1  35 

DO  35  K  =5 , 256 

0136 

Y(K )=A<1 >*JRAY( J  ,K  >+A<2)*JkAY< J ,K-i >+A(3)*JRAY( J ,K-2>  + 

0  137 

1 

A<  4 ) *JKAY<  J ,K-3 ) +A<5)*JRAY( J ,K-4 )+B( 1 )*Y(K- 1 >  + 

0138 

’  B(2)*Y(K-2)+B(3)*Y(K-3)+B(4)*Y(K-4) 

0  1  39 

35 

CONI INUE 

0  140 

DO  47  K - 1  ,256 

0141 

IMAX-  I  F  l  X ( Y ( K  )  ) 

1142 

I  MAX- I  ADS ( IMAX ) 

0143 

l  hAX- MAX ( J  ,  I  ) 

0144 

IF  ( IMAX  .  LE  .LMAXXJD  IO  47 

0145 

MAX ( J  ,  I  )  = I  MAX 

0146 

47 

CONT INUE 

0147 

IMAX=MAX ( J , I ) 

0148 

4 1> 

DO  4  0  K  =  1  ,  4 

0  149 

XI  (  J  ,K  )  -  .T R  A  Y  (  J  ,  252 +  K  ) 

0  150 

Yi  ( J  ,K  )  =Y ( 252  +  K ) 

0151 

40 

CONT I  NOE 

0  152 

DO  38  K - 1 ,256/NTH 

0  153 

KRAY  (  K  )-  -  (Y(  ( K  - 1  )  #NTH  + 1  )  +  0  .  5  ) 

0  154 

38 

CONTINO) 

0155 

L  N  O  T  H  -  2  5  f  >  /  N  1  H 

0  1  5  6 

t 

CAI  L  DUR  I  (  IF  TRK  ,1  NO  I II  ,KRAY  ,  1  ,  J  ,  1  ) 

0157 

CALI  STPT< IF TRK , KRAY , LNGTH , I , J , IMAX) 

0  158 

28 

CONT I  NOE 

0  159 

CAl.L  PT  APE  ( 1 2  ,  -2 , 0  ) 

0160 

NRECS=NR ECS/NT H 

0  161 

L.AF<F  (  3y  )  -NR  ECS 

0162 

DO  44  K - 1  , NR  ECS 

0  163 

L  ABF (K+39)=MAX( 1  ,K) 

0  1  64 

44 

CONT  INOI 

0  1  65 

CALL  DUR  HIE  IRK,  0,IAUF(1), 0,0,1) 

0166 

WRITE  < 1  ,  1 1  3  )  N  R  F  C  S 

1 1  67 

113 

FORMAT  (/"NO.  OF  256  SAMI'l  t  RECORDS  TRANSFERt  D  15  "  ,13) 

0  168 

STOP 

216621 


U  1  69 

o  r'o 

0171 

0172 

0173 

0  I'M 

0  175 

0  176 

0  177 

017B 

0  1  79 

OlfaO 

0  181 

0182 

0183 

018  4 

0  1  85 

0186 

0  1R7 

0  1  88 

0189 

(1190 

0191 

0  192 

0  19  3 

0  194 

0  195 
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Program  Listing  A-2-1  FLOCT  (cont.) 


PAGE  UUU4  FLUCT  (  F  1  N 4  —  RELEASE  24177b — JULY  ,  1971) 


i ono  uk i n  < i  , 900 > 

90  0  FORMAT (/ "FILE  NO  I  FOUND!!") 

GO  TO  5 

Him  WRITEC  1  ,901) 

901  FORMAT ( 7 "E  OF  ENCOUNTERED,  CHECK  NO.  OF  RECORDS!!") 

CAL  l  P  TARE  (  12  ,  —  3 ,  0  ) 

GO  TO  5 

1  0  02  WR IT E  M  , 9  02) 

902  FORMAT </ "*TAPF  FUROR"/) 

CAL  L  P  T  APE  (  1  2,  --2, 0  ) 

CO  TO  5 

1  003  NCHA  =  IIJRDS(  12  >/2S6 
WRITEC 1 , 9  0  3 ) N  C  H  A 

903  FORMAT  (  /  "WRONG  NO.  OF  CFIANNE  IS!  !  "/  "NO  .  OF  CHANNELS  IS  ",1 1  ) 
UR  ITEM  ,904) 

904  FORMA  T ( /  "DO  YOU  WISH  TO  CORRECT?  ") 

R  E  A  D (  1  ,  1 06) IANS 

IF ( IANS . EQ . 2HN0 ) GO  TO  29 
UR  I  1F( 1  , 1 07) 

Ri  AD< 1  .1 OBINCHA 
LAliF  (  V/  )  -NCI  IA 
CO  TO  29 

1  0  04  WR] TE< 1  , 907) 

907  FORMAT ( /"END  OF  TAPE") 

CAI  L  FXI  C<  3,414*0 

GO  TU  5 

END 
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APPENDIX  A-2  (cont. ) 

SUBROUTINE  DWRI  (IFTRK,  LENGTH,  NAME,  NUM,  ICHA,  I  FILE ) 

Descri ption 

Subroutine  DRW  I  ([)isc  WRI  te)  is  a  subroutine  designed  to  write 
data  files  from  memory  into  one  of  three  disc  files.  In  addition  a 
header  record  is  created  on  disc  containing  the  file  name,  number  of 
channels,  actual  sampling  rate,  and  number  of  data  records  contained 
in  the  file.  The  three  disc  files  numbered  1,  2,  and  3  are  created  to 
store  input  data,  processed  data,  and  di spl ay-formatted  data  respec¬ 
tively.  The  program  listing  and  disc  organization  associated  with 
DWRI  is  shown  in  Figure  A-2-2.  Each  disc  file  consists  of  96  disc 
tracks  with  the  data  stored  in  order  of  channel  number.  The  maximum 
capacity  of  each  of  the  three  disc  files  is  eight  channels  of  data 
each  containing  96  256-word  records. 

Subroutine  input  parameters: 

1.  IFTRK,  the  first  available  free  track  on  disc. 

2.  LNGTH,  the  length  of  each  data  record. 

3.  NAME,  the  name  of  the  memory  array  containing  the  data  file. 

4.  ICHA,  the  number  of  a  particular  channel  to  be  written  onto  disc. 

5.  IFILE,  to  specify  disc  file  1,  2,  or  3. 
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FAGF  0  u  «J  1 


(FTN4-- RELEASE  24177B — JULY,  1 9 7 1 > 


0  0  0  1  F  T  N  4  ,  L 

0  0  02  SUBROUTINE  DWK  I  (  IFTRK  ,LNGTH  ,NAME  ,  NUM  ,  ICHA  ,  IFILE ) 

0  0  03  DIMENSION  NAME (  1  ,  1 ) 

0004  IF (L NOIH . NE . 0 ) GO  TO  1 

0  0  OS  L -  1 3S 

0  0  0 ITOF=IFTRK+( IFILE- 1 ) *33 

0  0  07  IOFF  =  0 

0  0  08  I  UR DS  =  0 

0009  GO  10  2 

00  10  1  I  NLIM-NUM-  1 

00  11  L  =  L.NG'IH 

0012  I TNUM-L* 1 NUM 

0013  1  URDS-MOD  (  I  TNI  IM  .  1  28  ) 

0014  I RMDft-lTNUM- IURDS 

001S  ISTRS=IRMDR/128 

00  16  I  OF F  =  MOD< I  SIRS , 48) 

0017  ITOF  =  (  ISTRS-IOt  !")/48 

0  0  18  II  OF  = IT  OF  +  1  -M  ICHA-1 ) *4  + 1 F  TR K  +  ( I F I LE- i  ) *33 

0  0  1  9  2  CALI  BURN  (NAME  ,L,  2,  ITOF  ,  IOFF  ,IURI)Ei) 

0020  RETURN 

0021  END 


**  NO  ERR  URL* 


!  RA  K 

SECTuR 

1  FILL 

CHANNEL 

1 1  TKK 

0 

37 

38 

39 

1 

FI LE  NAME  0*36 

37  =  NO.  OF  CHANNELS 

38  =  ACTUAL  SAMPLING  RATE 

39  =  NO.  OF  RECORDS 

IF  IKK  t  1 

0 

] 

3 

- 

A 

- 

5 

2 

9 

3 

13 

4 

29 

8 

IE'l'KK  +  31 

0 

2 

Figure  A-2-2  Program  listing  of  DWR I  and  table  showing  disc  file 

organi zation . 


APPENDIX  A-2  (cont.) 


SUBROUTINE  STPT  ( I FTRK,  NAME,  LNGTH,  NUMR,  ICHA,  MAX) 

Descri ption 

Subroutine  STPT  ( STore  PloT)  scales  the  amplitude  of  the  values 
in  a  data  file,  generates  a  new  file  suitable  for  fast  display  on  the 
computer  CRT,  and  stores  this  new  file  on  disc  file  3.  Because  the 
vertical  span  of  the  CRT  is  only  256  points  and  because  a  maximum  of 
8  channels  may  be  simultaneously  displayed,  each  data  file  is  scaled 
to  have  a  maximum  amplitude  of  32.  To  display  a  data  point  on  the 
system  CRT  a  16  bit  word  is  output  to  the  D/A  converter.  The  most 
significant  eight  bits  of  this  word  define  the  vertical  position  of  a 
display  point  and  the  least  significant  8  bits  define  the  horizontal 
position.  Therefore,  after  the  data  values  have  been  scaled  the 
amplitude  information  is  shifted  into  the  most  significant  eight  bits 
and  values  correspond!' ng  to  the  location  of  the  data  points  in  the 
256-word  record  are  added  to  the  least  significant  8  bits.  The  scaled 
display- formatted  data  file  is  stored  in  disc  file  3  using  subroutine 
DWRI.  The  subroutine  is  given  in  Program  Listing  A-2-2. 

Subroutine  input  parameters: 

1.  I FTRK,  the  first  available  free  track  on  disc. 

2.  NAME,  the  name  of  the  memory  array  containing  the  data  file. 

3.  LNGTH,  the  length  of  each  data  record. 

4.  NUMR,  the  number  of  a  particular  record  to  be  stored  on  disc. 

5.  ICHA,  the  number  of  a  particular  channel  to  be  stored  on  disc. 

6.  MAX,  the  maximum  amplitude  in  the  data  file. 


. 


' 


APPENDIX  A-2  (cont. ) 


Subroutines  required: 


1 .  Subrouti ne  DWRI 


Program  Listing  A-2-2 


STPT 


CuOl 
0  0  02 
0  0  03 
0  0  04 
OOOS 
0  006 
0  0  07 
0  00B 
0009 
OOi  0 
0011 
0012 
00  13 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 


PAGE  UOOl  ( FTN4--RELEASE  24177B  —  JULY, 


FTN4  ,  L 

SUBROIU]  NE  S  I PT ( IFTRK , NAME ,  L.NGTH  ,  NIIHR  ,  ICHA,HAX) 
DIMENSION  NAME ( 1 ) 

MX=MAX 

IF (MAX . EQ . 0 )MX=32767 
N-256/LNG1 H 
NN=NUMR-i 
I P  =MOD ( NN , N ) 

DO  10  1=1 , LNG (H 
I R Y  =  NAME (  I) 

NAME ( J  )=FLOAT( IRY)*32/MX 
IR Y  =  NAME  < I ) 

IF ( IRY  .LE .31 ) GO  TO  2 
NAME (  I )  =  3 1 

2  IF ( IRY . GT .-31 >GO  TO  1 
NAME ( I >=-31 

1  NAME ( I ) =NAME ( I ) *256+1 +1 P *LNGTH-1 
10  CONTINUE 

CALL  D  UR  I  (  IF  IKK  ,  L.NGTH  ,  NAME  ,  NUMR  ,  ICHA  ,  3  ) 

RETURN 

END 


1971 ) 


**  NO  ERRORS* 


APPENDIX  A-2  (cont.) 


SUBROUTINE  DREA  (IFTRK,  LENGTH,  NAME,  NUM,  ICHA,  I  FILE) 

Descri ption 

Subroutine  DREA  (Disc  REAd)  transfers  data  records  from  disc  to 
memory  arrays.  It  assumes  the  disc  file  has  been  formatted  by  DWRI 
as  described  above.  The  main  difference  between  DREA  and  DWRI  is  th 
direction  of  data  flow.  The  subroutine  is  given  in  Program  Listing 

A-2-3. 


Subroutine  input  parameters: 

IFTRK,  LENGTH,  NAME,  NUM  ICH/1  and  IFILE  as  described  for  sub¬ 


routine  DWRI. 


Program  Listing  A-2-3 


DREA 


0  0  0  1 
0  0  02 
0  0  03 

100  4 

0  0  0S 
0  0  0b 
0  0  07 
0  0  08 
0  009 
0010 
0011 
0012 
0013 
0014 
OOlS 
001b 
0017 
00  IB 

0  0  19 
0  020 
0021 


PAGE.  00  0  1  (FTN4  —  RELEASE  24 177B--JULY  ,  1971) 


FTN4 , L 

SUBROUTINE  DREA ( IFTRK , LNGTH , NAME , NUh , ICHA  ,  IF  ILE  ) 
DIMENSION  NAME ( 1 ) 

IF  <  L  NGTH . NE . 0 )GO  TO  1 
L  =  1  35 

ITOF  =  IFTRK  +  (IFIL.E-i  >*33 
IOFF  =  0 
IURDS=Q 
GO  TO  2 

1  INUM  =  NlJM-i 
L=t.NOTH 

I T  NUM  =  L* I NUM 
IURDS=MOD(  ITNIIM  ,128) 

IRMDR=ITNUM-I WRDS 
ISTRS-IRMDR/128 
IOFF=MOD( ISTRS  48) 

I T OF MISTRS-JOrP  )/48 

ITOF  =  ITOF+l  +  (  ICHA  -1  >*4  +  IFTRK  +  (  IFILF.-1  )*33 

2  CALL  BREAD( NAME  ,  L  ,2  ,  ITOF  ,  IOFF , IURDS) 

RETURN 

END 


**  NO  ERRORS* 


17C 


APPENDIX  A-2  (cont.) 

SUBROUTINE  RDPT  (IFTRK,  NAME,  NUMR ,  ICHA,  LOC) 

Description 

Subroutine  RDPT  (Rea£  PloT)  transfers  di spl ay- formatted  data 
from  disc  to  memory  arrays.  The  length  of  the  display  record  is  fixed 
at  256  words.  RDPT  utilizes  subroutine  DREA  to  accomplish  the  disc 
read  operation  and  then  modifies  the  most  significant  8  bits  as  a  function 
of  vertical  location  (LOC).  Each  data  channel  may  be  displayed  in  one 
of  eight  vertical  positions  depending  on  the  value  assigned  to  parameter 
LOC.  A  value  of  1  positions  the  data  channel  at  the  top  of  the  CRT 
while  8  positions  the  channel  at  the  bottom.  The  subroutine  is  given 
in  Program  Listing  A-2-4. 

Input  parameters: 

1.  IFTRK,  NAME,  NUMR,  and  ICHA  as  previously  defined  for  STPT  and 
DWRI . 

2.  LOC,  location  of  the  data  channel  on  the  CRT. 


. 
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OOOi 
0  002 
00  03 
0  0  04 
0  0  05 
0  0  06 
0  007 
0  0  08 
0  0  09 


**  NO 


Program  Listing  A-2-4  RDPT 


PAGE  OUOl  ( FTN4 — RELEASE  24  177B— JULY  ,  1971) 


FTN4 , L 

SUBROUTINE  RDPT  ( IFTRK , NAME , NUMR , ICHA,LOC> 
DIMENSION  NAME ( 1 ) 

CALL  DREA( IFTRK  ,  256  ,  NAME  ,  NUMR  ,  ICHA , 3 ) 

DO  10  1=1,256 

NAME ( I  )=NAME(I)  +  (223-(L0C-l)*64)*256 
10  CONTINUE 
RETURN 
END 


ERRORS* 


APPENDIX  A-3 


PROGRAM  TDSPL  (Time  DiSPLay) 

Descri ption 

Program  TDSPL  displays  from  one  to  eight  channels  of  data  on  the 
HP  21  MX/E  system's  CRT  and  facilitates  either  low  speed  or  high  speed 
scanning  of  the  complete  data  file.  The  sense  switches  control  the 
program  mode  of  operation.  If  sense  switch  15  is  set  the  record  will 
slowly  move  across  the  CRT.  If  sense  switch  14  is  set  the  file  is 
rapidly  scanned  on  a  256-word  record-by-record  basis. 

The  data  to  be  displayed  must  be  previously  placed  on  disc  by 
FLOCT.  Additional  information  about  the  program  is  provided  by  the 
flow  chart  of  Figure  A- 3-1  and  Program  Listing  A- 3-1 . 

Required  subroutines: 

1.  Subroutine  GRID  displays  the  grid  on  the  CRT.  Program  Listing 
A-3-2  is  a  listing  of  subroutine  GRID. 

2.  Subroutine  RDPT 


3.  Subroutine  DREA 


. 
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Figure  A- 3-1  Flow  chart  of  TDSPL. 
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0012 

0013 
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001S 

0016 
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0019 

0  020 

0  021 

0  022 

0  02  3 

0  024 

0025 

0  0  26 

0  0  27 

0028 

0029 

0  0  30 

0  031 

0032 

0033 

0034 

0035 

0036 

0  037 

0  0  38 

0  0  39 

0  040 
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0  0  44 

0045 
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0047 
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0052 
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0055 

0056 
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Program  Listing  A-3-1  TDSPL 


PAGE  0  0  0  1  (F TN4  — RELEASE  241778  —  JULY  ,  1971) 


FTN4 , L 

PROGRAM  TDSPL 

DIMENSION  JR AY (256, 4) ,IGRID(2480  )  ,IRAY(256) ,LABF(13S)  , ICH(4) 
1  , K  R  A Y  (  256 , 4  ) 

EQUIVALENCE  ( I GR ID ( 1 457 )  ,  JR  AY  ) 

CALL  E  XEC ( 1 7 , I FTRK  ,  ILTRK  ,  IS1ZE) 

CALL  GR ID ( IGR ID  ) 

CALL  DREA< IFTRK  ,  0  ,LABF  ,  1  , 1  , 1 ) 

1  DO  10  1=1,4 
I CH ( I ) =  0 

10  CONTINUE 

DO  11  1=1,1024 
JR  A  Y  <  I  )  =  0 

11  CONTINUE 

WR I TF( 1 ,100) 

100  FORMAT (/"INPUT  THE  REQUIRED  FOUR  CHANNELS!  ") 

READ (  1  ,*) ( I CH ( J )  ,J  =  i ,4) 

NRECS=LABF (39) 

DO  26  1  =  1  ,4 
I  CHA=  I  C.H  ( I ) 

IF  ( I CHA  .  EG) .  0  )  GO  TO  15 
26  CONTINUE 
15  NCHA- I ~ 1 

I.TH-1  456  *  NCHA#256 
IFST  =  0 

CALL  EXEC  (  2 , 0  2  0  1 1 1 B  ,  I  GR  1 D  ,  l.TH ) 

DO  20  1=1 , NRECS-i 
14  IF ( IFST , GE , 0 ) GO  TO  12 
CALL  DElAY(iOOO) 

CALL  E  XEC ( 3 , 1 1 b ) 

CALI  EXEC ( 2, 0  20 11  IB, IGR ID, LTH) 

12  DO  21  J=i , NCHA 
IC= ICH ( J ) 

CAI  L  RDPT(  IFTRK  ,  IRAY  ,  I  ,  IC  ,  J  ) 

DO  25  K  =  1 ,256 
JR AY ( K  ,  J)=IRAY(K ) 

25  CONTINUE 

21  CONTINUE 

IF ( NR  ECS . EQ . 1 ) GO  TO  30 
DO  19  J= 1 , NCHA 
ICHA=ICH( J ) 

N-  I  ♦  1 

CALL  RDPT(IFTRK , IRAY ,N , ICHA , J ) 

DO  19  K  =  1 ,256 
K R  AY  (  K  ,  J  )  =  IR  A Y  (  K  ) 

19  CONTINUE. 

30  DO  24  L.~i  ,256 

22  I ST= I SSW (15) 

IFST= I SSU ( 14) 

I F (  IFST . LT . 0 )GO  TO  20 
IF  (  1ST  .I  T  .  0  )GO  TO  22 
DO  23  J  =  2 ,256 
DO  23  K=i , NCHA 
J  R  A  Y ( J  ~  1  ,  K )  =  JRAY ( J  ,K >-l 

23  CONTINUE 

DO  24  K  =  1  , NCHA 


- 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0  067 
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Program  Listing  A-3-1  (cont.) 


PAGE  0002  TDSPL  (FTN4  —  RELEASE  24177B— JULY, 


JR  AY ( 256 , K )  =K  R  AY  ( 1. ,  K  )  -  L +256 
24  CONI 1NUE 
20  CONTINUE 

WRITE< 1 , 102) 

102  FORMAT (/"DO  YOU  WISH  TO  REPEAT?  _"> 
READ ( 1 , 103) IANS 

CALL  EXEC (3,1  IB) 

103  F OR HAT ( A2 ) 

IF (IANS.EQ.2HYE)G0  TO  1 

STOP 

END 


1971  ) 


Program  Listing  A-3-2  GRID 


PAGE  00  0  1  ( FTN4--  RELEASE  24177B--  JULY,  1971) 


0  00  1 

FTN4 ,1. 

0  0  02 

SUBROUTINE  GRID(IGRID) 

0  0  03 

DIMENSION  I GR I D ( 1 396 ) 

0  0  04 

1 

DO  10  1=1,256 

0  0  05 

1GRID ( I  )  =  ( 1-1 ) *256+ 1 5 

0  0  0b 

10 

CONTINUE 

00  07 

DO  11  1=256,316,4 

00  08 

DO  11  J= 1 , 4 

0009 

IGR I D ( I + J ) = ( I -256) *1024 +12+ J 

00  10 

11 

CONTINUE 

0011 

DO  14  K=1 ,4 

0012 

Ii=32i+(K-1 )*268 

0013 

12=11+239 

0014 

DO  12  1=11,12 

0015 

IGR I D ( I >  =  ( 223-64* ( K - 1 ) ) *256  +  1 -305- ( K-i ) *268 

00  16 

12 

CONTINUE 

0017 

I1=S60+(K-1 )*268 

0018 

12=11+28 

0019 

DO  13  1=11,12,4 

0  02  0 

DO  13  J=1 , 4 

0  021 

IGR ID ( I  +  J  )  =  (220-64*(K-l )+J)*256  +  B*( I -558- (K-i ) *268) +32 

0022 

13 

CONTINUE 

0023 

14 

CONTINUE 

0024 

RETURN 

0025 

END 

**  NO 

ERRORS* 
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APPENDIX  A-4 

SUBROUTINE  PLOT  (X,  Y,  YG,  LABEL,  NPLOT,  NTR,  IT) 

Descri ption 

Subroutine  PLOT  produces  an  X-Y  plot  on  the  high  speed  printer 
of  from  one  to  four  data  arrays.  The  resolution  of  the  plot  is  16  units 
per  inch  in  the  horizontal  direction  and  12  units  per  inch  in  the 
vertical  direction.  The  span  of  the  plot  is  50  by  100  units.  When  a 
single  channel  of  data  is  to  be  plotted  the  subroutine  provides  the 
option  of  a  normal  plot  or  a  bar  graph.  A  flow  chart  of  subroutine 
PLOT  is  provided  in  Figure  A-4- 1 .  A  typical  output  for  one  data  array 
is  shown  in  Figure  A-4-2  while  a  typical  bar  graph  is  shown  in  Figure 
A-6-2.  Subroutine  PLOT  is  given  in  Program  Listing  A-4- 1 . 

Required  input  parameters: 

1.  X,  an  n  x  4  array  containing  the  horizontal  coordinates  for  all 
data  points. 

2.  Y,  an  n  x  4  array  containing  from  1  to  4  sets  of  vertical  coordinates. 

3.  YG,  the  maximum  vertical  amplitude  of  the  plot.  If  this  value  is 
zero  the  plot  is  scaled  such  that  the  maximum  vertical  amplitude 
corresponds  to  the  maximum  data  amplitude. 

4.  LABEL,  an  array  containing  up  to  25  alphanumeric  characters  con¬ 
stituting  the  vertical  label  for  the  graph. 

5.  NPLOT,  the  number  of  data  points  in  the  arrays  to  be  plotted. 

6.  NTR,  the  number  of  traces  to  be  plotted. 

7.  IT,  a  parameter  which  determines  if  a  single  data  array  is  to 
produce  a  normal  or  a  bar  graph. 


i  .  4  , 


APPENDIX  A-4  (cont.) 


If  IT  =  <j>,  a  normal  graph 
If  IT  =  1 ,  a  bar  graph  i s 


is  produced, 
generated . 
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Figure  A-4-1  Flow  chart  of  PLOT 
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Figure  A—  4—1  Flow  chart  of  PLOT  (cont.). 
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Figure  A-4-2  Typical  single  channel  output  of  subroutine  PLOT 
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OOOi 

0  0  02 

0  0  03 
00  04 

FTN4  , 

L 

SUBROUTINE  PLOT  (X  ,  YY  ,  YG  ,  LABEL  ,  NPLClT  ,  NTR  ,  IT  ) 

DIMENSION  X( 128)  ,YY( 128 , 4  )  ,Y ( 1 00 ,4) ,IL( 100 ) , LABEL (25) ,18(5) ,L1 (2 
DO  25  J=i ,NTR 

0  0  OS 

DO  25  1=1,100 

0  0  0  6 

1 1  =  1 FIX ( NP  LOT  #1/10  0 .  ) 

0  0  07 

Y ( I , J ) =YY ( I 1+ ( J-i ) *NPLOT ) 

0  0  0  8 

25 

CONTINUE 

0  0  0  9 

DO  26  1=1,100 

0010 

II=IF IX(NPLOT*I/100 . ) 

0  011 

YY ( I  )  =X ( 1 1 > 

0  0  1 2 

IF ( 1  . LT . 100 . AND . I . NE . 1 ) YY ( I ) =YY ( I )  +  ( X ( I I  +  i ) -X ( 1 1 ) ) * ( NPl  OT* 

0  0 1 3 

.1  1/1  00.  -II) 

0  0)4 

26 

CONTINUE 

0015 

DO  27  1=1,100 

0  016 

X ( I )=YY ( I ) 

0  017 

27 

CONI INUE 

0  0  1 8 

UR  ITS:. ( 1  , 101 )  15446B  ,  2H1 1 

0  0  1 9 

WRITE. (4  , 1 12 )  1HU 

0  0  2  0 

WRITE.  ( 1  ,101)  15446B  ,2.1418 

0  02.1. 

URITEd  ,  112)1  HV 

0  0  22 

ITAB=15461B 

0  023 

IS(  .1  )=lHx 

0  0  24 

IS(2)=lHo 

0  02  S 

I S  ( 3 )  - 1 H+ 

0  026 

I S  <  4 )  =  1 H  " 

0  02.7 

IS ( 5 ) =7740  0B 

0  0  28 

WRITE ( 1  , 1 0  0 ) II AB , I TAB , J TAB 

0  0  2.9 

10  0 

FORMAT  (  2. OX  ,  A2 , 4X  ,  A2 , 6X  ,  A2. ) 

0  0  3  0 

WRITEd  ,  10i)15446B,2H10 

0  031 

WRITEd  ,  1 1 2 )  1 H I) 

0  0  32 

10.1. 

FORMAT < 2A2,  "  ") 

0  033 

WRITEd  >101)  15446B  ,  2Hk2 

0  0  34 

WR I TE (1,112)1 HS 

0  035 

YMAX=~  1  .  E.+36 

0  036 

ym:i:n=i  .  e+36 

0  0  37 

XMAX=-i .E+36 

0  038 

XMIN=i . E+36 

0  0  39 

DO  50  1=1,100 

0040 

IF ( X ( I ) ,GE. XMAX)XMAX=X( I ) 

0  0  41 

IF (a(I) . LE . XMIN ) XMIN=X ( I ) 

0042 

5  0 

CONTINUE 

0  0  43 

DO  51  I=i,100*NTR 

0  044 

I F ( Y ( I ) . GE. YMAX)YMAX=Y(I) 

0  0  45 

IF  (Yd)  .  LE.YMIN)YMIN=Y(I) 

0  0  46 

51 

CONTINUE 

0047 

IF ( YG . GT . YMAX ) YMAX=YG 

0  0  48 

NMIX=0 

0049 

NP  LX  =  0 

OOSO 

XMX=XMAX 

0  05.1. 

XMN=XMIN 

0  0  S  2 

XM=AMAX1 (ABS(XMAX) ,ABS(XMTN) ) 

0  0  53 

IF(XM.GT .9999. )GO  10  11 

0  054 

13 

I F ( XM . GE .1000.) GO  10  14 

0055 

NMIX=NMIX+1 

0  OS6 

XMX=XMX>K  1  0  . 
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0  057 

XMN=XMN*i 0 . 

0  058 

XM=XM*10 , 

0  059 

GO  TO  13 

0  060 

1 1 

NPLX=NPL..X+i 

0  0  6 1 

XMX=XMX/10 . 

0062 

XMN=XMN/10 , 

0  0  63 

XfUXM/i  0  . 

0  064 

I F ( X M , G T , 9 9 9 9 , ) GO  TO  11 

0  0  6  5 

14 

MXX=IF IX(XMX) 

0  0  66 

MNX=IFIX(XMN) 

0  0  67 

DGX= ( XMX-XMN ) /I 0  0 , 

0  0  68 

NMIY-0 

0  069 

NP LY  =  0 

0  07  0 

YMX=YMAX 

0  0  7  i 

YMN-YMIN 

0  072. 

YM=AMAXi ( ABS ( YMAX ) , ABS ( YMIN ) ) 

0  073 

IF ( Y  M  ,  G  T  .9999 . )GO  10  1 

0  074 

3 

IFtYM.GE. 1000 . )GO  ID  4 

0  0  75 

NMIY=NMIY+1 

0  0  76 

Y Pi X = 1  0  .  *YMX 

0  077 

YMN=YMN*i 0 . 

0  0  78 

YM=YM*10 . 

0  079 

GO  TO  3 

0  0  80 

1 

NPLY=NPL.Y+i 

0  0  8 1 

YMX=YMX/i 0  . 

0  082 

YMN  =  YMN/1 0  . 

0  083 

YfUYM/10  . 

0  034 

IF ( YM . G  T . 9999 .  ) GO  TO  1 

0  0  85 

4 

MXY=IP IX(YMX) 

0  0  8  6 

MNY=IFIX(YMN> 

0  087 

DGY = ( YMX-YMN ) /5 0 , 

0  0  88 

UP 

FORMAT (A) ) 

0  089 

NT  =  i  0 

0  0  90 

IT AB=440  OB 

0  0  9 1 

IF ( NPL  Y  .  GT  . 0>WRITE( 1 , 124 ) I TAB , ITAB , NPLY 

0  0  92 

124 

FORMAT ( 2A1 , "#E#*+ " , 12/ ) 

0  0  93 

IFCNMIY , GT . 0 ) UR ITF ( 1 , 125 > ITAB , ITAB , NMIY 

0  0  94 

125 

FORMAT (2A1, "♦E**-" , 12/ ) 

0  0  95 

I RN-6400B 

0  096 

DY=(YMAX-YhIN)/50 . 

0  0  97 

DO  20  1=1,50 

0  0  98 

WRITE ( 1 ,104 ) ITAB 

0  0  99 

LL  =  MOD ( I  ,2) 

0 .1  0  0 

URITEd  ,  101)  15446B,2Hk0 

0  i  0  i 

UR  TTE ( 1  ,  104) 1HS 

0  10  2 

IFtLL.EQ  .  0)  URITEd  ,  1  04  )  LABEL  <  1/2  ) 

0103 

IF  (  LL  .  NE. .  0  )  GO  TO  17 

010  4 

I.F  (LABEL.  (1/2)  ,  EQ  .  0  )  UR  ITE  (  1  , 1 13) 

0  1  0  S 

17 

IF  (L.L..  .  NE  .  0  )  WRITE  (  1  ,  .1. 13) 

0  1  0  6 

113 

FORMAT (”  ") 

0107 

URITEd  ,  10 1 )  15446B , 2Hk2 

0108 

URITEd  ,104)  1 H8 

0109 

104 

FORMAT ( A1 , 

0110 

URITEd  ,104)  ITAB 

0111 

I G=I P IX ( MXY- (I — 1 >*DGY) 

Oil? 

I GG=MOD d-1 ,10) 

,  1971) 
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0113 

IF  <  IGG  .  EQ  . 0 ) WRIT  E( 1 ,105)  IG  ,ITAF< 

0114 

105 

FORMAT < 15, Al,"  " ) 

0115 

IF ( IGG , NE , 0 )GO  TO  21 

0116 

WRITE* 1  ,  116) 

0117 

116 

FORMAT ( " +  ") 

0  118 

DO  10  J=i , NT 

0 1 1 9 

WRITE ( 1 ,103) 

0120 

1  03 

FORMAT  (" - +  •') 

0121 

10 

CONTINUE 

0122 

WRIT  E( 1  ,  1 0 4  )  I R N 

0123 

WRITE  C  1 , 123) ITAB, ITAB , ITAB 

0124 

123 

FORMAT ( 3A1 , "  ") 

0125 

21 

IF ( IGG . EQ . 0 ) GO  ID  5 

0126 

WRIT  F.  <  1  ,  106)  I TAB 

0127 

106 

FORMAT ( Al ," :  " ) 

0128 

WRITEC 1  ,  104) IRN 

0129 

WR IT  E ( 1  , 1 23 ) 11 AB , I T AB , IT AP 

0130 

5 

DO  70  K=i , 1 00 

0 1 3 1 

II...  (K)=1H 

0132 

70 

CONTINUE 

0133 

Yi=YMAX~(  I - .1  )  *DY 

0134 

Y2=YMAX-I*DY 

0135 

DO  9  J“ 1,100 

0136 

GO  TO (7, 8) IT 

0137 

7 

DO  60  K=NTR , 1 , -i 

0138 

IF (Y( J,K) .LE.Yi . AND . Y ( J , K ) . GT . Y2) 1L( J )=IS(K> 

0139 

IF  ( Y  <  J  ,  K  )  .  EQ  .  Y2 )  IL(  J  )  =  1.8  ( K  ) 

014  0 

60 

CONTINUE 

0141 

GO  TO  9 

0  142 

8 

I F  (  Y  (  J  )  .  GT  .  Y 1 )  1 L  ( J  )  =  I S  ( 5  ) 

0143 

9 

CONTINUE 

0144 

IF  ( IGG  ,  EQ  .  0  )  WR  I TE  ( 1  , 1 07  )  IL. 

0.1.45 

IF ( IGG  .  NE . 0 ) WR IT E  < 1 , 1 1  4 ) IL 

0146 

107 

FORMAT  (  100A1 ) 

0147 

1  14 

FORMAT (100A1 ,  "  :  “  ) 

0148 

IF (  IT , EQ . 1 ) GO  TO  2  0 

0149 

DO  19  J  =  1  , 1  0  0  0 

0150 

DO  19  K=i,75 

0151 

CONTINUE 

0152 

19 

CONTINUE 

0153 

20 

CONTINUE 

0154 

WRITE  (1,115)11  AB, 1TAB,MNY  ,  ITAB 

0155 

115 

FORMAT  (2A1  ,15, Al,"  +  ...") 

0156 

DO  30  1=1, NT 

0157 

WRITE*!, 103) 

0158 

30 

CONTINUE 

0159 

WRITE ( 1 , 112) IRN 

0160 

WRITE* 1  , 112) IRN 

0161 

WRITE* 1  ,  108) ITAB 

0162 

1  08 

FOR  MAT  (  A 1,"  ...") 

0163 

DO  40  1=1 , NT+i 

0164 

IG=IFIX(XMN+(I-1 )*DGX*10  .  ) 

0165 

WRITE* 1  , 109) IG 

0166 

109 

FORMAT <5X, 15, "  " ) 

0167 

4  0 

CONTINUE 

0168 

IF (NPLX.GT.O ) WRITE <1,110 )NPLX 

,  1971) 
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01 69 

017  0 
0171. 
0172 
0173 
0174 
0175 
0176 

0  1 77 

J  F ( NMIX . GT  .  0)WRITF.(1  ,111) NMIX 

110  FORMAT  ("  "  ,  12) 

111  FOR  MAT <"  *E**-",I2) 

URITEd  >101)154466,  2Hk  0 

UR  J  TE  (1,112)  .1.  HEi 

URITEd  ,101)154466, 2H16 

UR  I TE ( 1 ,112) 1 HD 

RETURN 

END 

>K  X  NO 

ERRORS* 

71361 
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APPENDIX  A- 5 
PROGRAM  FILTR  (FILTeR) 

Descri ption 

Program  FILTR  applies  an  eighth  order  bandpass  filter  to  from  one  to 
eight  channels  of  data  stored  on  disc.  The  eighth  order  filter  is 
obtained  by  cascading  two  identical  fourth  order  filters  which  have 
been  obtained  from  a  bilinear  transformation  of  a  second  order  Butter- 
worth  equation.  The  difference  equation  for  the  filter  is  obtained 
from  a  subroutine  BPASS  for  which  the  listing  is  given.  Figure  A- 5-1 
is  a  flow  chart  of  program  FILTR.  The  program  is  given  in  Program 
Listing  A- 5-1 . 

Required  input  parameters: 

1.  Upper  and  lower  frequencies. 

2.  All  other  required  information  such  as  sampling  rate,  number  of 
data  channels,  and  record  length  is  obtained  from  the  data  header 
record  on  disc. 

Required  subroutines: 

1 .  DREA 

2.  DWRI 

3.  STPT 

4.  BPASS  -  see  Program  Listing  A-5-2. 
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Figure  A-5-1  Flow  chart  of  program  FILTR 
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0  0  0  1 
0  0  02 
00  03 
0  004 
0  0  05 
0  0  06 
0  007 
8008 
0  0  09 
0010 
0011 
0012 
0013 
00  14 
00  15 
0016 
0017 
0018 
00  19 
0020 
0021 
0022 

0023 
0  0  24 
0025 
0026 
0  027 
0  028 
0029 
0  030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0  038 
0  039 
0040 
0  041 
0042 
0  043 
0  0  44 
0045 
0  046 
0047 
0  048 
0  049 
0050 


FTN4 , L 

PROGRAM  FILTR 

DIMENSION  L.ABF  (135)  ,A(5>  ,B<5)  ,Y(256)  ,  IX (256)  ,  1X1(8),  Yl(8) 
EQUIVALENCE  (1X1(5), IX), (Yi(5),Y> 

CALL  EXEC( 17, IFTRK , ILTRK , ISIZE) 

CALL  DREA( IFTRK ,0 ,LABF ,1 , 1 , 1 ) 

SR  =  LABF ( 38 ) /60 , 

NR  ECS  =  LABF ( 39 ) 

NCHA=l.ABF  ( 37 ) 

WRITE  <  1  , 1 00 ) 

100  FORMAT ( /"INPUT  LOWER  AND  UPPER  FREQUENCIES  IN  CYCLES/MIN. 
READ ( 1  ,*  )LF , IMF 
ALF=LF/60 . 

AHF=IHF/60 . 

CALL  BP  ASS  ( SR  ,  AL.F  ,  AHF  ,  A  ,  EO 
DO  5  K  =  1  ,  4 
1X1  <  K ) =  0 
Y1 ( K ) =0 
5  CONTINUE 
DO  10  L.  =  l  ,2 
DO  10  1  =  1,  NCHA 
DO  10  J= 1 , NRECS 
M  AX=  0 

CALL  DREAdFTRK  ,256,  IX,  J  ,  I  ,1) 

DO  11  K  =  1  ,4 

Y(K)=A(1)*IX(K )+A(2)*IXi (3+K  )+A(3)*IXi (2  +  K )  + 

1  A(4  >  *1X1 ( 1+K )+A(5)*IXi (K  )+B(i )*Y1 (3  +  K  )  + 

2  B( 2 ) *Yi (2  +  K )+B(3)*Yl ( 1 +  K  )  +  B ( 4 ) *Yi ( K ) 

11  CONTINUE 

DO  12  K  =5 , 256 

Y(K)=A(i)*IX(K)+A(2)#IX<K-i>+A(3)#IX(K-2)+ 

1  A(4)#IX(K-3)+A(5)*IX(K-4)+B(l)*Y(K-l)+ 

2  B(2)*Y(K-2)+B(3)*Y(K-3)+B(4)*Y(K-4) 

12  CONTINUE 

DO  13  K=i ,4 
I X 1 ( K  )  =  I X ( 252  +  K ) 

Y 1  ( K ) =Y ( 252  +  K  ) 

13  CONTINUE 

DO  14  K  =  1 ,256 
IX ( K ) =Y ( K  > 

IM  =  I ABS ( I  X ( K  )  ) 

IF ( IM.Lt . MAX ) GO  10  14 
MAX=IM 

14  CONTINUE 

CALL  DURI ( IFTRK  ,256 , IX , J , I , 1 ) 

I F ( L , EQ . 1 ) GO  TO  10 

CALL  STPT( IFIRK , IX, 256, J , I ,MAX> 

10  CONTINUE 
STOP 
END 


**  NO  ERRORS* 


Program  Listing  A-5-2  BPASS 


PAGE  0  0  0  1  (FTN4  —  RELEASE.  24 177B-- JULY  , 


0001  FTN4.L 

0002  SUBROUTINE  BP ASS ( SR , F 1 , F2 , A , B ) 

0003  DIMENSION  A(S),B<5) 

0004  T=1./SR 

000S  Ui=6.2B32*Fi 

0006  U2=6.2832*F2 

0007  D=COS ( (U2-W1 )*T/2 . )/SIN( (U2-U1 )*T/2 . ) 

0  0  08  D1=SQRT (TAN(Ul*T/2  .  ) *T AN ( W2*T /2 ) ) 

0009  D2  =  2 . * AT  AN ( D 1 ) 

0010  E=2*C0S  <  D2 ) 

0011  Z0=D**2+i . 41 42136*D+i , 

0012  Zi=-2  .  *E*D**2- 1  . 4142136*E*D 

0013  Z2=(E**2+2 . )*D**2-2 . 

00  14  Z3  =  -2  .  *E*D**2+1  .  4142136*E*D 

0015  Z4=D**2-1 . 4142136*0+1 . 

0016  A  <  1  >  =  1  ./ZO 

0017  A  (  2 )  —  0  . 

0018  A ( 3  )  =-2  .  /Z  0 

0019  A ( 4 )  =  0 . 

0020  A(5)=A(1) 

0021  B(i)=-Zi/Z0 

0022  B ( 2 ) =-Z2/Z  0 

0023  B ( 3 ) =-Z3/Z0 

0024  B(4)=-Z4/Z0 

0025  RETURN 

0026  END 


**  NO  ERRORS* 


1971) 
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PROGRAM  DCAC 
Descri ption 

Program  DCAC  determines  the  average  or  dc  value  of  each  channel 
of  a  data  file  on  disc,  and  then  removes  this  value  by  subtraction. 
The  program  is  given  in  Program  Listing  A-5-3. 

Required  input  parameters:  None. 


Subrout 

ines  required: 

1 . 

DREA 

2. 

DWRI 

3. 

STPT 

0  0  0  i 
Q  0  i)2 
0  0  03 
n  o  0  4 
0  0  OS 
0  006 
0  0  07 
0  0  0  6 
0  0  09 
0010 
00  i  1 
00  12 
0  013 
0014 
00  IS 
0  0  16 
00  17 
0  018 
0  0  1  9 
0020 
0  021 
0  0  22 
0023 
0024 
0  0  2S 
0  0  26 
0  077 
0  028 
0029 
0  0  3  (i 

00  31 
0  0  32 
0  033 
0  0  34 
0  0  3S 
0036 


Program  Listing  A-5-3  DCAC 
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ETN4 , 1 

PRUU  Ah  DC AG 

DIMENSION  l  ABF ( 133)  ,Y<256)  , IX (256) ,A(B) 
CALI  EXECt  17,  IF  IKK  ,  ILTRK  ,  J SIZE) 

CAI  I  DR  F  A(  IFTRK  ,0,1  ABF  ,1,1,1) 

SR  -  LABI  (  38  )  /()  0  . 

NR  EEC  l ABF (39) 

NCMA-L ABF i37) 

DO  12  1=1 , NCHA 

Dll  12  J  =  1  ,  NR  ECS 

CALL  DR E A ( IFTRK ,?B6, IX , J , 1,1) 

DO  13.  K  =  1 , 2S6 
A(I)=A(I)-»IX(K) 

12  CONTINUE 

DO  S  K  =1  ,  NCUA 
A  (  K  )---A(X  )/<NRECS*256  ) 

S  CONT I NUt 

DO  10  1=1, NCUA 
DO  10  J= 1 , NRFCS 
MAX  =  0 

CAI  L  DREA( IFTRK , 2 5 fa, IX, J, 1,1) 

DO  1 1  K  =  i  ,256 
Y(K )  =  IX(K  )-A( I) 

11  CONTINUE 

DO  14  K  —  1 , 256 
IX(K)=Y(K) 

IM-IABS( I  X ( K )  ) 

IF  < IM  .  L  T  .  MAX ) GO  TO  14 

MAX - TM 
14  CONT  INI  IE 

CALL  I)WR  1  i  I  FI  RK  ,256 , 1  X  ,  J  ,  I  ,  1  ) 

I F ( L  .  F  Q .  1  >  GO  TO  10 

CAL  L  STP1 ( IFTRK , IX ,236, J , 1 ,MAX) 

10  CONTINUE 
STOP 
END 


**  NO  ERRORS* 
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APPENDIX  A-6 

PROGRAM  PSPT  [Power  SPecTrum) 

Descri ption 

Program  PSPT  produces  single  or  multiple  power  spectra  for  up  to 
8  channels  of  data  stored  on  disc.  The  complete  data  file  may  be 
divided  into  n  segments  with  power  spectra  produced  for  each  segment. 

The  power  spectra  are  obtained  from  the  application  of  the  fast  fourier 
transform  to  each  256  word  record  in  a  data  segment.  The  power  spectrum 
for  each  segment  is  an  average  of  the  power  spectra  obtained  from  the 
records  in  the  segment.  The  magnitudes  of  the  power  spectra  are  output 
on  the  teleprinter  using  subroutine  PLOT.  The  flow  chart  of  program 
PSPT  is  given  in  Figure  A-6-1  and  a  typical  output  is  shown  in  Figure 
A-6-2.  The  program  is  given  in  Program  Listing  A-6-1. 

Required  input  data: 

1 .  Number  of  channel s . 

2.  Whether  or  not  absolute  magnitudes  are  required  as  an  output. 

3.  The  number  of  segments  into  which  the  data  record  is  to  be  divided. 

4.  The  number  of  times  the  output  data  is  to  be  smoothed  using  a 
moving  average  filter. 

5.  The  number  of  plot  points  to  be  surpressed  at  the  low  frequency 
end  of  the  spectrum. 

Subroutines  required: 

1.  IFFT  -  the  fast  fourier  transform  subroutine  written  by 
Dr.  Z.  Kol es ,  U .  of  A. 


2.  DREA 


.  ' 
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3.  PLOT 

4.  SMTH  -  a  moving  average  filter.  See  Program  Listing  A-6-2. 
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Figure  A-6-1  Flow  chart  of  program  PSPT 
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Figure  A-6-1  Flow  chart  of  program  PSPT  (cont.) 
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Figure  A-6-2  Typical  single  channel  output  from  program  PSPT 
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Program  Listing  A-6-1  PSPT 


PAGE  0001 


(FTN4— RELEASE  24177B  — JULY,  1971) 


FTN4  ,  L 

PROGRAM  PSPT 

DIMENSION  LABE (135)  ,ICH(8) ,Y(128,8) ,ISINES<65) ,IP<4)  , 

1  JR  AY (256) , IR ( 256 )  ,IM(256)  ,NFL0U(8) , YM ( 8 ) , X ( 128 ) ,IESC(9)  , 

2  LABEL.  ( 50  ) 

DATA  LABEL/18#1H  , 1HF , 1HF , 1HT  ,  1H  , 1HM  ,  1HA , 1HG , 1 HN , 1HI , 1HT , 1HU , 
1  1HD.1HE/ 

DATA  IESC/15446B ,2Hk2 ,2HS  , 2H18 , 2H1 1 , 2HV  ,2HD  ,2Hk0,2H16/ 

CALL  EXEC( 17,IF1RK ,1LTRK , ISIZE) 

CALL  DREAdFTRK  ,0,LABF,1  ,1,1) 

IWID1 H=  1 0  0 
1  WRIT  L( 1 , 100 ) 

100  FORMAT (/"INPU1  1  HE  REQUIRED  CHANNELS?  _M) 

READ ( 1  ,*  ) ( ICH( 1 ) , 1  =  1 ,8) 

DO  10  1=1,8 

I F ( 1 CH ( I ) . EQ . 0 ) GO  TO  11 

10  CONTINUE 

11  NCHA= I - 1 
NRECS=LADF ( 39 ) 

SR=LABF ( 38 ) /6  0 . 

13  UR ITE ( 1 ,101) 

101  FORMA! (/"DO  YOU  WAN!  AbSOLUTE  MAGNI1 UDES?  _"> 

READ ( 1  , 1 02) ISC 

102  FORMA! ( A2) 

IF ( ISC . EQ . 2HNO)GO  TO  12 
IF ( I  SC . NE . 2HYL ) GO  TO  13 

12  WRITE(1 ,104) 

104  FORMAT (/"NO  OF  SLGMEN1 S?  _"> 

READ ( 1  ,* )NSEG 

IF ( NSEG . LE . 0 ) GO  TO  12 
IF(NSEG.GT,  NRECS  )NSEG==NRECS 
ISEG=NRECS/NSLG 

14  URIlEd  ,105) 

105  FORMAT(/“SMQO!H  N  11MES?  _"> 

READ (  1  ,*)NFILT 
IF(NFILT.LT. 0)GO  IO  14 

15  URI1 E ( 1  , 113) 

113  FORMA! (/"SUPPRESS  PLOT  POINTS?  _") 

READ ( 1  ,* )NSUP 
I F ( NSUP . L  T . 0 ) GO  TO  15 
NPLOT  =  128 

SEGL= ( 60 . /LABF (38 ) )*2S6*ISEG 

LREC=SEGL*NSEG 

LSEG  =  SEGL. 

MSEG=SEGL/60 . 

MREC  =  L.REC/60  . 

ISECR=MQD(LREC  ,  60  ) 

ISECS=MOD(LSEG ,60 ) 

URI1E(1  , 107)MREC , ISECR 

107  FORMAT (/"RECORD  LENGTH*" ,X , 13 ,X , "MINUTES" ,X , 13 ,X , "SECONDS" ) 
URIlEd  ,  108) MSEC,  ISECS 

108  FORMAT (/"SEGMENT  LENGTH* " , X , 13 ,X , "MINUTES" ,X , 13 ,X , "SECONDS" ) 

DO  29  M= 1 , NSEG 

IF ( NCHA . GT . 4 ) GO  TO  17 
I TIMES=NCHA 
GO  10  18 
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Program  Listing  A-6-1  PSPT  (cont.) 


P  ALE  0UQ2 

♦ 


PSPT 


(FTN4-- RELEASE  24 177B — JULY ,  i97i) 


* 

i 


0057 

17 

ITIHE.H-4 

0058 

18 

N=i 

0059 

I SND=i 

0060 

I AGAIN=1 

0061 

20 

DO  23  1  =  1  ,1) 

0062 

YM ( I ) =0 . 

0063 

DO  23  K  =  i  ,128 

0  064 

Y ( K , I ) =0  . 

0  06S 

23 

CONTINUE 

0066 

YMAX  =  0  . 

0067 

NT  =  0 

0068 

DO  35  I =N , I T I MES 

0069 

DO  22  J=i , ISEG 

0070 

NFLOW ( I )  =  0 

0071 

I  SET  =  1 

0072 

CALL  DREAdPTKK  ,256,JHAY,  J  ,ICH(I)  ,1) 

0073 

34 

DO  19  K  =  1  , 256 

0074 

IM  (  K  )  =  0 

0075 

IR(K)=IFIX(JRAY<K)*(i.+SIN(<K-i)*6. 2832/255.-1 .5708) )/( 16 

0076 

19 

CONI  1 NUE 

0077 

NT=N1 +1 

0078 

CALL  IFF! ( IR , 1 M  ,  I5INES ,256  ,NT  ,0 , IFL) 

0079 

IF  (  IF  l. .  EQ  .  0  )GO  10  21 

0080 

NFLOU  <  I  )=NFLOU( I )  +  i 

0081 

ISET=ISFT*2 

0082 

GO  TO  34 

0083 

21 

SET=ISET**2 

0084 

DO  22  K  =  1  , 128 

0085 

Y ( K , 1 ) =Y ( K , I )  +SET  # ( FLOA1 ( I R ( K ) ) **2  +  FL0AT ( IM ( K ) )**2) 

0086 

OO 

t— 

CONI  1 NUE 

0087 

Y(  1 , 1 >  =  0  . 

0088 

IF(NFILT.LE.0)GO  TO  24 

0089 

DO  25  K  =  1  ,NFIL1 

0090 

25 

CALL  5MTH ( Y ( 2 , 1 )  ,  127) 

0091 

24 

IF ( NSUP . LT . 1 ) GO  TO  26 

0  092 

DO  27  K  =  i  ,  NSUP 

0093 

27 

Y(K  , I >  =0 . 

0094 

26 

DO  35  K  =  1  , 128 

0  095 

Y ( K ,1 ) = Y ( K ,I)/ISEG**2 

0  096 

IF ( Y (K , I ) .GE . YM( I ) )YM( I )=Y(K , I ) 

0097 

35 

CONI 1NUE 

0098 

DO  28  I=N , ITIMES 

1099 

I F ( YM ( I ) . GE , YMAX ) YMAX=YM ( I ) 

0100 

28 

CONI  1 NUE 

0101 

IF  (  1SC.EQ.2HYDGQ  10  30 

0102 

YMAX  =  1 . 

0103 

DO  31  Id  ,8 

0104 

IF  (  YM  ( I )  .  EQ.  0  )  YM(  I  )  =  i  .E36 

0105 

DO  32  K  =  1  , 128 

0106 

Y  <  K  ,  1  )  =Y  ( K  ,  I  )  / YM ( I  ) 

0107 

IF  <  Y  (  K  ,1)  .  GT  .  1 )  Y  (K  ,  I )  =  1  • 

0108 

32 

CONI  1 NUE 

0109 

31 

CONI  1 NUE 

0110 

30 

S  =  FLC)AT  ( LABF  (38)  )/256. 

0111 

DO  33  1=1,128 

0112 

X( I )=( I-i )*S 
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Program  Listing  A-6-1 


PSPT  (cont.) 


0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 


PAGE  0003  PSPT  ( FTN4  — RELEASE  24 177B-- JULY ,  1971) 


33  CONI INUE 
NTR  =4 

I F ( NCHA . LE . 4 ) NT  R  =NCHA 
60  IT=i 

IF (NTR . EQ .  1  ) IT  =  2 

CALL  PLOT(X  ,  Y(512*(  ISND-i  )+l  >  , 1 . 5 , LABEL , 128 , NTR  ,11  ) 
URIlt(l  ,109)H, 1 SND  ,  ( LAbF ( I ) , 1  =  1 ,36) 

109  FORMAT  ( 18X  ,  "SEGMENT  ",  12 , 11  PART  "  ,  1 1  , 1  OX  ,  36A2  ) 

IF  (NIK  .EQ  .  1  >GCJ  10  40 
URI1E(1,110)(1CH(I),I=1,4),21000B 

110  FORMA! <30X,"CH. ",12,“  x " , 8X , "CH . " , 12 , "  o  "  ,  8X  ,  "  CH . " , I 
1  8X , "CH . " , 12, "  " , A1 ) 

URI1E(  1  ,111) (NPLOU(I) , I =N , I T I MES ) 

111  F  OR  MAT ( 1 4X , "  OVERFLOWS:  ",4(I3,12X)> 

WRIT E(1 ,112)(YM(I) , I=N , ITIMES) 

112  FORMAT ( 1 4X , “  SCALE  FACTORS:  "  ,  4 ( El  0 , 5 , 5X ) ) 

4  0  IF(N1R .EQ. 1 )WR1TE(1 ,151 )ICH( 1) ,NFLOW( 1 ) ,YM( 1) 

151  F0RMAT(28X, "CHANNEL  NO.  ",  1 1 , 4X ,  “OVERFLOWS :  "13, 4X, 

1  "SCALE  FACTOR:  ",E10.5> 

♦  N  =  5  ' 

III MES-NCHA 

1AGAIN=-1 AGAIN 
ISND=ISND  H 

IF( 1AGAIN.LT . 0 .AND. NCHA. GT.4)G0  TO  20 
29  CONTINUE 

WRI1E( 1 , 150 ) IESC( 1  )  , IESC ( 5 )  ,IESC(6) 

150  FOR  MAT ( 3A2 ) 

STOP 

END 


**  NO  ERRORS* 
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Program  Listing  A-6-2  SMTH 


PAOE  0  0  0  1  (  FTN4  —  RELEASE  24 177B — JULY  ,  1971) 


SUBROUTINE  SMTH(X,NPTS) 
DIMENSION  X ( 3 ) 

A=X ( 1  ) 

B=X ( 2 ) 

C  =  X  (  3  ) 

X( 1  )= .5* ( A  +  B) 

DO  20  J=2,NPTS-1 
X< J)=(A+2.*B+C>*.25 
A=B 
B=C 

20  C  =  X  ( .1  +  2  ) 

X ( NPTS ) = . 5* ( A+B ) 

RETURN 

END 
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APPENDIX  A-7 
PROGRAM  PHSE  (PHaSE) 

Descri ption 

Program  PHSE  uses  zero-crossi ng  techniques  to  determine  cycle-to- 
cycle  frequencies  and  relative  phase  shifts  of  any  two  data  channels 
stored  on  disc.  The  frequencies  and  relative  phase  shift  are  plotted 
as  a  function  of  cycle  number  using  the  subroutine  PLOT.  Linear  regres¬ 
sion  analysis  is  applied  to  the  arrays  of  phase  angles  and  frequencies 
to  produce  equations  describing  changes  in  these  parameters  as  a  function 
of  cycle  number.  A  flow  chart  of  program  PHSE  is  shown  in  Figure 
A- 7-1  and  a  typical  output  is  shown  in  Figure  A-7-2.  The  program  is  given 
in  Program  Listing  A- 7-1 . 

Required  input  data: 

1.  The  numbers  of  the  two  channels  to  be  processed. 

2.  All  remaining  data  is  obtained  from  the  data  header  record  on  disc. 

Subroutines  required: 

1 .  DREA 

2.  PLOT 

3.  LFIT  (Unear  FTT)  -  a  subroutine  to  do  linear  regression 

analysis.  See  Program  Listing  A-7-2. 


. 


202 


Figure  A-7-1  Flow  chart  of  program  PHSE. 
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Figure  A- 7-1 


Flow  chart  of  program  PHSE  (cont.)- 
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Figure  A-7-2  Typical  output  from  the  program  PHSE 
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Program  Listing  A-7-1  PHSE 


P  AC.i_  OOOi  (  FTN4 — RELEASE  24177B  —  JULY,  1971) 


FTN4 , L 

PROGRAM  PIISE 

DIMENSION  1 F  SC ( 9  )  ,LABF(135) ,JRAY(255>  ,  KRAY (256)  ,Y(100>3)>Fi(10n) 
1  F  2 ( 1 0  fl ) ,  P  H  ( 1 0  0  >  ,  LABEL (25)  ,  X ( 1 0 0  )  , CPI  I  ( 1 0 0 ) 

DOUBLE  PRECISION  AF1 , SD1  ,  AF2  , BD2 , AP H , SD3 , ACPH , SD4 

DATA  I  ESC/ 15446B  ,  2Hk2 , 2HS  , 2H18 , 2H1 1 , 2HU  ,2HD  , 2Hk 0 , 2H16/ 

CALL  EXFC(17,IF1RK,IL,IS) 

CALL  DREA( IFTRK , 0  ,LABF, 1 , i , 1 ) 

SR  =  1  ABF ( 36 ) /60 . 

DT  - 1  .  /SR 
NRFCS=l  ABF  (  39  ) 

wrii r ( i ,  ioo ) 

100  FORMAT ( /“INPUT  1  HE  REQUIRED  CHANNELS!  .") 

RF  AD ( 1  , * ) I CH 1 , 1CH? 

1=0. 

111  =  0 
T  1  =  0  . 

T  2  ~  0  . 

1P  =  0  . 

J=1 
K  =  1 
L  =  1 
IP  =  0 

DO  1 O  M-  1 , NK ECS 

CAL  L  DRE A ( IF  TRK , 256 , JR AY , M , ICH1  ,1) 

CALL  DREA( IF  TRK , 256 , KRAY , M , 1CH2 , 1 ) 

11  =  1 

*  IF  (H  .  NE  .  1  >GO  10  11 
DO  12  1=2,256 
I SN1  =  1  S  1 GN ( 1 , J  R  A  Y ( T - 1 ) ) 

TSN2=1  S  l  f.N  <  1  ,  JR  A  Y  ( I  )  ) 

IF ( I SN2  .  NE . ISN1  .AND.JRAY(I) . GT . 0 ) GO  10  13 
1?  CONTINUE 

13  1  1  =  (  F  LOAT  (  JR  AY  (  l  )  )/  (FLOAT  (  J  R  AY  (  I  )  ) -FLOAT  (  JR  AY  (  i-1 )  )  )  )*l)1 
I  1  =  1 

11  DO  2  0  I  -  II  ,256 

IE  (  I  .  NF  1  )GU  TO  14 
I SN i  =  1 5 1 GN ( i , I C 1  ) 
f,n  TO  15 

14  ISNl=lSIGN(i , JR AY ( I ~i ) ) 

15  ISN2=TSIGN( 1 , JRAY(I) ) 

IE (1SN2.NE.ISN1 . AND , JR AY ( J ) .G1 .0)00  10  16 
GO  TO  17 

16  IF  (  1  .  f  ■  T  .1  )  It!  1  -.1  RAY  (  I  - 1 ) 

IF  (  11 1  . EQ . 0 )GU  10  IB 

1  =T-D1 * (FLOAT ( JR  AY ( 1 ) ) /(ELOAT ( JR AY ( I ) ) -FLOAT ( IC1 ) ) ) 

F 1 ( K  )  =  5  0  .  /(T-1 1  ) 

IF(M  (K  )  .GT  .30  .  ) GO  10  70 
AL'F  1  =F  1  (K  ) 

IF  <  K . OT  .  3) AOF 1= ( E  1  ( K-l ) +  E 1 <  K  -  2 )  +  F 1 (K-3) )/3 , 

I F  (  E  1  <  K  )  .  GE  .  ( 1 . S#AUE‘  1 )  )  CO  10  70 
IE  =1 

IF <F 1 (K  )  .LT . ( . 67*A0Fi ) )GO  10  18  » 

K  -K  + 1 

IE  (K  .01  .  100 ) GO  10  25 
18  T 1 =T 
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Program  Listing  A- 7-1  PHSE  (cont.) 


F'  AGE 

ii  0  0  2  PHSE  (FTN4 — RELEASE  24177B— JULY,  1971) 

0057 

7  0 

T=T  ♦  DT*( FLOAT (  JkAYd  )  )  /  (  FLOAT  (  JR  AY  (  I  )  ) -FLOAT  (ICi)  )  ) 

0  056 

I T1  =  1 

0  0  59 

17 

1F(  I  .  NF;  .  1  )GO  TO  19 

0  060 

ISN1=1SIGN( 1 , IC2 ) 

0  061 

uo  to  i 

0062 

19 

ISNl=IsiGN(i , K  R  A  Y  ( I - 1 ) ) 

0  0  63 

1 

I SN2= 1 S I GN ( 1 , KR AY ( I ) ) 

0064 

IF (  ISN? . NE . ISN1 . AND , KRAY ( I ) . GT  .  0  ) GO  TO  2 

0065 

GO  TO  3 

0  0  66 

o 

IF  (  I  .(.1  .  1  )I(;2-KK’AY<  I-i  ) 

0067 

T  -  1-1)1  *  (FI  OAT  (KRAY  ( I )  )/<M  OAT  (KRAY  (  1 )  ) -FLOAT  ( IC2)  )  ) 

U  068 

IF (  r? . fc  Q . 0 )GO  10  4 

0  0t.9 

F-  8  ( L  )  =  6 0  .  /  ( T-l  8  ) 

0070 

1F(F21L).GT.30.  ) 1  =  T  +  DT* ( F LOAT ( KR AY ( I ) ) / ( F LOAT ( K R AY ( I ) >- 

0  0  71 

1 

F  LOAT ( IC?) ) ) 

0  0  72 

IF  (  F  2  (  L  )  .GT..50.  )GO  10  3 

0  0  73 

AUF  2  =  F  8  (  L  ) 

f  074 

IF(I  .Gl  .  3  )  AUF2  =  ( K?  ( L-l )  +F2  ( L  -2 )  +  F2  ( 1.-3  > ) /3  . 

0075 

IF  (F  2(1  )  .  GE  .  (  1 . 5*AUF2>  >  T  =T  +  TD*  (  FL.OA  I  (KRAY  ( I  )  )  /  (  FLGA1  (KRAY  ( I )  ) 

0  076 

1 

-F  L  OAT ( 1C2  )  )  ) 

0  0  77 

IF ( F 2 ( 1  )  .GE.  (1  . 8*AUF2) )LO  TO  3 

J  0  78 

I F  ( F  2  ( 1  )  .LT  .  (  .  67*AUF2)  )GO  TO  4 

0  0  79 

l  -1  +1 

0  UffO 

IF<  1.  .  l,T  .  1  0  0  >GO  10  25 

nom 

4 

T2-T 

0  0  8  3 

T-T4D1  *  U  l  OAT  (KRAY  (  1  )  )/(  F  l.OAT  ( KRAY  (  1  )  > -M  OAT  (  TCP) )  ) 

0083 

TP  T2-T1 

0  0  84 

IF  (  IP  .  F  U  .  II  )GO  TO  3 

0  085 

PH(J )-<TH*Fl (K-l >/A0 . ) *86 0 . 

0  0  86 

IP^O 

0  087 

J-.J  +  l 

0  0  88 

IF (J . GT .100)00  T O  25 

0  0  89 

J 

i  =  r  +  d  r 

0  090 

20 

CONTINUC 

0  091 

IC) - JR A Y ( 256) 

0  0  92 

IC2=KRAY (256) 

1)0  93 

1  0 

CONT  I NUE 

0  094 

■jc 
»-  _> 

IF  (FI  (K  )  .GT  .  3(1  .  )P  1  (  K  )  =30  . 

0095 

IF(F  2(1  )  .  GT  .  3(1  .  )F2(L  )  =  30  . 

0096 

AF1 =0 . 

0  097 

SD 1  -  0 

0  098 

l  =1.-1 

0  099 

K-K  -1 

0100 

J  -  J  - 1 

0101 

DO  8  0  I  1  ,K 

0102 

AF  1-AF  1+F  1  (  I  ) 

0103 

SD1  -SD1  t  FI  ( i  >**;/ 

U  1  04 

3  0 

CONI  1 NUE 

0105 

CALI  1  F I T ( F 1  , K , A01 ,Ail  ) 

Olio 

AF  1  =AF  1  /K 

0  1 C  7 

SD1-5M /K 

0  1  08 

SI)1=DAI:S(SD1-AF1**2) 

0109 

5D1  -I'SQRT  (  SD  1  ) 

0110 

AF  2-0 . 

0111 

SD2=  0 . 

0112 

DO  4  0  l  1,1. 

* 
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Program  Listing  A-7-1 


PHSE  (cont.) 


PACE  00  Uw  PHSL  (FTN4 —  RELEASE  24 1 77B  —  JULY  ,  1971) 


(■  113 
t)  1  14 
01  IS 
0  116 
0117 

u  i  ia 

Ui  19 
0120 
0121 
0  1  22 
0123 
LI  2  4 
0125 
0  126 
0  12'’ 
0  128 
0  129 
0130 
0  131 

0  1  52 
0135 
0134 
0135 
0136 
U  1  37 
0  13B 
U  1 3V 
0  140 
H  1  4  1 
1)142 
0143 
0  144 
014S 
i  1  46 
0147 
0  148 
0  149 
1)150 
0151 
0  152 
0  15  3 
0154 
01S5 
0  156 
0157 
0  158 
0159 
0160 
0161 
0  162 
0163 
0164 
0165 
0166 
0167 
0  1  68 


AF2=AF2+F?( I ) 

S1)2  =  5L'?  +  F2<  I  )  **2 
40  CONTINUE 

CALL  LF1T<F2,L,A02,A12> 

AF2=AF2/L 

SD2=SD2/L 

SD2=DABS ( SD2-AF  2**2 ) 

«  SD2=DSQNT (SD2) 

SD3 -  0 . 

DO  50  1=1  ,.T 
APM-Af'H+PH<  I  ) 

LD  5=  SD3+PH ( I )**2 
50  CONI  INUE 
APH=APH/ J 
SD3  =  SI>3/J 

SD3=DA8S<SD3-APH**2> 

SD3=DSQRT(SD3> 

A  O  P  H  -  0  . 

SD4=U , 

DO  28  1=1,. I 
CPH ( I )=PH( I ) 

IF ( APH . L  T .90 .  . AND  .PH( I > .GT .270 . )OPH( I )=PH(I >-360  . 
IF  (  APH  .  01 . 27  0  .  .  AND  .  PHC I )  .  LT  .  90  .  >C.PH <  I  >  =PH ( I )  +360  . 
ACPH=ACPH+CPH(  I  ) 

SD4=SD4+CPH< I )**? 

7  1F(PH(1  )  .OF  . 0)UU  10  2D 
PH ( I  )  =P H ( I ) +360 . 

OO  TO  7 
28  CONI INUE 

CALL  l F 1 T ( CPH , J ,A03,A13> 

ACPH=ACPH/J 

ED4=SD4/J 

SD4  =  I)ABS(SD4-ACPH**2) 

SD4=DSQRT ( SD4 ) 

F 1 MX  =  0 . 

F2MX=0 . 

FMX=0 . 

PHhX=0 . 

DO  59  1=1,J 

IF (PH( I )  . GF .  PFlrt  OPHrtX=PH( I ) 

59  CONI INUE 
rtXPH=300 
F rtAX  =  0  . 

IF  (PHrtX  .GT  .45,  )CO  It)  55 
PHrlX-45. 

GO  TO  86 

55  IF  (  PHhX  .  GT  .  90  .  )GO  I  O  Sfl 
P  HhX  =9  0 . 

GO  TD  56 

58  IK  ( PHrtX  .  GT . 180 . >00  1  O  57 
PHrtX  180 . 

GO  IO  56 

57  IF  ( PHrtX  .  GT  .  360  )  GO  10  54 
PHHX =360 . 

CO  in  56 
54  rtXPH-=MXPH  +  i  00 
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Program  Listing  A-7-1 


PHSE  (cont.) 


PAGE  0u04  PhSE  (FTN4 — RELEASE  24177B — JULY,  1971) 


ul69 
0170 
0171 
0172 
0  1  73 
0  l  7  4 
0  1 7S 
0  1  7ti 
1)177 
0  1  78 
0179 
0180 
0  1 B 1 
U  1&2 
'1181 
U  1  84 
0  1  85 
u  1  86 
0  1  8*^ 
U  1  88 
0  187 
0  190 
0  191 
0  192 
0193 
U  1  94 
019S 
0  196 
0197 
0  1  98 
0  1<?9 
L?U  0 
02  01 
02  02 
02  0  3 
020  4 
0205 
0  20  6 
0207 
0  2  0  8 
0  2  0  V 
0  210 
0  711 
0  7 1 8 
0213 
0214 
0215 
0  2  lo 
021  7 
0  7  18 
0219 
J  2  8  0 
0  221 
0  277 
077  1 
0  784 


IF(PHMX.GT,MXPH)GO  TO  54 
PHMX=MXPH 
56  DO  60  I  -1  ,K 

IF (FI (  I  )  .GE.F 1MX  )F 1MX=F1 ( I ) 

60  CONTINUE 

Dl)  61  1-1, L 

I F  (  F  7  (  I  )  .  GF. .  F2MX ) F 2MX=F2 <  I) 

61  CONTINUE 

IF(F2HX . GE.F 1  MX) GO  TO  6? 

FMX  =  F 1 h  X 
GO  TO  66 

62  FMX=F2MX 

66  FMAX=F'MAX  +  i  0 , 

IF (F M< . GT . F MAX ) GO*  1  O  66 
F  MX  =  F  MAX 

6  i  DO  64  1  =  1,100 

P H ( I )  =  PH( 1  ) *F  MX/PHMX 
4  CONTINUE 

Nx-hAXO ( J , K , L ) 

DO  65  1=1,100 
X  (  I  )  —  I 

Y  (  I  ) =F 1 (  I ) 

Y  (  I  -•  1  0  0  )  =F  7  (  I  ) 

Yi 1+2*1  0  0)=PH( J ) 

65  CONTINUE 

104  FORMAT  <  3A2) 

185=41378 

CALL  Pl.OT  ( X  ,  Y  ,F  MX,  LABEL,  10  0,3,1) 
tURITKl  ,  110  )  1E5C<  1  )  ,  IFSC,  (9)  ,  I ESC (7) 

1  1 1)  FORMAT  <  3A2  ,  "  "  ) 

UR  IT  E (  1  , 1 05) (l  AUF< I ) , 1  =  1  ,36) 

105  FORMA i ( 1 6X , 36 A2) 

UR  T  TE ( 1  ,  101  )  I  CHI  ,  F  MX , Al  1  ,  IBS , SD1 
101  FUR  MAT ( 16X, "CM.  NO . " , 1 1 ,  "  xxxx,  GPH  MAX.  =",FS.l, 

1  "  CYL . /MIN . " ,7X  ,  "AUG.  FREQ.  =  ", F  6 A2 , F  4 . 7 , 

2  “  CYC . /MIN .  "  ) 

UR  IT  F ( 1  , 1 03 ) TCH2  ,  FMX , AF  2  ,  IBS  , SD2 
11)3  F  ORMAT  <  16X  ,  "CH  .  NO  .",11,"  oooo,  GPH  MAX.  =",F5.1, 

1  "  CYG  ./MIN.  "  ,7X,‘‘AVG.  FREQ.  =", F6 A7 , F 4 .? , 

2  "  C  TC . /MIN . " ) 

107  FORMAT  (  16X,  "PH.  ANG  .  (  "  ,  1 1  ,  "  urt  "  ,  1 1  ,  "  >  +  +  *+,  CI'H", 

1  "  MAX.  = " , F  5 . 1 ,  "  DEG . " , 3X , "AUG .  ANGLE  = " , F  6 . 1 , "  +  " , 

2  A2.F4. 1  ,  "  DEG. "  ) 

UR  I  IF  (  1  ,  102)  ICH7 , 1  CHI  ,PHMX  ,  APH,  IBM  ,  51)3 
URITEd  .  106)ACPH,IBS,SI)4 

10m  F  ORMAT ( 54X , "CORRECTED  AUG .  ANGIE  =  "  ,  F6 . 1 , "  +  " , A7 , F4 . i 
1  *•  DIG.") 

UR  IT  F  (1  .107)  IC.Hl  ,A()1  ,A11 

107  FORMAT ( 1 6X ,"  f  "  ,  II,"  =  " , F 8 . 3 , "  +  " , F  8 . 3 , " * X  "  ) 

UF  Tlf  ( 1  , 1 08) ICH7  ,AU2,A12 

10!)  FORMAT  (  16X  ,  "f"  ,  1  1  ,  "  =  "  ,  F  8 , 3  ,  "  +  "  ,  F  8 . 3  ,  H*X  "  ) 

Uk  I  TL  (1,1  09)A0  3,A13 

1(1°  FORMAT  (  16X  ,  "PH  =  "  ,  F 8 . 3  ,  "  +  "  , F  FT  .  3  ,  " * X "  ) 

UR  I  1 L  (1  ,  1  10  )  I  r  r  C  ( 1)  ,  IF  SC  (  9  )  ,  IF.  SC  (7  ) 

UR  I U  ( 1  , 1 04 )  I  fc  SC  <  1  )  ,  I E:  SC  (  5 )  ,  I  ESC  (  6  ) 

STOP 


0  0  0  1 
OGQ2 
0  0  03 
0  U  04 
0  0  05 
0  0  06 
D0  07 
0  0  0  ft 
0  009 
00  1  0 
0  0  1  i 
0012 
0013 
0014 
0015 
OGlti 
0017 


Program  Listing  A-7-2  LFIT 


PAGE  OCiOl  ( 0 TN4  —  RFLEASt  24177B —  JULY  ,  1971) 


FT  N4  ,  L 

SUBkCJtJT  I NE  L  F  i  I  <X,NUM,A0  ,A1> 
DIMENSION  X(l) 

SX=0 . • 

SN=  0  , 

SNS=  0  . 

S  X  N  =  0  . 

DO  10  1  =  1  ,  NIJM 
5X  =  SX  +  X  < 1 ) 

SN=SN+I 

SNS=SNS+I**2 
SXN=SXN+X( 1 )* l 
10  FONT  1 NUF 

A0  =  (f)X*UNS-fiN*‘iXN)/(  Nl)M*SNS-SN**2) 
A 1  —  <  NUH*SXN-SN*SX)/<NUM*SNS-SN**2> 
RFTUKN 
END 


**  NO  ERRORS* 


210 


APPENDIX  A-8 

PROGRAM  XCORR  (Cross  CORRel ation) 

Descri pti on 

Program  XCORR  computes  the  cross-correl ation  of  two  data  channels 
stored  on  disc.  If  the  two  specified  channels  are  the  same  an  auto¬ 
correlation  is  obtained.  The  normalized  correlation  versus  delay 
time  is  plotted  on  the  teleprinter  using  subroutine  PLOT.  Also  if 
the  correlation  function  is  cyclic,  its  frequency  and  relative  phase 
shift  are  determined  from  the  function's  zero  crossings.  The  flow  chart 
of  XCORR  is  given  in  Figure  A-8- 1  and  a  typical  output  is  shown  in 
Figure  A-8-2.  The  program  is  given  in  Program . Li sti ng  A-8- 1 . 

Required  input  data: 

1.  The  numbers  of  the  channels  to  be  processed.  If  the  two  channel 
numbers  are  the  same,  the  autocorrelation  of  that  channel  is  produced. 

2.  All  remaining  data  is  obtained  from  the  header  record  on  disc. 

Subroutines  required: 

1 .  DREA 


2.  PLOT 


211 


Figure  A-8-1  Flow  chart  of  program  XCORR 
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Figure  A-8-2  Typical  output  of  program  XCORR. 
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Program  Listing  A-8-1 


XCORR 


171558 


pa  :l 

Oi)  Gl  (F1N4  —  RELEASE  24177b  —  JULY, 

1971 ) 

■  051  1 

FTN4  i 

t  0  J2 

PR  OLKA."  XCORR 

f;  o  i)  3 

DIMENSION  IX(4,'24>  ,IY(4224)  ,f.(257)  ,X(257)  ,LAbl: 

1  ( 25 )  , 

1)0  04 

1 

LAbP  ( 135)  ,  IF:SC(9) 

GCOS 

EQIJIVAI  FNCE  (  Y  ,  1  Y  )  ,  (X  ,  IX  ) 

H  C  Oe, 

DATA  1 E  SC/ 1  5446  b  ,  ?Hk2  ,  ?HS  ,  2H  l  R  ,  2H  1  1 , 2HU  ,2HI) 

,  2Hk  0 ,2H16/ 

0  G  >5 7 

DATA  l  AbFI  /?*  1  H  ,  1HC  ,  1  HO  , 1  HR , 1  HR  ,  IMF  ,  1  HI ,  1  HA  ,  i HT  ,  1 H I  ,  1HO  , 

1  0  0  H 

1 

1HN . 1 H  , 1HM , 1HA , 1HG , 1HN , 1H1 ,  1H T  ,  1HU  ,  1HI) ,  i  EIF  / 

0  .j  0  9 

r.Al  1  EXE  C( 17, IF) RK , ILTRK ,  I  SIZE) 

0  0  1  0 

CAI  l  erf:  A(  IF  IKK  ,0,1  A  HE  ,1,1,1) 

0  G  1  1 

Sk  -l  AHf  <  3B  )  /60 

Gut;-' 

NKlLB-  1  AF  P  (  39  ) 

g  n  i  3 

NCllA-l  AHF  (37) 

l)  1  1  4 

DO  15  1-1,4224 

■If  15 

I  X  (  I  )  --  0 

0  1  1  *J 

1  Y  (  l  )  =  0 

OL  17 

15 

lIIN  1  INK! 

U  0  1  fl 

DO  19  1-1,257 

0019 

0  (  I  )  -  0  . 

■  ■0  "0 

19 

CONI  I  NHL 

ohm 

1 

UK  I  TE(  1  ,  1  011  ) 

H  0  22 

1 014 

FORMAT (/"INPUI  THE  KF  QtJiRFD  CHANNt.l  5?  “  ) 

0  0  23 

READ  (  1  ,  *  )  10H1  ,  10H? 

0  0  24 

IF ( 1CH1 . GT . NEHA ) GO  1 0  1 

0  0  25 

IF  (  T  PH?  .  ITT  .  NEHA  ) GO  Tl)  1 

l)  026 

WRITE  (1  ,10  1) 

0  0  27 

1  0  1 

FORMAT < /“INPUT  MAXIMUM  TIME  SHTFT  IN  SECONDS? 

" ) 

1)0  2d 

K  1  AI)  (  1  ,  *  )  T  MAX 

0  0  79 

N - I F I X  < TMAX#SR/2  .  >*2 

On  3  0 

IF  (N.l  E.  12E1  )G(i  TO  2 

0031 

N-- 1  28 

0  0  72 

TMAX-N  '  5R 

0  0  3  3 

UK  IT  F  ( 1  ,  1  02  )  rnAX 

U  0  34 

l  02 

r  UR  MAT  (/"MAXIMUM  TIME  SHIFT  REDULLI)  TO  "  ,  F  7 , 1  , 

“  SECONDS") 

0  0  35 

l.  NET  H-4  u  96 

0  0  3* 

IS  i  -1  29 

0  03  7 

1ST  - 1 29 

0  0  3K 

nKF  m  -  nf  1) i  NRFE!> ,  Ik ) 

0  0  3? 

N  I  Inir-(  NKF.CS-NREM)/16 

0  0  4  0 

If  (  N  Fl  ME  Nf  .  0  )  00  IT)  3 

0041 

4 

l  NL,  F Ef-Nkl. M*2L  6 

0  042 

N  F  T  M  c  -  1 

0  043 

Nf<  EM  •  H 

0  0  44 

3 

DH  2  0  .1-1  ,  N  (  I  Ml 

0045 

IREC  - 1  J  1 ) *  1 6  + 1 

004fe 

CAI  1  I)  <1  A(  TFTKK  ,1  NGTH,  IX  ( 1ST  )  ,  TREE  ,  IEFI1 ,1 ) 

5  0  47 

f  ALi  D K .  A (  IF  TKK  ,  LNG  1  H  ,  I Y  (  TSY  >  ,  I KL  E  ,  J  EH 2 , 1  ) 

0043 

DO  1i)  l  =-N  N 

C  049 

DO  it  K  1 . 1  NkT  H  M  23 

0  0  5  0 

II v=T+K 

0  051 

IF  tJlx.it  .  u  )  Cl)  1  i)  16 

0052 

I M A  *  - 1  f*  0  H  +  1 2  B 

0  n  S3 

IF  (  l  IX  .  r,  I  .  IMAX  >GO  Til  1 6 

0  054 

C(N  +  1  +  t)=C(N-»  i  +  I  )+EI  OAT  (IY(K)  )*FLCiAI  (  IX  (UK)  ) 

0  0  55 

16 

CONI  TNIJE 

0056 

10 

CONT  I  NIIE 

MTIft 
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Program  Listing  A-8-1  XCORR  (cont.) 


•'AuE  0  0  02  XCORR  (FTN4 —  RELEASE  ?.A  1 77B  —  Jt.IL  Y  ,  1971) 


..OS'7 

DO  ? 0  K  = 1 , 128 

0  Jc  6 

IX(K  >  = l X ( L NGTH+K ) 

(  i  :  c,9 

l Y (K  )-  1  YILNGTIH  l<  ) 

i  (]  y>  >| 

?li 

l  ON 1  INUE 

0  061 

It  (NR1  M  .  NF  .  0  )GO  TO  4 

10  6.? 

N=2»N 

0  0/1 

6 

Y  MAX  =  0 

0  0  64 

DO  13  1  =  1  .  N-t  1 

0  0  f,s 

AY  =ABS ( C (  I ) ) 

0  :J .  >  6 

IMAY.GI.I  MAX)  1  M  A  X  =  1 

0  Of  7 

IF  ( AY . G 1  .  YMAX ) YMAX  =  AY 

0068 

1  3 

COM  I  MJE 

0  0  69 

DO  14  1=1 ,N+1 

0  0  70 

l.  v  l  ) =  C  i  1  )  /  YMAX 

0  0  7  1 

X  (  [  )  =  (  l  -N/2  .  -1  )  /SR 

0  0  7? 

i  1 

CON  T 1  NOE 

0071 

NP  Fk  =  0 

0  0  74 

HE  R--0  . 

0  0  7S 

I>1  R  -  1  . 

0076 

Du  21  1-1  N 

0  0  77 

F.N1  =SJCN<  1  .  ,  C(  1  )  ) 

0078 

6N,*--SI(,Nt  1  .  ,l:(  1  «  1)  ) 

0  079 

1F(SN1 .EQ.SN3)6Q  IO  21 

0080 

1KC.  (  l  +  l  )  .1  1  .(.(1  )  )  D 1  k  =•  •  - 1  . 

0  UE  1 

1  1=  (('.(  1  +  1 )/  UK  I  +1  )-C<  I  >  >  )*DIR/RR 

0082 

r,D  1  (  22 

0  u  63 

21 

COM  I  N1IE 

0  0  84 

0  -J 
&  <L 

Oil  23  J  =  1  +- 1  N 

0  0  8  s 

SM  =  S  1  ON  (  1  .  1.  (1)) 

0  0  6  o 

SN2- S I  UN ( 1  .  ,C(  1  *  1  )  > 

0  JH7 

11  <  6N1  L  ij.nN?)Gu  TO  23 

0  0  46 

U  <  SNE  .  NF  .  I)  1R  )l,0  TO  23 

0089 

t?=uk  j  +  i  >/<r<  .1  + 1  )-c(j  ) )  )*dir/sk 

0090 

PERM  1  +  (  J-I ) /SR- 1 8+PF fi 

0091 

NPF.R  =NPE  R  + 1 

0  092 

GO  TO  21 

0  09  $ 

?  3 

CONI  1 NLE 

0094 

GO  TO  ?S 

0  0  95 

24 

1= J 

0  096 

T1  =1  2 

0097 

LO  IO  22 

0  0  98 

?* 

F - ( FLOAT ( NPE R ) /PER  >*6Q . 

0099 

1 MX=( FLOAT ( I  MAX )-N/2 .-1 )/SR 

0100 

P FIA'8 - 3 •> 0  .  *1  MX*NPER/PLR 

0101 

RATI  0  =  1  MAX*  1 1)  0  .  /  (  2S6*I  ABF  (39  >*6  0  .  /I..ADF  (36)  ) 

0102 

CAM  Pi  O T  (  X  ,  f. ,  0  .  ,  LABE  L  ,N+1  ,  1  , 1 ) 

0  1C  3 

UR  I  11  (  1  ,104) 

0104 

104 

FORMAT  <  42X , "1 I  Ml  ( SECONDS ) " ) 

0  10S 

URIlE-t  1  ,10S5  I  CHI  ,  ICH2,  ( L  ABF  (  I  )  ,  1  =  1  ,36) 

0  106 

10S 

FORMAT ( /20X ,  “CH .  “  ,  1 1  ,  "  RELATIVE  TO  OH.  "  ,  1 1 

0107 

UR  11 E  ( 1  , 107) RAT  IO  F.PHAS 

U  1  06 

1  0  7 

FORMAT  (  2  (IX  ,  SHIFT  3X,"ERTU.  =M,Fb. 

0109 

1 

3X  ,  "PH ASF  =  "  , F  7.2,  "  DF.  GREFS"  ) 

0110 

UR  11  F  <  1  ,  1  0Q ) YMAX 

3111 

1  UR 

F ORF A1 ( 2UX , "MAX  1  MUM  CORRI  1  AT  1  ON  AMPLITUDE  IS 

0112 

UR  UK  (  1  ,103)1  r!it;  (  1  )  ,  ItSC(S)  ,  IESC(  6) 

CYCLLS/M  IN 


"  ,  E 1  2  .  6  > 


' 


Program  Listing  A-8-1  XCORR  (cont.) 


0  113 
0  114 
011S 


*4  NO 


Pr-Gl  U  003  XCUKK  (FI  N4  —  R  EL-E  AST  E4177B — JULY 


103  FfJkMAT  1 3A?  ) 
STOP 
END 


ERRORS* 


,  1971) 
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APPENDIX  A-9 

PROGRAM  MAIN 

Program  MAIN  is  a  master  program  which  permits  the  researcher  to 
process  up  to  eight  data  files  using  any  or  all  of  the  previously  des¬ 
cribed  programs.  The  input  files  are  read  from  magnetic  tape  and  are 
specified  on  the  basis  of  their  file  numbers.  The  file  length,  number 
of  records  to  be  skipped  at  the  beginning  of  each  file,  and  the  sample 
data  reduction  rates  are  input  for  all  eight  files.  The  programs  can 
be  applied  twice  to  the  same  set  of  data  files.  Program  MAIN  provides 
a  list  of  the  available  programs  and  requests  the  order  in  which  these 
programs  are  to  be  applied  to  the  data  in  each  of  two  processing 
applications.  An  input  order  of  zero  indicates  that  a  particular  program 
is  not  to  be  applied.  Program  MAIN  outputs  each  program  name  and 
awaits  two  sequence  numbers,  one  for  each  processing  pass  through  the 
data  files.  The  following  is  a  typical  output  followed  by  the  required 
sequence  numbers: 


PROGRAM: 

1  st . 

2nd 

FLOCT 

1, 

1 

DCAC 

0, 

2 

FILTR 

0, 

3 

PSPT 

2, 

0 

PHSE 

0, 

5 

XCORR 

0, 

4 

The  above  set  of  sequence  numbers  will  result  in  programs  FLOCT 
followed  by  PSPT  being  applied  to  all  specified  data  files  (a  maximum 


■ 


. 
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of  8  files)  during  the  first  processing  pass.  During  the  second  pass 
programs  FLOCT,  DCAC,  FILTR,  XCORR,  and  PHSE  will  be  applied  to  the 
data  in  the  order  listed.  Program  FLOCT  will  always  be  first  in  the 
sequence  unless  only  one  file  is  to  be  processed  and  it  already  exists 
on  disc. 

Required  input  parameters: 

1.  Number  of  files  (a  maximum  of  8). 

2.  Specific  file  numbers. 

3.  File  lengths  (number  of  256-word  records). 

4.  Skip  sequence  (number  of  records  to  be  skipped  at  the  beginning 
of  each  file). 

5.  Two  sample  reduction  rates  (one  for  each  pass). 

6.  Input  sample  rate. 

7.  Number  of  channels. 

8.  Number  of  segments  for  fourier  analysis. 

9.  Maximum  time  delay  for  cross-correlation. 

10.  Digital  filter  cutoff  frequencies. 

11.  Two  sequence  numbers  for  each  program. 

Programs  and  subroutines  required  are: 


1. 

FLOCT 

9. 

STPT 

2. 

DCAC 

10. 

RDPT 

3. 

FILTR 

11  . 

PLOT 

4. 

PSPT 

12. 

I FFT 

5. 

PHSE 

13. 

SMTH 

6 . 

XCORR 

14. 

BPASS 

7. 

DREA 

15. 

LFIT 

8. 

DWRI 

. 


Program  MAIN  is  given  in  Program  Listing  A-9-1 


■ 

' 


m rm 
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0  0  0  1 
UU02 
0  0  03 
0  0  0  4 
0  0  0b 
0  0  06 
0  0  07 
00  08 
0  0  09 

n  o  10 
001 1 
00  12 
i)  U  1  ! 

II IJ  1  4 
001b 
0016 
0  0  17 
0018 
0  0  19 
0  0  20 
0  0  21 
0  U  2? 
0  U  2  i 
u  074 
0  0  2b 
00  6 
0  0  77 
0  0  7  S 
0  0  79 
0  0  30 
0031 
0  0  32 
00  13 
0  0  14 
0  0  3b 
0  0  1 6 
0  0  J7 
0  0  lb 
u  0  39 
0  0  4  0 
0  0  4  1 
0  0  42 
l  04  1 
f  0  4  4 
0  0  4b 
0i46 
0  0  4  7 
0  04H 
0  0  4v 

on  so 
o  0  b  1 
00S2 
0  0  b  3 
0  G  b  4 
OObb 

171737  uos6 


Program  Listing  A- 9-1  MAIN 


PAGE  0  0  0  1  1 i- TN4---RFLEASE  241778  —  JULY  ,  1971) 


FI N4  ,  L 

PROGRAM  MAIN 

DTMtNSION  I SLQ  (  6  >2)  ,JF  II  P.  <  B )  ,NS<2)  ,L.0CT(3)  ,  IDCAC(v->)  ,  IF1LTC3)  , 

1  1PSPT  (3)  ,  IPHSt.<3>  ,  IXCOE  (3) 

I.OMMdN  l  AHF  <  13b)  ,  IF  11  F . ,  1CH1  >  ICM2,NRC(8>  ,NSK (8)  ,  NIH  ,  NBEO  ,  FL  ,PH 
DATA  l  0CT/2HI  0,7HCT  ,?HS  /  ,  TDCAC/PHDG  ,  7HA0 , 2Hb  /  ,  IF  II  T/2HF  I , 

1  7HI  1  , 2Hb  / , IPHSE/?HPH,2HSE ,2H5  /  ,  IXC0R/2HXC , 2HOR , 2Hb  /, 

1  I  POP  I /2HPS , 2HPT  ,  2H5  / 

UR  I  TP  (  1  ,10  0) 

100  F  l  IK  MAT  (  /  "NO  .  OP  PILES?  ) 

K  h  A  D  (  1  ,  *  )  N  F  1 1  EE* 

ur  nr  (i  ii3) 

II  1  I  III  MAI  "1  Ill'll  I  llir  Ml  1(4)  (I  KG!  I  "> 

I  I  i,l  M  mil  III  Il),|  1  , HI  11  IS) 

UR l  I  I  <1  ,10  3) 

101  FORMAT  (  /  "  INP I )  I  RECORD  LENGTHS !  .") 

R l  AD (  1  ,  #  )  ( NRO ( I )  ,  1  =  1, N  FILES) 

UK  1  I  P ( 1  ,102) 

102  I  (IKMA  I  (  /  “RF  CORD  SKIP  SEQUENCE !  !  ") 

Rl  AIM  1  ,  *  )  ( NSK  (  J  )  ,  l-l  ,  NF  ILFS) 

UR  1 1 P  (  1  ,106) 

1J6  FORMAT (/"INPUl  1W0  SAMPLE.  REDUC1  ION  RATES!!  "  ) 
k  L AD ( 1  ,  *  )  NS ( 1 )  , NS ( 2 ) 

UR  11  I  (1,114) 

114  FORMA r< /"INPUT  SAMPLE  RATE!!  „") 

RL AD( 1  ,  *  )SR 
UR  f  TE  l  1 , 1  lb  ) 

lib  FORMAT  (/ "INPUl  NO.  OF  CHANNELS!!  “> 
t  READ ( 1 , * ) NO 
UR  I  I  P  ( 1  , 121  ) 

171  P  OR  MA I  (/"INPUT  NO.  OF  SEGMENTS !  !  „") 

R  L  AD ( 1  ,  *  >  NSG 
UKJTL(1 , 122) 

1  72  IORMA1  ( /"INPUl  MAXIMUM  1IME  DP  LAY  '  !  ") 

Rl  AD ( 1  , *  )  T MAX 
UR  (1  F  <  1  ,  1  04  ) 

104  FORMA1  </“  I  Mi'll  1  1  WO  SEQUENCE  NUMHPRF.  OPPOSTTP  F  ACM  PROGRAM'!"///: 

UR  I  I  F'  (  1  ,10S)  1S446H  ,  2H(1D,  1  S4468,2H(IA 
1  US  F  i IRMA F ( 2A2 , "PROGRAM :  "  ,  6X  ,  "1ST .  ,?ND  .  "  ,  PA?// ) 

UR  11 1  (1  ,107) 

107  F  Ok r)AT  (  1>,  ,  "FLCJOT"  ,  1 1 X  ,  "  ") 

RLADL1  *)ISEQ(1,1),  1SEQ  (1,2) 

UR  ITE (  1  ,108) 

1  i  a  f  OR  MAI  (  1  X  ,  "Dl-AC"  ,  12X  ,  "  ") 

Ri  A  I)  ( 1  ,*)  I SEU  (  2 , 1  )  ,  I  SL'Q  (  2 , 2  ) 

UR  I  TF  (  l  ,109) 

t  09  F  OKMA1  (  IX  ,  "HI  IK  "  ,  1 IX,  "  "  ) 

READ  (  I  :*)1SEU(.4,  1 )  ,  1  SE.Q  (3,2) 

UR  F  TC ( 1  ,110) 

110  FORMAT ( IX , "PSP  I " ,  1?X  ,  "  "  ) 

REAIH  1  ,  *  )  I  SEIJI  4 , 1  )  ,  I SEQ  (  4,2) 

OR ITE (  1  ,111  ) 

III  f  ORMAT  (  IX  ,  "PFISF  "  ,  12X,  "  ") 

R  FAD  ( 1  .  *  )  1  Sb  O  ( b  ,  1 )  ,  ISE  Q  ( b ,  2 ) 

UR  II P  ( 1  ,112) 

11?  P  Ok MA  1(1*.,  "XLORR  "  ,  i  IX  ,  "  ) 


uim 
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Program  Listing  A-9-1 


MAIN  (cont.) 


PAGE  0  0  0 MAIN  (FTN4 —  RELFASF  241  77B— JULY,  1971) 


0057 
GO  SB 
0  0  59 
0  0  60 
0  061 
0062 
0  063 
U  0  64 
0  065 
0066 
0  0  67 
0068 
0069 
U  0  7  0 
0  0  71 
0  .i  72 
0  0  73 
U  074 
0  0  75 

0  07ft 

0  077 
0  0  78 
0  079 
0  OHO 
0081 
G  0  82 
0  0  83 
0  0  84 
0  0  85 
0  086 
0  0  87 
0  0  88 
0  089 
0  0  9  0 
0  0  91 
0  0  92 
0  09  3 
U  0  v  4 
0 .1  v  5 
0  0  96 
0  U97 
0  098 
0  0  9  V 
0  1  II  0 

o  i  o  i 
0  10  2 
0  1  111 
1104 
0  105 
0  106 
0  1  0"? 
0  1  l)  8 
0  1  09 
0110 
0  111 

1 7 1 7  3  8 n ' 1 2 


RF  A1X  1  ,  »  )IBEQ(6 , 1  )  ,  1 SEQ  (6,2) 

UR  n  F<  1  ,  125)  15446B  ,2811 ,2H0 
IF  II.  I  -0 

D(l  10  l  =1  ,  NF  II  FS 
DO  ?IJ  J-l  ,2 
N  TH-M6 ( J  ) 

FL  =  1  . 

FH  =  SR*  .67/(2  .  *NTII) 

I  '8  20  K  =  1 , 6 
DO  20  1=1,6 
15=15 FQ(L,J) 

I F  (  1 8  .  N  E  .  K  )  G  f.l  10  2  0 
GO  1D( 1  ,2,3 ,4,5 ,6)  ,L 

1  I  FI  l  t  =  JFll.F .(  I  ) -IF  ILL 
ir:m=i 

NSEG  = 1  FIX ( SR > 
f  ai  l  fx eg (8, inn  ) 

NCHAX  A HF  (  37  ) 

1  I  (  Ml  I  I  .  NGHA  >NC;HA=NC 
IF  II  F  =  .IF  II  F  (  1  ) 

GO  II)  20 

2  GAI  L  F  XLC ( H  ,  1  DCAC) 

OR  I TE< 1  ,118) 

118  FORMAT  ( /15X  ,  "DCAC  GOMPL.fc.TED "  ) 

GO  10  20 

3  CAl  l  F  Xt7:<8  ,  ]  Fit  T  ) 

UR  I  n  <  1  ,  1 3  V)FL  ,F  H 

119  F  OR  MAI  ( / 1 5X  ,  "F  1 1  Ik  COMPUTED  WITH  BAND  PARS  "  }  fc  1  2.3,  "  10  ", 

1  F  12 .3  ,  "  CrCI.FR/M  IN  .  " ) 

GO  10  2(1 

4  NSEG--NSG 

UR  HE  (  1  ,120) 

120  FORMAK/15X  ,  "PSP  r  INITIATE  I)  !  !  "  ) 

FH=LALF <  38 )*  . 67/2 . 

1 SEG=LAKF  (  39)/NSEG 

SFGI  =(60. /I  AHF(.i8)  )*856*iSFG 

l.Rl  OSF.Gl  *N5F  G 

LSI  G  =  RFGL 

nSFG  =  SEGI  /6 0  . 

MRF  C  =  l  KEC/60 . 

TSfcr:R=MOD(l  REG, 60) 

1SI  f;S  =  MOD(l.  SI  G  ,  60  ) 

UR  I T  F  <  1  .  123 ) MR  EC , 1SLCR 

1  23  FORMAT  (  / 1  5X  ,  "RF-LOR  D  IFNG1H  =  "  ,  X  ,  1  1 ,  X  ,  "M  INlITtS  "  ,  X  ,  1 3 ,  X  ,  "RFC  .  "  ) 
UR  11  F  (  1  ,  124  ) h S E fl ,  ISF  CS 

124  FORMA!  ( /15X  ,  "SFGMENI  LFNGIH  =  "  ,  X  ,  13  ,  X  ,  "M 1  NOT  fc.8  "  ,  X  ,  1 3  ,  X  ,  "BF.C  .  "  ) 
NCHA  =  L  AbF ( 37 ) 

1F(NC.I  1  ,NCHA)NLHA  =  NC  . 

f  DO  21  N=i , NCHA 
ICHi=N 

I  AI  L  F>  fc.C  <  8  ,  IPSPT  ) 

uR I  I  E (  1  ,  125)  1 54  4  6b, 2H11 ,2HV 

125  FORMAT (3A2) 

21  LONTINUF 

GO  III  2 U 

5  UR  I1E(  1  ,  116) 


Program  Listing  A-9-1 


MAIN  (cont.) 


PALL  OIJiU  MAIN  (FI  N4 —  RELEASE  2 


0113 

1  1  8 

KIKP1AI  (  /1SX  ,  "PHSE  INITIATED! !“ ) 

0114 

DO  2?  N  =  1 , NCH A 

011S 

ICHl - N 

Oil  8 

iCH  >-Hi 1 

Oi  17 

IF  (N.f.Q  .  NCHA  )  1CH2=1 

01  18 

CALI.  F  xt  C  (  8  ,  I  PHSF  ) 

0119 

Uk  11 F ( 1  , 12S) 15448b ,2H1 1 ,  2HV 

0120 

O') 

CON  HNIJE 

0  l  p  1 

MJ  10  2  0 

0  12? 

i1-) 

1  EMP  =1  L 

0  123 

F  L  =T MAX 

0124 

WRIIECi , 1 17) 

0  l  25 

117 

F  OK  MAI  </15X,  "XCfJKR  INITIATED!  !  “) 

0126 

N  =  1 F  JX(TMAX*l  AHF  <  38  1/120  .  >*2 

0127 

I  F  (  N  .  0 T  .  1  28  )T  MAX -128  .  *80  ,  /I.  AHF  (  3B) 

0128 

UR  n  F( 1  ,128) TMAX 

0129 

128 

F  OkMAT  < /1SX  ,  "MAXIMUM  TIME!  SHIFT  IS 

0  130 

DO  23  N- 1 , NCHA 

0131 

I  C  H  1  =  N 

0132 

I CH2-N+ 1 

li  i  13 

IF  (N.f  IJ .  NCHA)  ]CH?  =  i 

0134 

CALI  E XEC(8 , IXCOR ) 

0  1  3S 

WR  I  1  F ( 1  ,  125) 15448b ,2H1)  ,2HV 

013m 

5  3 

CONI  1 MU  £ 

0137 

Fl  =TEhP 

0  t  38 

20 

CUN  I 1 NUF 

0  l  39 

l  0 

CONI  1  NHL 

f  1  40 

SI  UP 

0  14  1 

END 

*♦  NO  ERRuftS* 


4177b — JULY , 


"  >F7.1) 


1771  ) 
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APPENDIX  B 


The  modified  Krebs-Ringer  solution  used  in  this  work  had  the 
following  composition: 


Sodi urn 

144.0  mM/1 iter 

Potassi urn 

5.4  mM/1 iter 

Calcium 

2. 5  mM/1 i ter 

Magnesi urn 

1 . 2  mM/1 iter 

Chloride 

128.8  mM/1 i ter 

Bi carbonate 

22.0  mM/1 iter 

Phsophate 

1  . 2  mM/1 i ter 

G1 ucose 

10.1  mM/1 iter 

